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SUMMARY OF MATHEMATICAL MODELS FOR A CONVENTIONAL 
AND VERTICAL JUNCTION PHOTOCONVERTER 

By 

John H. Heinbockel* 

1. ONE-DIMENSIONAL CONVENTIONAL PHOTOCONVERTER 
The first model considered is for a one-dimensional analysis of a con- 
ventional photoconverter as described in References 1, 2. The geometry for 
this conventional n/p photoconverter is illustrated in Figure 1. We assume 
a generation rate 

g(x) = <|>o (1-Re) 0 exp (- a X) (1.1) 

and solve for the minority carriers in the n-region. These minority car- 
riers must satisfy the diffusion equation 

°p P no> - - n -'- Pn0> - -9(*> d.2) 

d X 2 Tp 

subject to the boundary conditions 

D — (P n - PI = S n (P n - P n J at X = 0 (1.3) 

p ,j x n no P n no 

and 

P - P =0 at X = X.. (1.4) 

n n 0 j 

The solutions of the equations (1.2) (1.3) (1.4) gives the photocurrent 
density at the junction edge. 

J p = -^p ^ < p n - p no) at x = x j t 1 - 5 ) 

♦Professor, Department of Mathematical Sciences, Old Dominion University, 
Norfolk, Virginia 23508. 






which evaluates to 


where 


q<l>o (1-Re) a L p 

-S 2 " 1 


( f l/ 2 - ctLp Exp (-aXj)) 


and 


fi 


S L 

- + S 

P 


-exp (-aXj) 




S L 

f 2 • f(i, _!£, x,/l) 
0 J p 


f(a, b, c) = a cosh (c) + bsinh(c) 


The diffusion equation for minority carriers in the p-region i 

by 


(n - n ) 

D Jl ( n - n ) - -E 521 > -g(x) 

P dx* P P° 


with boundary conditions 


n - n =0 at X = X .+w 
P PO J 


-°n d7 <n P " "pc’ * S " (n P ' "po 1 “ X ‘ H ' 

This produces the photocurrent density 


( 1 . 6 ) 

(1.7) 

(1.8) 

(1.9) 
given 

( 1 . 10 ) 

( 1 . 11 ) 

( 1 . 12 ) 


3 



at 


(1.13) 


J = q 0„ — (ri - n _) 
n n dx p p0 


X=X -+(d 

sj 


which evaluates as 




d-Re)^ n 


0 2L 2 
n 


exp (-a(Xj+o))) [oL n - f3/f«»] 


where 


S L S L H-X.-u 

: 3 = [ 0 L - JL1] exp [-a(H-X--io)] + f (— , 1, L_) 

D„ J D„ L„ 


S L H-X.-m 

U = f (d — , 

0 L 

n n 


where f is defined by equation (1.9). 

The photocurrent from the depletion region is given by 


-aX . -au> 
J 


J. = q<f> 0 (1-Re) e J [1-e ] 


and the total current due to the light source is 


J = J + J + J , 
sc p n dr 


For the dark currents we have the injection current J. .and the 

inj 

recombination current J which are given by 

rec 


(1.14) 


(1.15) 


(1.16) 


(1.17) 


(1.18) 
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(1.19) 


qV 

J inj = J o [«» ( F ) - 1 ] 

with 


Jo = 


qD n. 2 S L 

• r\ -j n 

r i 


X. / S L 

p -1- f C _LE . l, JL) /f(i, -±1 , J.) 

N_, k D L 7 0 L 

P 


X. 

_J_ 

‘P 


qO n. 2 S L H-X.-u. / S L 

+ _2. _ L f(JLH , l, _J_)/f(i, JUL 

L N 0 L / D 

nan n n 


H-Xj —to) (1.2U) 

J 


and 


J 


rec 


frqn^to 

Vv^ 


sinh fill 
2KT 

q(V bi -V)/kT 


(1.21) 


The dark current is J n - J. . + J and the total current is J = 

D inj rec 

J sc" J D* 

The effect of a series resistance R and shunt resistance R , can be 

s sh 

included in the model by considering the equivalent circuit diagram of 
Figure 2. 

Nominal values for the conventional solar cell parameters are given in 
the computer program "conven" listed in the Appendix A. In this computer 
program the absorption coefficient as a function of energy is obtained from 
interpolation of numerical input data. Empirical formulas approximating 
lifetimes, mobilities, intrinsic carrier density and bandgap energy as 
functions of temperature are included in the program. 
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2. ONE-D IMENS ION VERTICAL JUNCTION PHOTOCONVERTER MODEL 
The geometry for the vertical junction photoconverter is illustrated in 
Figure 3. We consider a narrow strip at a depth n below the surface y = A. 
Following References 3, 4, we treat this strip as a conventional photo- 
converter and then sun the results over all n, 0 < n < A. 

In the n-region we have the diffusion equation 

°P ^ ' P no> - = ' 9tA - y) ’ (2 ' 1) 

subject to the boundary conditions: 


P 

n 



= 0 


at X=X .+X 
J n 


( 2 . 2 ) 


•D n _ (P -P n J = S n (P -P n J 
p d x n no P n no 


at 


X = B 


(2.3) 


From this equation we obtain the current density at X ,+X given by 

J n 


J p • 0 °p - ( p n - p no) X**J + *n 


(2.4) 


which simplifies to 


o„ . . g(A - y , qL , 'VVWW 


P «V L P' S p> ( B - X j- X n) /L p ) 


(2.5) 


In the p-region we have the diffusion equation 
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fjURE 3 ' Vertical Junction photoconverte 



( 2 . 6 ) 


d 2 

D (n -n ) - (n -n )/t = -g(A-y) 

n dx 2 p po p p0 n 


which is subject to the boundary conditions: 


n -n = 0 at X=X.-X 
p po j p 


(2.7) 


D„ — - (n„-n„„) = S„ (n -n„ rt ) 
n dx p po n P P° 


at 


X=0 


( 2 . 8 ) 


This produces the current density at X.-X given by 

J P 


J = -q D — (n -n„) at 

n “ n ' n nrw 


dx 


P po' 


X=X 4-X 


J P 


(2.9) 


which simplifies to 


J n - -9(A-y) 



V f <V 

f < D n /L n- S n- < x j- X p )/L n) 


( 2 . 10 ) 


The total current is obtained by integrating these current densities over 
the depth of the photoconverter. Assuming the photoconverter is 1 cm in 
depth we have: 


A A 

I = / J (n) dn + / J (n) dn = I + I (2.11) 

bL g P 0 n P 


The dark current is given by I . = I * + I * where 

3 J d p n 
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Using P. = B Wo (1000) (watts) as the input power and ¥ = P m ,„ as the 

i n OUt ITiaX 

maximum power obtained from the I-V curve, we calculate the efficiency as 


« 100 P. IP . 
ff in out 


(2.14) 


Nominal values for the vertical junction photoconverter are given in the 
computer program "VJSCP" listed in Appendix B. We use the same empirical 
formulas approximating the lifetimes, mobilities, intrinsic carrier density 
and bandgap energy as a function of temperature as listed in the Appendix B. 
For both computer programs listed in Appendix A and Appendix B the tempera- 
ture of the photoconverter is assumed controlled by the use of a heat pipe. 
The current voltage relationship is obtained from the equivalent circuit 
diagram of Figure 2. 

3. THREE-DIMENSIONAL CONVENTIONAL PHOTOCONVERTER 

The geometry for the three-dimensional conventional photoconverter 

model is illustrated in Figure 4. We let R n = {(x,y,z) | x | < a , | y j <b , 

0<Z<Z .} denote the n-region and let R = {(x,y,z)l |x|<a, jy j <b, Z.+w<Z<H} 

J P I J 

denote the p-region. We assume that the light impinges upon the surfaces 
z=0. The equation for the diffusion of the minority carriers in the n- 
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X 




photoconverter 



material is given by 


L(P -P ) = V 2 (p -p ) - (p -p )/L 2 = -g(Z)/D , x,y,zeR (3.1) 
n no n no n no p p n 

where L() is the linear operator v 2 () - ( )/L 2 and g(x) is the generation 

P 

rate defined by (1.1). Boundary conditions for the above equation are: 


°p ^ < P n- P no> “ S p< Vno> at Z=0 • Vno* 0 at Z=Z j (3 ' 2 > 

"^p — ^ P n' P no^ “ ^ P n -P no^ at *” a * a7 ^ P n~ P no) = So (P n -P n0 ) 

at x=-a (3.3) 

" D p ~ ( P n" P no^ ~ ^°^ P n“ P no^ at y=t) ’ % ( p n _P no) = So ( p n _p no) 

at y=-b (3.4) 

where we have assumed that the surface recombination velocity is So on the 

sides and S on the surface z=0. 

P 

Similarly, the equation for the diffusion of minority carriers in the 
p-material is given by 


L(n -n ) = V 2 (n -n ) - (n -n )/L 2 = -g(z)/D , x,y,ze R (3.5) 

p po J p po' p po' n ' n’ p 


with the boundary conditions 


n p - n po = 0 at z=z j +w 



(n -n ) = S (n -n ) 
v p po' n v p po 


at Z=H 


(3.6) 


12 


D n W^W 1 

at x=-a (3.7) 

3 

n ( n p- n p 0 )-Si( n p- n po ) 

at y=-b (3.8) 

where Si, S p are respectively the surface recombination velocities of the 
sides and bottom surface z=H. 

The solution to the equations (3.1) and (3.5) which satisfy the 

corresponding boundary conditions is given in terms of Green's functions 

&i(x,y,z; x 0 ,y 0 ,z 0 ) for the region R n and G 2 (x,y,z; x 0 ,y 0 z 0 ) for th e region 

R . 

P 

The Lagrange identity associated with the linear operator L( ) is 

uL( v) -VL(u) = U V 2 V - W2 U . (3.9) 

By Green's second formula we can write 

///vol[“L(v)-VL(u)j dxo = ///volt uv 2 V-Vv 2 u] d-r o = // s ( uVV-Wu)«vdao (3.10) 

where d-ro is a volume element, dao is an element of surface area and v is a 
unit normal to the surface. Let V = Gj(x,y,z; x 0 ,y 0 ,z 0 ) be a Green's 
function and let U = P n -P no denote the solution to (3.1) with vol = R n , then 
the integral (3.10) can be written 
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/// R (P n -Pno> L < G i) d *o -111 G > ( 

n n 


-9(Z 0 ) 

D_ 


■) *0 


// s [(Pn-Pno) vS i- G i 7 ( p n- P no^‘ vdo ° (3.11) 

For L(Gi) = 6(xo-x) 6(yo-y) <5(zo-z), 6 = Dirac delta function, the above 

gives 


Z. 

J 


P "P = / 

n no o 


/ b / a 

-b -a 


-a<po ( 1 -Re) 



D 


P 


Gj (x,y,z; x 0 ,y 0 ,z 0 ) dx 0 dy 0 dz 0 


(3.12) 


where the boundary conditions on Gi are chosen such that the surface inte- 
gral vanishes. The surface integral in (3.11) can be written for the region 

R as follows: 
n 


/ { Zj [< P n- p no) C-^} ( 

-b u 3xo 


+ / f j t( p n- p no) ( + — ) ( 

-a 0 3yo 


* / n< p n- p no>U*— ) - G l ( 

_b u n v- 


3Gj 

3Xo 


a b 


3G, 


+ ! f J (( p n - p no> (-— ) - G > C 

-a 3yo 


* L /„ [<V p no> (-rf) ' G i ( 

-a -b 3 z 0 


-3 (P -P 
' n 

„o> 

3Xo 


3 (P - p 

) 

n 

no 

ayo 


-3(P -P 

nn ) 

n 

no 

3Xo 


-3( p -P 

) 

n 

no 

3yo 


-3(P -P 

n ) 

n 

n o ■ 


■)] dy 0 dzo 
xo=-a 


■)] dx 0 dz 0 
y 0 =b 


•)] dy 0 dz 0 
xo = a 


■)] dx 0 dz 0 
yo=-b 


(3.13) 


_)] dxodyo 
8Z 0 zo=0 
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+ / /„[ < VW (^) - G > ( 


3(P _P ) 
v n n f 


-a -b 


3Zr 


3 z 0 


-)] dx 0 dy Q = 0 


Z »“ z j 


The above equation gives the boundary condition for the Green’s function and 
we obtain that Gj must satisfy 


+ ii§L + liiL - Gj /L 2 = 6(x 0 -x) 6(y 0 -y) 6(z 0 -z), x 0 ,y 0> z 0 eR n (3.14) 

3x 0 2 3yg 2 3 Z 0 2 


with boundary conditions: 


0 p 3G]_ = SqGi at Xq = _ a> 3G^_ = Sq6i at 


x n =a 


3x r 


3X r 


D„ = SoGj at y 0 * -b, -U 0 = SqGj at y 0 =b 1 


(3.15) 


3yo 


3yo 


0 = S Gj at z 0 = 0, G x = 0 at z 0 = z- 

h 3z ^ J 


We perform a similar analysis for the region and find the Green's func- 


tion G 2 must satisfy 


111. + _ g 2 / l 2 = 6(x 0 -x) fi(y 0 -y) s(z 0 -z), x 0 ,y 0 ,z 0 eR (3.16) 

3X 0 2 3yg 2 3Z 0 2 

with boundary conditions 
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- 0 . 


3.G 2 

- s x g 2 

at 

x o 

= a. 

D n 

3G 2 

= Sj Gq at xq - - a 

X 

o 





3*0 


3 G 2 

= SjG 2 

at 

y 0 

= b. 

0 

n 

3 G 2 

= SjG 2 at y 0 = -b 

3y 0 





3y 0 


3 G 2 

" S n G2 

at 

2 o 

= H, 

g 2 

= 0 

at z 0 = z.+o) 

3 Z 0 






J 

> 


(3.17) 


and the solution for the minority carrier density in the R^ region is given 

by 


H b a 


-a<p 0 (l-ReJe *^ 0 


n p n Po ^Z-+u -b -a D 

j n 


G 2 (x,y,z; x 0 ,y 0 ,z 0 )dx 0 dy 0 dz 0 (3.18) 


From the relations (3.12) and (3.18) we can obtain the current 
densities 


J p ■ -«°p ~ ( P n- P P0> at l -h J n * <#>„ ^ VV at Z ’ Z j + " (3 - 19 > 


To find the Green's functions Gi , G 2 for the regions R n and Rp we examine 
the eigenfunctions and eigenvalues associated with the operator L(). We let 
u=u(r) satisfy the eigenvalue problem 


d 2 u 

dr 2 


+ A 2 u = 0 


(3.20) 


with boundary conditions 
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du du 

D — + SqU = 0 at r=h, D - S,u = 0 at r= -h (3.21) 

dr dr 


The eigenfunctions of this problem are: 


u = u n (r) = cos (A n r) = u (r; D, S 3 , h, x n ) (3.22) 

where A = A satisfies the equation 
n 

(Ah) tan (Ah) = Sah (3.23) 

n n D 

We let ? n = A n h satisfy tan(c n ) = S 3 h/K n then is characterized by the 
intersection of the curves y = tan(c) and y = S 3 h/D^ (See Figure 5). Note 
that for large values of n, c approaches rnr, and that the above .functions 
u^Cr) are orthogonal on the interval (-h, h) 


and 


h 0 , m*n 

(u , u ) = / u (r)u (r)dr = { , 

n’ m u n m 1 « u n ii 2 , m=n 


m 


(3.24) 


where 


"V 2 


h + cos 2 (A h) = h + 
DA* 


S 3 h2 

D 


cos 2 (c n ) 



(3.25) 


We will use these functions in the construction of our Green’s function 
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FIGURE 5. Graphical scheme for determining eigenvalues 
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solutions for G x and G 2 . 

Define the eigenfunctions associated with the region R^ as 

U (x) = cos(A x), -a <x<a, m=l,2,3... 

m m 

with 


■ , l2 , s oa2 cos2(5 J 
= V’ ,U n ,’ 2 * a + -T- --- 


«ra 2 


and where 


5 m tan c m = S 0 a/D p , m-1,2,3, 


Also define the eigenfunctions 


V m (y) = cos(u m y), -b<y<b 


with 


”m ' %»’ ,v m ' 2 = b + 


S 0 b2 cos2 ("m) 

°p V 


where 


n m tan = S o b/O p’ 

There will be a similar set of eigenfunctions associated with the region R . 

P 

★ ★ 

We denote these eigenfunctions by U (x) and V (y) and note that S 0 is 
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replaced by S x in the calculation of these eigenfunctions. We now assume a 
solution for the Green's functions Gi , G 2 as: 


co 00 


Gi = E E A.- -j ( zo ) U-j(xo )V.- (yo ) 
1=1 j=l J J 


(3.26) 


00 00 


G 2 = E E B . . ( Zq ) U.(x 0 )V.(y 0 ) 
i=l j=l J J 


(3.27) 


as the functions and V.. satisfy the boundary conditions and 


U n (x) = cos (x p h) = u(x; D , S 0 , a, * n ), -a<x<a 


U n (x) = cos (A x) = u(x; D n , S ls a, A n ), -a<x<a 


V m (y) = cos (u m y) = u(y; 0 D , So, b, %), -b<y<b 


P’ 


V m (y) = cos (u m y) = u(y; 0 n , S lt b, ym ), -b<y<b 

where it is assumed that So, Si are different from zero. 'Substituting 

(3.26) and (3.27) into (3.14) and (3.16) respectively, we find that A mn and 

B must be chosen to satisfy the conditions: 
mn 


» 1 U i (x)V i (y)6(zo-z) 

A.. - -JL A.. = J , 0< z 0 <Z 

1 J Y . . 2 1 J n 1 1 11 2 d \i n 2 J 


(3.28) 


1J 


nU,ir n V- im 
3 
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dA. . 

*<<(!,) =0, 0„ - S„ A„ ■ 0 »t z, • 0 


TJ J' 


Pdz, P ’J 


with 


Y. . 2 = L 2/(1 + L 2 (x? + y 2 )) 

1 J P PI J 


11 1 * 1 
B,, - _ 1 . B • • = J 


U- (x)V i (y)S(z 0 -z) 


’J a. -2 1J 
1 J 


n U* n 2 n V*n 2 
J i 


, Zj+03 < Zg<H 


dB. . 

B- .(Z.+w) = 0 , D n _li+ S n B.. = 0 at z 0 = H. 


ij J 




with 


iJ 


i 2 • L n 2 /d + L n 2 (V + »I*H 


These equations have the solutions 


U J -(x)V 1 -(y)y 1 (z 0 )y 2 (z) 


A ij(zo) =1 


ii U ^ n 2 ii V .(i 2 Aj 

Uj(x)V.(y)y 1 (z)y 2 (z 0 ) 

II U. II 2 II V. II 2 A, 

J i 1 


, 0<z 0 z 


, z< z o < Z . 


where 


yi(z) = cosh ( z/T i j ) + ( SpY i j/D p ) sinh (z/^j) 


sinh (Z .-z)/y- •) 

y 2 (z) > £ ££_ 

s1nh<z j /1r fj) 


( 3 . 29 ) 


( 3 . 30 ) 
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and 


Also, 


Y u sinh ‘W’ 


B u ■ 


C U*(x)V*(y)y t (2)y 3 (2 0 ) 

mj*n 2 iiV*ii 2 a 2 

J ' 

U*(x)V*(y)y 3 (z)y 4 (z 0 ) 

ll U*ll 2 fl V - II 2 A, 

V. J T 


, Zj +to<z Q <z 


, z<z 0 <H 


(3.31) 


with 


y 3 (z) = sinh (Z-Z-j-* 0 )/^) 

■ S ct . . 

y 4 (z) = cosh ((H-z) / CT . .) + n . sinh ((H-z) /cr . .) 

’ J n ' J 

n 

a 2 = - — 

CT. . J 


From the equations (3.12) (3.18) (3.19) (3.26) and (3.27) we obtain the 
current density 


OO 00 


J n = 4<xq<t>o ( 1 - R) Z £ y 


m-l n=l 


mn 


U (x)V (y)sin(x )sin(u )G 
n m J na 'mb' mn 

,u n , 2 , V 2 


(3.32) 
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with 


^mn 


a 2 y 2 -1 
' mn 


S P T 1.1 /0 P <,,1f mn- e 7p Y i/y° sM ¥ y mn> tsinh 7 /y mn7 Y mn e J 


*< Z j> 


Observe that G has a removable singularity where ay =1 and in this case 
mn mn 

we can write 


mn 


ay 


mn+1 


i+iyij/yn 


1 - e 


-Z ,/y (ay -1) 
j mn v mn 


ay, -1 

mn 


-e 


- z j 


y >( z j) 


The current density J n is given by 


u n< x ' V m< y > V n <V> sin (“m b > H »n 


OO 00 


J = 4aq<J) 0 (l-R) j; j; 


m 1 n=1 llvVnvV 

n in 


* * 
J m 


A n ^(Zi+w) 


with 


(3.33) 


mn a 2 a 2 -1 
mn 


u S a . . 
e -« H ( n ij 


a a 


mn 


-a ( z ,+w) S a . . 

) + e J ( - H . -V , a -1) sinh(H7a ) 
J - mn 1 v mn 


D 


mn ' 


-a(Z -+w) 


+ e 


f°mn a ' S n °ij /0 n ) cosh ( H ' /o <J 


with H' = H-zj-w. The total light current is = Jp+O n . 
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DARK CURRENTS 


The equations describing the dark currents associated with the n and p 
regions may be expressed 


72 < Vno> - * 0 

P 


with boundary conditions: 


P = P -P = P 
n no no 


qV i 

exp ( — =-) -1 
KT 


at z=z- 

J 


3P 

D — - S P = 0 at z = 0 
p 3z p 


3 P 3 P 

-D S 0 P = 0 at x - a, D — - S 0 P = 0 at x= -a 

p 3x p 3x 


3 P 3 P 

-D __ - S 0 P = 0 at y = b, 0 n _ - S 0 P = 0 at y= -b 


ay 


ay 


and 


v 2 (n -n„ 0 ) 
' p po ' 


<w> 


= 0 


(3.35) 


with boundary conditions 


n = n -n = n 
p po po 


qV i 

exp ( — i) -1 
KT 


at Z=Z.+w 

J 
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-D - S n = 0 at z = H 
n 3z n 


-D !!!.- S,n = 0 at x = a. 0 l!l - S,n = 0 at x=-a 
n 1 ’ n — 1 

' 3x n 3x 


3n 


3n 


-D — - Sjn = 0 at y = b, 0 - Sjn = 0 at y=-b 

3y 3y 

with V. the junction voltage. 

J 

We use the relation (3.11) with g=0 and obtain the Green's second 
formula 


/// (P _|3 no) dT o - P n "^no ”. J ' 1 ^i 

R n i-l 


(3.36) 


where I ^ , i =1 ... .6 are the surface integrals defined by (3.13). Applying 
the boundary conditions from (3.34) the relation (3.36) reduces to 


/ P no 

* d 

qV J 

exp ( -1 

3Gj 

■KT 

3Z 0 


- 



Zq =Z -: 


dx 0 dy 0 

j 


(3.37) 


and the dark current opposing the photo current is given by 


J , = -q D grad (P -P„ ) • k 
pd ^ p 3 ' n no' 


(3.38) 


z “ z j 


A similar analysis for the region produces the electron density 
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b a 

w [ b [ a -v 


‘ IV, 

exp (_i.) -1 
KT 


3G 2 

3Z 0 


dx 0 dy 0 (3.39) 


and the dark current opposing the photocurrent is given by 


J nd = -qD n gnd (n p -n pQ ) • (-k) 


z=Z,+w 

J 


(3.40) 


The equations (3.38) and (3.40) simplify to 


J pd ” 4abq D o P nn( e J 'I)” 


q Y kT 11 . , U n (x) Vy>Vi (z j> sin <V )s1 "<% b > 


p no 


m n NU n ii 2 nV m ii2 yi ( Z .) (x p a) ( % b) 


(3.41) 


and 


* * i 


J nd s 4abq D n n pof e J - 1 ) r E 


qV j /kT n r r U n (x)V m (y)y 1+ (Z j+ w)sin(x n a)sin( V|n b) 


(3.42) 


m ° "uVnvV y 4 (Z.+w) U*a) ( p *b) 


where P = n 2 /N , and n = n?/N 
no i d po i a 

The total dark current density is 


J dark = J pd * J nd ‘ °» ( exq £r> - 1 ) (3.43) 

K I 

We use the recombination current density from (1.21) and calculate the 
total currents by integration 


26 



(3.44) 


b a 


sc 


= / / J. dxdy 

-b -a L 


I = J • 4ab 
rec rec 


b a 

Io = / / Jo dxdy 

-b -a 


b a 

*dark ~ [ a d dark dxd ^ 


(3.45) 


(3.46) 


(3.47) 


The current voltage relation is obtained from the equivalent circuit 
diagram of Figure 2. The computer program "SC3D" which utilizes the above 
equations is given in Appendix C. 

4. THREE-DIMENSIONAL VERTICAL JUNCTION SOLAR CELL 
With reference to the geometry illustrated in Figure 4, we assume the 
photoconverter is illuminated from the side x=a. For this assumption, we 
use the generation term g(a-x) where g(x) is defined by equation (1.1). 

Also we must modify the Green's functions G 2 , G 2 used in the conventional 
photoconverter. Note that Gi , G 2 were expanded as double Fourier series 
which were symmetric in both x and y. For illumination from the side x=a we 
must remove the syimietry in the x-direction. Consequently, for the n-region 
and p-region we construct Green's functions g 1? g 2 having the following 
form: 


J ” Vm^o) W £^ z o)yy) w ^ z ) li; 3m& 

m=l £=1 'V 2 , V 2 p mi 


(4.1) 
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★ * ★ ★ 


00 oo 


g 2 = e r 

m=l £=1 


y*o) M t (z o> v »<*>*i t < 2 >*4 01 t 


,v m , 2 *V 2 %l 


(4.2) 


where 


W £ = sin (n £ (Z j -z)), 0 < z < Zj, £=1,2,3. 


(4.3) 


and 


W £ = sin (n £ (z-Zj-w) ) , Zj+io<z<H, £=1,2,3. 


(4.4) 


are eigenfunctions with eigenvalues determined by the roots of the equations 


-0 (n Z.)/(S Z.) = tan (nZ.) 
P * J P J l J 


and 


(4.5) 


" D n (n £ H ' )/S n H,) = tan * H' =H-Zj-ni (4.6) 


respectively. These eigenfunctions have the norm squared 


i 2 = I Z. - I 
2 J 4 


sin(2n.Z .) 
£ J 


(4.7) 


nW*n 2 


i H , . i sin < 2 "t K '> 


(4.8) 


In (4.1) and (4.2) V m , V m have been previously defined in (3.26) and (3.27) 
Also in (4.1) and (4.2) we have 
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1 So 


2 r . . 
u 


p. . = (J_ + I Li) sinh (i!_) + !5o. cosh (£i_) 


• 2a 


ij 


F ij D ' 


r. . 


r ij 


★ 

ij 


1 , Sl ‘ 


e . . 


e . . D 2 

ij 


pi = (_L_ + _ ii) sinh ( 2a ) + ^LL cosh (£i_) 

D 


■2a 


e ij 


n 


3mn 


e . .2 = I 2 

» ij L n 

/(l - L n 

-yi(X)y 2 (Xo), 

-a<Xo<X 

-yi(Xo)y 2 (X), 

X<Xo < a 

->4(X)y3(X) , 

-a<Xo<X 

-^4 (Xo )ys (X) , 

X<Xo < a 


e ij 


★ ★ 


*4U = 


y^X) = cosh (i-Jl) + ii sinh (LLl^.) 

r D r 

mn p mn 


Y 2 (X) = yi (-X) 


y 3 (X) = cosh (±i) + lilii sinh (— ) 


£ . . D 
i j n 


e . . 
1J 


*(X) = ys(-x) 

In the n-region we can represent the solution to the boundary value 


(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 

(4.16) 

(4.17) 
problem 
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(3.1) by 


1. w 

J b d n(dX) 

P n" P no = / / / — 2 — 9i (x,y,z; x 0 ,y 0 ,z 0 ) dx 0 dy 0 dz 0 


(4.18) 


and in the p-region we can represent the solution to the boundary value 
problem (3.5) by 


H b a , . 

n D" n DO = iff 92(x.y»2; x 0 ,y 0 ,Zo) dx 0 dy 0 dz 0 (4.19) 

iu -b -a D_ 

J+U) 


where the source term has been replaced by g(a-Xo) which models the illumi- 
nation from the side X Q = a. From these equations we use (3.19) to obtain 
the photocurrents. The computer program "VSC3D" describing the vertical 
junction model is given in the Appendix D together with representative 
graphic output. 
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APPENDIX A 

BASIC COMPUTER PROGRAM FOR 
ONE-DIMENSIONAL CONVENTIONAL PHOTOCONVERTER 


A-l 



10 
20 
30 
40 
50 
SO 
70 
SO 
90 . 
100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 . 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
6 1 0 
620 


OPTION BASE 1 

f **************************************************************** ***** 

I 

I PROGRAM CONVEN 


CONVENTIONAL SOLAR CELL PROGRAM SET UP FOR Si SOLAR CELL CHARACTERISTICS 


REM CONVENTIONAL SOLAR CELL PROGRAM CONSTANT TEMPERATURE 

t 

I «*•***•»•*■**■»■»*■*••*•**■»•*■**■•*■**•#••*■*■»■***•<*•*■*■*•*•*•*■*••* ■*•*■*■*■*•■*■* •*■***■»*■«■■*•*■***■*■*•*»•*■#* •*■*■*•■*■*•■* 

DIM X<50> ,Xx<50> , Aa<50 , 4) , Z<50) , Vol ts< 1 00) , Amps (100) 

DIM Esi(20> ,Asi<20> 

t ********************************************************************* 


♦ 

I READ IN Si Icon ABSORPTION COEFFICIENT DATA 

Nsi=20 ! THIS IS NUMBER OF DATA POINT PAIRS TO READ 

♦IF YOU CHANGE THIS NUMBER CHANGE DIM STATEMENT ABOVE 

FOR 1=1 TO Nsi 

READ Esi<I) , Asi < I> ’READ IN ENERGY <EV) . AND ABSORPTION COEFF CM-1 

NEXT I 

DATA 1.1,1,1.15,10,1.24,100,1 .378,350 . 1 .46,950 , 1 . 77 , 2000 
DATA 2.067,4500 ,2.48, 1E4, 2.76, 2E4, '3. 1 , 5E4 .3 . 26 , 1 E5 . 3 . 5 , 1 E6 . 3 . 6 , 1 . 1 Eb 
DATA 4.1 ,8E6, 4.2,2. 5E6 , 5,1 .9E6, 5.2, 2E6. 6, 1 .5E6. 8,1 ,3E6, 10,1 .1E6 
GINIT 

GRAPHICS ON 
GCLEAR 

WINDOW 0,100,0,100 
CSIZE 9 
LORG 5 
MOVE 50,50 

LABEL "SPECTRAL RESPONSE CURVE" 

CSIZE 5 
MOVE 50,40 

LABEL "IT TAKES A MOMENT TO CALCULATE” 

! SOLAR CELL CONSTANTS 
f CONSTANTS FOR MOBILITY Man 

A0=65 . 02 ! FOR Si; FOR GaAs A0=5076.11 ; FOR Ge A0=458.68 

A 1 =5 . 72E-9 IFOR Si;FOR GaAs A1=6.03E-8 ; FOR Ge A1-7.92E-9 

A2=2 . 49E+9 IFOR Si;FGR GaAs A2=2.13E6 ; FOR Ge A2=6.08E7 

A3=2 . 5 IFOR Si ; FOR GaAs A3-1.Q ; FOR Ge A3=1 .66 

ICONSTANTS FOR MOBILITY Mup 

B0=54 . 46 IFOR Si ; FOR GaAs B0=37.33 ; FOR Ge B0=3140.27 
B1 =6 . 76E-9 IFOR Si; FOR GaAs B 1 =8 . 51 E-3 ; FOR Ge B1-1.79E-7 
B2=2 . 93E+9 1 FOR Si; FOR GaAs B2=6.05E7 ; FOR Ge B2-1.18E9 
B3=2 . 7 IFOR Si; FOR GaAs B3-2.1; FOR Ge B3=2.33 
ICONSTANTS FOR BANDGAP AS FUNCTION OF TEMPERATURE 
Eggo-1.16 IFOR Si; FOR GaAs 1.522 : FOR Ge .741 

Alphao=7 . 02E-4 IFIR Si; FOR GaAs 5.8E-4; FOR Ge 4.56E-4 

Betao=1 1 08 . IFOR Si; FOR GaAs 300: FOR Ge 210 
ICONSTANTS FOR INTRINSIC CARRIER DENSITY 
N ial =4 . 9E+1 5 
Nia2=1 .61 


B=7. 1E+15 
Xk k o= 1 .2E-3 
Xkk 1 = .95 
T0=300 

Kb=8 . 6 1 72E-5 
0=1 .6E-19 


ICONSTANT FOR LIFETIME CALCULATIONS 
ICONSTANT FOR ABSORPTION SHIFT 
ICONSTANT FOR ABSORPTION SHIFT -SEE BELOW 
I FIXED REFERENCE TEMPERATURE 
I BOLTZMANN'S CONSTANT EV/K 
I ELECTRON CHARGE 


t r***-**-**-*"*-*-*-***--*-****-*-*-**-*-*^ ■*■■*■»*■**■*■■*■*■*■*■**■*■*•*•*•*•*■*••*■**•*■**•*•*,*•*•*■*■*■■» *■*•*■*»■* 


j 

I 


ADJUSTABLE SOLAR CELL PARAMETERS 


♦ 


A- 2 


L • --J 




40 Parameters: 

50 Sn=0 . E+0 
60 Sp=0 . E+0 
70 Na= 1 .25E+1 7 
80 Nd=5 . 0E+1 9 
90 Tn0= I .86E-4 
00 T P 0=3 . 52E-3 

10 Mo- 1.000 


Hh2= 1 00 

Kback = 1.5 
l p ipe= 1 00 

Lambda= 1 . 06 
Re= . 05 
Rs = 0 . 

Rsh= 1 . 0E+6 
Epso= 11.6 
H=4 . 500E- 1 
X j=5 . 0E-4 

f *•***■»•*-*•*■■*•*•*■* 


CHANGE THE VALUES BELOW IF DESIRED 

! SURFACE RECOMBINATION VELOCITY CM/SEC 
{SURFACE RECOMBINATION VELOCITY CM/SEC 
! DOPING DENSITY P REGION CM-3 
! DOPING DENSITY N-REGION CM-3 
! CONSTANT FOR LIFETIME CALCULATION 
! CONSTANT FOR LIFETIME CALCULATION 
1 LASER POWER KW/CM2 
! Wo DETERMINES TEMPERATURE IN CELL 
’SEE BELOW 

1HEAT TRANSFER COEFFICIENT TO BACK SURFACE 
’UNITS OF WATTS/CM2-DEG C 

’WATTS/CM-DEG C THERMAL CONDUCTIVITY BACK SURFACE 
’TEMPERATURE DEG C OF HEAT PIPE (ASSUMED TEMP 
’OF HEAT PIPE LOWER TEMPERATURE 
! LASER WAVELENGTH IN MICRONS 

’REFLECTION COEFFICIENT FOR ILLUMINATED SURFACE 

* SOLAR CELL SERIES RESISTANCE IN OHMS 
’SHUNT RESISTANCE IN OHMS 

♦ Si PERMITIVITY COEFFICIENT (DIMENSIONLESS) 

! SOLAR CELL HEIGHT IN CM 

’JUNCTION DEPTH IN CM 

-*-*-*• •*«■-***■■*■*■ 


» BEGIN CALCULATIONS FOR THE SPECTRAL RESPONSE 

t 

* * •**•*•» ■*•*■**.***■•*•*■•*•*•*•*■*•* •*■*•**•*■* •*■*•*-*•*•*•»•*•*■*■*■*•*■*•*•* •»•*■*•*■»■*-****•***•* ■**■** ■*•*■•* ■*•*•*•*■*■*•**■*■*■ 
’ESTIMATE ABSORPTION COEFFICIENT ALPHA 

T=273+T P ipe+1 000*Wo*( 1 -Re)/Hh2 ’ESTIMATE SOLAR CELL TEMP DEG K 
Ener = 1 .2402/Lambda ! ENERGY IN EV FOR WAVELENGTH LAMBDA 

Egt=Eggo-Al P hao*T*T/ (T+Betao) ’BANDGAP FOR THIS TEMPERATURE 
IF (Ener<=1 .7) THEN Xkko=( 1 .3*Ener-1 >*1 . OE-3 
Deie=2*Xkko*( T-TO) 

E=Ener+Dele 

Alpha=FNAbsorp(Nsi .Esi (*) ,Asi<*) ,E) ’APPROXIMATION OF A1PHA 

! CALCULATE TEMPERATURE OF SOLAR CELL AT LASER SURFACE INTERFACE 
Yy 1 =Alpha*H 
Yy2 = 0 

IF ( Yy 1 <220 ) THEN Yy2=EXP(-Yy ! > 

T = Tpipe+( Wo*1 000* ( 1 -Re) /Kback >*( (Kback /Hh2) - 1 /Alpha) *( 1 - Yy2) 

IF (T<Tpipe) THEN T=Tpipe 
T=T+273 ! CONVERT TO DEG K 

Kbtg=Kb*T ’ BOLTZMANN'S CONSTANT TIMES TEMPERATURE 


-Yy2> 


Eg t=Eggo-Alphao*T*T/( T+Betao) 
Tn=TnO/( 1+Na/B) 

Tp=TpO/( 1 . +Nd/B) 

Mun=FNMan(Na, T , AO ,A1 ,A2,A3) 
Mup=FNMup(Nd , T , BO , B1 ,B2.B3) 
Dn=Kbtg*Mun 
Dp=Kbtg*Mup 

Phio=Wo*Lambda*5 . 03306E+21 I 


! BANDGAP CALCULATION 
LIFETIME CALCULATION 
LIFETIME CALCULATION 
’ MOBILITY CALCULATION 
MOBILITY CALCULATION - 
DIFFUSION COEFFICIENT Ci 
DIFFUSION COEFFICIENT C 


-SEE BELOW FOR 

Mup=FNMup(Nd , T , BO , B1 .B2.B3) ’ MOBILITY CALCULATION -SEE BELOW FOR FN 

Dn=Kbtg*Mun * DIFFUSION COEFFICIENT CALCULATION 

Dp=Kbtg*Mu P ! DIFFUSION COEFFICIENT CALCULATION 

Phio=Wo*Lambda*5 . 03306E+21 ! PH0T0NS/CM2-S ' 

Ni=(Nia1 >*(Nia2)*(T'1 .5)*EXP(-Egt/(Kbtg*2) ) ’INTRINSIC CARRIER DENSITY 
Vbi=Kbtg-*LOG(Na*Nd/(Ni*Ni ) ) ’ BUILT IN VOLTAGE 
W2=(2.*Epso*(8.85E-14)/Q)*Vbi*(1 ./Nd+1 ./Na) 

W=SQR(W2) ’SPACE CHARGE REGION WIDTH 

Hp=H-Xj-W ! DEPTH OF BOTTOM SECTION (P REGION) 

Lp=SQR(Dp-*Tp) » DIFFUSION LENGTH CALCULATION 
Ln=SQR(Dn*Tn) ! DIFFUSION LENGTH CALCULATION 

Dl= . 075 ! WAVELENGTH INCREMENT FOR SPECTRAL RESPONSE CURVE 


’SPACE CHARGE REGION WIDTH 
’ DEPTH OF BOTTOM SECTION (P REGION) 

’ DIFFUSION LENGTH CALCULATION 
! DIFFUSION LENGTH CALCULATION 
! WAVELENGTH INCREMENT FOR SPECTRAL RESPONSE CURVE 


A-3 


1220 

1230 

1240 

1250 

1260 

1270 

1230 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 

1760 

1770 

1780 

1790 

1800 

1310 


| ********************************************************************** 

I 

I PLOT SPECTRAL RESPONSE AND THEN PAUSE TO LOOK AT IT 
} 

! ************************************************************************ 
FOR J-1 TO 50 

X < J ) = 1 . 1 1 5+< J- 1 >*D1 ! THIS IS ENERGY IN EV 

Xx(J)=1 .2402/X<J) ’THIS IS WAVELENGTH IN MICRONS 
Ener=X( J) 

Engo=X( J) 

IF (Ener> 1.7) THEN Xkk=Xkko 

IF <Ener<1.7> THEN Xkk=<Xkk 1 *Ener-l )*] . OE-3 

Dele=2*Xkk*< T-TG ) 

Ene=Ener+Dele 

Alpha=FNAbsorp<Nsi ,Esi<*),Asi(*> , Ene) 

1PRINT T,Ener .Alpha, Dele 
Aal =Sp*Lp/Dp 
Aa2=Xj/Lp 
Y1 =Alpha*Lp 
AaO=Aa1 +Y1 
Y2=Y1/(Y1*Y1-1 ) 

F j=FNJ<Aa1 ,1 . , Aa2) 

Z1 =Alpha*X j 
Z2=A lpha*< X j+W) 

Z3=Alpha*W 

Z4=Alpha*(Hp) 

IF (ZK200) THEN GOTO L53 
Y 1 z = 0 
GOTO L63 

L53 : Ylz-EXP(-ZI) 

L63 : IF < Z2<200 ) THEN GOTO L54 

Y2z = 0 
GOTO L64 
L54: Y2z=EXP(-Z2) 

L64 : IF <Z3<200> THEN GOTO L55 
Y3z = 0 
GOTO L65 
L55: Y3z=EXP<-Z3> 

L65 : IF (Z4<200) THEN GOTO L56 

Y4z = 0 
GOTO L66 

L56 : Y4z=EXP(-Z4> 

L66 : Y3=Aa0-Fj*Y1 z 

Y4=FNH(Aa1 ,1 . ,Aa2> 

Y5= Y 1 *Y 1 z 

Aa< J ,1 )=Y2*< < Y3/Y4)-Y5> ! TOP OF JUNCTION 

Aa3=Sn*Ln/Dn 
Aa4=Alpha*Ln 
Aa5=Hp/Ln 
Aa6=Aa4~Aa3 
Y6=Aa4/<Aa4*Aa4-1 > 

Y7=Aa6*Y4z 
F j=FNJ< Aa3 , 1 . ,Aa5> 

Y7 = Y7+F j 

Y8=FNH<Aa3 , 1 . ,Aa5> 

Aa( J , 2) =Y6*Y2z*( Aa4- Y7/Y8) 

Aa( J , 3 ) = Y 1 z*<1 . -Y3 z ) 

Aa< J , 4) =Aa< J , 1 ) +Aa< J . 2> +Aa< J , 3) 

NEXT J 

j ***************************************************************** 
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1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1 980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
21 1 0 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 

2390 

2391 

2392 


! GRAPHICS FOR SPECTRAL RESPONSE CURVE 

I 

t *»**•*•»*•■*■•****•*••«■■*••**•*•*•»■*•»•»«■*•*■*•»■•»•**■»■***•*••*•«■*■*•*■*■*■* ■»•**•*■*•*■**•*■■***■*■**■■**•*•*•*■** 

GINIT 

GRAPHICS ON 

DUMP DEVICE IS 701 !FOR LARGE GRAPHS USE 701 .EXPANDED 
WINDOW -.25,1 .60, -.1 , 1 .1 ! XMIN , XMAX . YMIN , YMAX FOR GRAPH 

FRAME ! DRAW A BOX AROUND GRAPH 

AXES .10, .10,0. ,0,5, 2, 6 ! XTIC , YTIC , ORIGIN , POSITION OF MAJOR TICS 

f *•*•***■*■*■*•*■»■*■»•»*■»•*■**•*■*•)*•**•*■»•»■»■*•*«•*■*■)»•*■*■»*■*• -»***r* If « -“**■*•■*■*•**■****-*■* 

t 

! NOW TO LABEL THE AXES— -X AXIS FIRST 

t 

f ********•*•*■*»*•»*•»•****-»•*****■***•»•*■»*■**■*■*■*■*•*■*•»■»■*••*•*■■*■*•*•■**•!*■*•*■*■*■*■■*■*■*■» *•*■* *•«••* 

CSIZE 4 
LORG 6 

FOR X_label=0 . TO 1.5 STEP .5 

MOVE X_label ,-.045 

LABEL X_label 

NEXT X_label 

LORG 8 

FOR Y_label=0 TO 1 STEP .1 
MOVE - . 03 , Y_label 
LABEL Y_iabel 
NEXT Y_label 
LORG 6 

f *K***»«***»»*»*it**«*j(«ji«*****************«**»*i(*««*********«#**« 

\ 

♦LABEL THE AXES WITH NAMES 


DEG ! THIS SETS ANGLES TO DEGREES 

LDIR 0 ! ZERO ROTATION FOR DIRECTION OF LABELS 

MOVE ( .75) ,-.01 

LABEL "WAVELENGTH < MICRONS > " ! WRITE LABEL IN CENTER OF AXES 
LDIR 90 ! GET READY FOR A 90 DEG ROTATION DF LABEL 

MOVE -.18,. 5 ! MOVE TO CENTER OF Y AXIS 

LABEL "SPECTRAL RESPONSE" 

| *-*-»*****»**#*****»**#*K*#*«»***#**-»-****-»»»*iHHH**##**«*»-»*»*i<-»»*#* 

j 

! PLOT THE SPECTRAL RESPONSE CURVES FOR 

f 

! CONTRIBUTION: FROM FRONT ; BACK : DEPLETION REGION: TOTAL 

f 

t **»*********»**********!(**#*•»*#*****■ S* ■*■*•*■*■»■*■*•■*■*■**-*■•*■*■■*■*■*•»■*•*■*•* r*** *•*■*■*■ 

FOR 1=1 TO 4 

FOR J-1 TO 50 
Z< J)=Aa< J , I) 

NEXT J 

MOVE X x O ) , Z ( 1 > 

FOR K= 1 TO 50 
DRAW Xx (K ) , Z(K) 

NEXT K 
NEXT I 
LDIR 0 
CSIZE 3.8 
LORG 3 

MOVE -.15, -.05 
LABEL "Si" 


010 
020 
030 
040 
050 
060 
070 
080 
090 
100 
1 1 0 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
290 
300 
310 
320 
330 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 


J01 =<Q*Dp/Lp)*<Ni*Ni/Nd)*<Z1 /Z2> 

J02= (Q*Dn/Ln>*< Ni*Ni/Na>*< Z3/Z4) 

J0=J01 +J02 

Z3=Alpha*Lp 

Z4=Alpha*Ln 

R1 =Q*PhiO*< 1 -Re)*Z3/<Z3*Z3- 1 ) 

R2=Aa1 + Z3 
R3=EXP(-Aipha*Xj) 

Jp=R1 *< ( <R2-R3*Z1 >/Z2>-Z3*R3> 

R4=EXP<-Aipha*<X j+W) ) 

R5=Q*Phio ~ < 1 -Re)*Z4*R4/<Z4*Z4-i > 

R6=FNCosh(X2) 

Yy 1 “Alpha*Hp 
Yy2=0 

IF ( Yy 1 <220 ) THEN Yy2=EXP< -Yy 1 ) 

R7=Z4*Yy2 

R8=Z4-<Aa3*<R6-Yy2) +FNSlnh(X2)+R7)/FNH( Aa3 , 1 . ,X2> 

Jn=R5*R8 

Jdr=Q*Phio*<1-Re)*EXP(-Aipha*Xj>*< I . -EXP< -Alpha*W> ) 

Jsc= Jp+Jn+Jdr 

{ **************** ************ ************»************-***-*•*■•**■*'** ■*•■*••*•■*■*■ 

| 

’ VOLTAGE GOES FROM 0 TO ITS OPEN CIRCUIT VALUE 

! AND WE CALCULATE THE CORRESPONDING CURRENT 

f *******^******4********* ******************** ********4** ******-j*-*-*'*-*-*'*-‘»-*‘* 

I = Jsc 

Begin: X1x=V+I*Rs 

Jrec=FNRec<Ni , W , XI x , Kb , T , Vbi , Tn , Tp> 

IF (Rs=0) THEN GOTO Sk iper 
I=FNN1 (Rs , V , T , JO , Jsc , Jrec , Kb , Rsh > 

GOTO End 

Sk iper : I = Jsc- JO* < EXP < V/(Kb*T> >-1 ) -Jrec-V/Rsh 

End: ! SAVE RESULTS FOR PLOTTING 

Vol ts< Icount >=V*1 000 . ’MILLIVOLTS 

Amps< Icount )= 1*1 000 . 1MILLIAMPS 

P=V*I ! POWER IN WATTS 

IF <P>Pmax) THEN Pmax=P 

IF ( Amps < Icount > >Max ) THEN Max=Amps( Icount) 

IF <Amps< Icount) <0) THEM GOTO Stop 
GOTO Start 

Stop: ’GRAPH CURRENT-VOLTAGE RELATIONSHIP 

Amps < Icount ) = 0 

? 

! GRAPH OF CURRENT VOLTAGE RELATIONSHIP FOLLOWS 

; 

f * *•*•*■»■*•*■*■**•*■<*••«■*•*■*■**■*••*•****■* ■*■*■■*•*■*■*••■*■*■*■** 

GCLEAR 

CSIZE 4 

Maxy=1 .2*Max 

Y 1 =DR0UND< . 1*Maxy,2) 

Dely =Maxy / 1 0 

WINDOW -150, 1 000 , - Y1 *1 . 6 , Maxy + . 2*Y 1 

FRAME 

DEG 

AXES 1 0 , <De ly > ,0,0, 10,2.6 
! LABEL AXES 
LORG 6 
LDIR -90 

FOR X_label= ! 00 TO 1000 STEP 100 
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660 MOVE X_labei.-12*Y1/16 
670 LABEL X_label 

680 NEXT X_label 

690 LORG 8 

700 LOIR 0 

710 FOR Y_label=0 TO Maxy STEP Dely 

720 MOVE 0 , Y_labe 1 

730 LABEL DROUND< Y_label ,2) 

740 NEXT Y_labe i 

750 LORG 5 

760 MOVE 350,-Y1*1 .3 

770 LDIR 0 

780 LABEL "MILLIVOLTS" 

790 MOVE - 100, Maxy /2 

800 LDIR 90 

810 LABEL "MILLIAMPS" 

820 CSIZE 4 

830 LDIR 0 

840 MOVE Vol ts< 1 > , Amps< 1 ) 

850 FOR I j= 1 TO Icount 

860 DRAW Volts<I j> ,Amps<I j) 

870 NEXT Ij 

880 CSIZE 4 

890 Xpo=300 

900 Ypo= . 7 

910 LORG 5 

911 MOVE 400 . ,95*Maxy 

912 LABEL "CONVENTIONAL Si PHOTOCONVERTER" 

920 MOVE Xpo , < Ypo)*Maxy 

930 LABEL "Rs=";Rs 

940 MOVE Xpo , < Ypo- . 04>*Maxy 

950 LABEL "T =" ; DROUND< T-273 , 4) ; "DEG C" 

960 MOVE Xpo, (Ypo-.08)*Maxy 

1970 LABEL "Sn=";Sn;" Ln=" ; DROUND<Ln ,5) 

198 0 MOVE Xpo, (Ypo- . 12)*Maxy 

5990 LABEL "S P =";S P ;" Lp=" ;DROUND<Lp ,5) 

1000 MOVE Xpo, (Ypo-. 16)*Maxy 

(010 LABEL "Na=” ; Na; " Mun=" ;DROUND<Mun ,5) 

1020 MOVE Xpo, (Ypo-.2)*Maxy 

(030 LABEL "Nd=";Nd:" Mup=" ;DROUND<Mup ,5> 

(040 MOVE Xpo , (Ypo- . 24)*Maxy 

(050 LABEL "Xj=";Xj 

(060 MOVE Xpo, <Ypo-.28>*Maxy 

(070 LABEL ”H =";H 

(080 MOVE Xpo , ( Ypo- . 32)*Maxy 

(090 LABEL "Rsh =";Rsh 

(100 MOVE Xpo , < Ypo- . 36)*Maxy 

(110 LABEL "Wo =" ; Wo 

(120 MOVE Xpo , < Ypo- . 4>*Maxy 

(130 LABEL "Pout=" ;DROUND<Pmax ,5) 

(140 MOVE Xpo , < Ypo- . 44>*Maxy 

(150 LABEL "Ef f = " ;DROUND< 1 00*Pmax/(Wo*1 000) .5) 

(160 MOVE Xpo , < Ypo- . 48)*Maxy 
(170 LABEL "Al P ha=" ; Alpha 
(180 MOVE Xpo , < Ypo- . 52)*Maxy 
(190 LABEL "Lambda=" ; Lambda 
(200 MOVE Xpo . (Ypo- .56>*Maxy 
(210 LABEL " Jsc=" ; Jsc 

(220 J ***-»***-S**-»***-»*-» ■*■**■■*-»**■**•*•*•*■*■**■* »*■*•*•* ******■**:■*■*-*■#* ****** ****** 

1230 ! 

(240 ! PAUSE LOOK AT GRAPH IF YOU LIKE IT TYPE — DUMP GRAPH 
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1250 
1260 
1270 
1230 
1290 
(300 
(310 
1320 
1330 
1340 
*350 
(360 
*370 
*330 
*390 
*400 
*410 
*420 
*430 
1440 
*450 
*460 
*470 
*430 
*490 
*500 
*510 
*520 
*530 
*540 
*550 
*560 
*570 
*580 
*590 
*600 
*6 1 0 
(620 
*630 
*640 
*650 
*660 
1670 
*680 
1690 
*700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1730 
i790 
1300 
1310 
1820 
1830 
1840 


t *■*■**-***-*-* + **■*■**#*****-* Kit*-*--******-*-* ******* **-*-**-******-*--*-**-* : ***-*-* ***-**■*■*■* 


TYPE GCLEAR 


TO CLEAR THE SCREEN 


e * *•** ***** * *■***-*■*•* **** ********* * 


! CLEARS THE GRAPHICS SCREEN 


****** 

PAUSE 
GCLEAR 
LORG 5 
CSIZE 9 

WINDOW 0,100,0,100 
HOVE 50,50 

LABEL "THE PROGRAM IS OVER" 

FOR 1=1 TO 3000 
NEXT I 
GCLEAR - 
END 

t *4* *"**-* -*"**"»*"*'*'»*"*’ ***"*'*«'*'*'*'*'* «*»***»*****«#****** *■»**■**■**•* •*•*■»■*■■**■*•*■■*■*■•**■* •*•*■•*■ 

def FNAbsorp ( Ns i , Esi ( *) , As i < *> ,Energ) 

1 = 0 

IF (Energ<Esi < 1 ) ) THEN Alpha=0 
IF <Energ<Esi ( 1 > ) THEN Stop3 
St3: 1=1+1 

IF ( I>Nsi ) THEN Last3 

IF (Esi (I)<=Energ) AND (Esi (1+1 )>=Energ) THEN RatioS 
GOTO St3 

Last3: Alpha=Asi (Nsi ) 

GOTO Stop3 

Rat io3 : Alpha=Asi < I)+<Asi(I+I > -Asi < I) >*(Energ-Esi ( I > ) /(Esi < 1+1 > -Esi < I ) > 

Stop3: ! 

RETURN Alpha 
FNEND 

def FNJ(A1 , A3 , A2) 

Y=A!*FNCosh<A2)+A3*FNSinh<A2> 

RETURN Y 
FNEND 

t ■*■*■*•■*■**■*•*■•* -+ ♦ ■if-*-*-*-***-**-*#'*-*-*.*-*-#"**-*-*-*-*-*--***-*-*--**-*-*-* *•**•»*■**•*■*•**■*•■»•* *■■*•*■*■•* 

def FNCosn<X) 

Y=<EXP<X)+EXP<-X))/2. 

RETURN Y 
FNEND 

j **■*■*«■■*■*•■**■**••**•■**•*•**■-*•*■*-*■**■*•*•*•*•**■» ■»■»•»•»*■»■»■•»«•■»»•» ■*■*■*•*■**•*■*■»■■*•* *■*■**■**■**••*■*■***■* 

def FNSinh(X) 

Y=<EXP(X)-EXP(-X))/2. 

RETURN Y 
FNEND 

def FNH(A,B,C) 

Y=A*FNSinh(C)+B*FNCosh(C) 

RETURN Y 
FNEND 

} *■**■* ■**■*■*■■»*•*■■*■*■*■* ■*■■*■*•■*■* *-*-*•**♦-*-*■**•*■■* *■*■■*■*■*•*■*■***■*•*■*■)**■■*■*•*■**•*•*-**•*■*•■**•**■*•*■**•**■* 

DEF FNMiin < Nd . T , AO , A1 , A2 , A3) ------ 

Bb1-A1*T*T/<N<T<2./3.)> 

D=Nd*Bb1 *< 1 . - . 5*Bb1 > 

X 1 = A0*< T' 1 , 5)/D 
X2=A2*EXP<-A3*L0G(T) > 

R= ! . /X2+ 1 ./XI 
RETURN 1 ./R 
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850 FNEMD 

860 f »******#*#*#*»******»*#***#*«»#***#**»#■»**#********«****#*#**#*##**#* 
870 def FNMup(Na. T . BO . B1 , B2 , B3) 

880 Xt=B2*EXP(-B3*L0G<T) > 

890 Bb 1 =B1 *T*T/<Na" <2 . /3 . ) > 

900 D=Na*Bbl*<1 .-.5*Bb1 ) 

910 X2-B0»<T A 1 .5>/D 

920 R-1 ./X2+1 ./XI 

930 RETURN 1./R 
940 FNEND 

950 i **#»***#*»#*****»*******************#*******«»»*#«*##*#*»#»#***#*#** 
960 def FNRec(Ni ,W,V,Kb,T,Vbi,Tn,Tp) 

970 J1-(i .6E-19>*Ni*W*PI*FNSinh<V/<2*Kb*T>> 

980 J2=S.QR<Tp*'Tn>*<Vbi-V) / <Kb*T> 

990 Jrec= J 1 /J2 

000 RETURN Jrec 
010 FNEND 

02 0 DEF FNN1 (Rs , V , T, Io , Isc , Irec ,Kb ,Rsh> ! NONLINEAR EQ SOLVER 

030 Ib=Isc-Irec-V/Rsh+Io 

040 Kbtg=Kb-»T 

050 U=.9999 

060 Starter: U1=1+Rs/Rsh 

070’ F=Ib*U*U1-Ib+Io*EXP< < V+U*Ib*Rs) /Kbtg> 

080 Fp=Ib*Ut+Io*EXP<(V+Ib*U*Rs>/Kbtg)*<Rs*Ib/Kbtg> 

090 U1=U-F/F P 

100 Error=ABS(U1-U) 

110. IF (Error<1 .0E-6) THEN GOTO Ender 

120 U=U1 

130 GOTO Starter 

140 Ender: RETURN U1*Ib 

150 FNEND 


A-9 



ORIGINAL PAGE 85 
OF POOR QUALITY 




FIGURE Al. Representative output from the program "CONVEN" 

i illustrating spectral response and current voltage 
; relationship for conventional photoconverter. 




APPENDIX B 

BASIC COMPUTER PROGRAM FOR 
ONE-DIMENSIONAL VERTICAL JUNCTION PHOTOCONVERTER 
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10 OPTION BASE 1 

?Q f ****-**--****-*-*'*-K-**-****r-*-***** J **-* : X-****-****-*** : *-** : *** : * : *-*-********-****--**-*- 

30 ! 

40 ! PROGRAM VJSCP 

50 ! 

60 ! VERTICAL JUNCTION SOLAR CELL PROGRAM 

70 ! 

SO REM VERTICAL JUNCTION SOLAR CELL PROGRAM CONSTANT TEMPERATURE 

90 ! 

1 00 ! ***■*■*■**•■»•*•***•*•**■****•*•*■*'»'*•***•»*'**■»■*■*'■*•*■*'»•*■■*■■*■*•*'**•**•*■*■*■**'*■*■*■***•*■*■*■'*•■**•* 

110 DIM Vo 1 ts< 1 00) , Amps< 1 00) , Yisc<201 ) . Yiaark <20 1 ) , Yrec<201 ) 

120 DIM Aa< 100,3) , Z< 1 00) ,X< 100) , Xx< 100) 

130 DIM Esi (20) ,Asi<20) 

I 40 i )•*#»*##**##***##«»****»*****#**#»**»***##***»######»*#»**»****#* 

150 ? 

ISO ! Si ENERGY <EV) , ABSORPTION COEFFICIENT <CM-1) DATA 
170 ! 

1 80 f k******##*******#********^**#*************^#********************* 

190 Nsi=20 

200 FOR 1=1 TO Nsi 

210 READ Esi < I) ,Asi<I> 

220 NEXT I 

230 DATA t .1 ,t t 1 ,t5,10, l .24,100*1 .378,350 , 1 . 46,950 , 1 .77,2000 

240 DATA 2.067,4500,2.48, 1E4. 2.76, 2E4, 3.1 ,5E4, 3.26, 1E5 

250 DATA 3.5,1E6,3.S, 1 .1E6, 4.1 .8E6, 4.2,2. 5E6, 5,1 .9E6, 5.2. 2E6 

260 DATA 6,1 .5E6, 8,1 .3E6, 10,1 .1E6 

270 ! *"*■■*■*■*■*■* -3** *•»*■»•»'»*■»•»•*■»■»■»•»* »■»■»*•»»•*•*»**•*•*■»■»•**■»*•*•*•*•* 

280 GINIT ! INITIALIZE GRAPHICS 

290 GRAPHICS ON 

300 DUMP DEVICE IS 701 

310 WINDOW 0,100,0,100 

320 LORG 5 

330 MOVE 50,50 

340 LABEL "SPECTRAL RESPONSE" 

350 ! **•*•»■*■*■***■* •*•*•*••*■*■•* ■*■**•*•*■*•*■»•»•»•*•»■* **•*•»■*•«•♦•*'»•*'*'!*•»•*■*•*•*•*•■***■*■*■»•*•■»■**•*•*■*■»••*•*•■*•*■•*•*■* 

360 ! 

370 ! SOLAR CELL PARAMETERS 

380 ! 

390 ? ********-***■*-**■******-**** ******** ************-** **■***-**■*■***-*■ ***-***■* -***-* 

400 Parameters: ! 

410 Wo= 1 . 0 ! LASER POWER IN KW/CM2 

420 ! TEMPERATURE CALCULATED BELOW BASED ON THIS POWER 

430 Re= . 05 'REFLECTION COEFFICIENT 

440 Lambda= 1 . 06 ! WAVELENGTH IN MICRONS 

450 Rs=0 . ! SERIES RESISTANCE IN OHMS 

460 Rsh= 1 . OE+6 ! SHUNT RESISTANCE IN OHMS 

470 Sn=0 . ! SURFACE RECOMBINATION VELOCITIES CM/S 

480 Sp = 0 . 

490 Na-1.25E+!7- ! DOPING DENSITIES CM-3 

500 Nd=5 . OE+1 9 ! DOPING DENSITY CM-3 IN N REGION 

510 B1 =7 . 1E+15 ! CONSTANT FOR Si LIFETIMES SEE Taun AND Taup BELOW 

520 Eggo-1.16 !FOR Si; USE 1.522 FOR GaAs ;USE .741 FOR Ge 

530 Ai P hao=7 . 02E-4 IFOR Si;USE 5.8E-4 FOR GaAsiUSE 4.56E-4 FOR Ge 

540 ■ Betao= 1108. IFOR Si;USE 300 FOR GaAs; USE 210 FOR Ge 

550 ! SOLAR CELL DIMENSIONS 

560 A= 1 . 0 ! CM Y=A IS TOP, Y=0 IS BOTTOM 

570 B= . 002 ! CM X=0 IS LEFT, X=B IS RIGHT SIDE OF SOLAR CELL 

580 Xj=.001 ! CM JUNCTION LINE AT X = Xj 

590 IF <Xj>B> THEN PRINT "PLEASE MAKE JUNCTION INSIDE OF CELL-THANK YOU" 
600 Kback = 1.5 ! WATTS/CM-DEG C THERMAL CONDUCTIVITY FOR Si 
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SPECTRAL RESPONSE CURVES 


Kb tg=Kb*T {BOLTZMANN'S CONSTANTS 

Egt=Eggo-Alphao*T*T/< T+Betao) ! BANDGAP 

Tn = Tno/( 1 +Na/B1 ) {LIFETIME 

Tp=Tpo/< 1 +Nd/B1 ) 

Mun=FNMun(Na,T.AbO,Ab1 ,Ab2,Ab3> {MOBILITIES 

Mup=FNMup(Nd , T , BaO .Bal ,Ba2,Ba3) 

Dn=Kbtg*Mun {DIFFUSION COEFFICIENTS 

Dp=Kbtg*Mup 

Phio=Wo*Lambda*5 . 0330E+21 ! PHOTONS/CM2-S 

Ni=NiaO*Nia1 *(T~1 . 5>*EXP < -Egt/ ( Kbtg*2> > {INTRINSIC CARRIER DENS 

Vbi=Kbtg*LOG(Na*Nd/(Ni*Ni>> {BUILT IN VOLTAGE 

W2= (2 . *Epso*(8 . 85E-1 4) /Q)*Vb i*( 1 ./Nd+1 ./Na) {SQUARE OF DEPLETION W 
W=SQR( W2) {DEPLETION WIDTH 

Lp=SQR<Dp*Tp) {DIFFUSION LENGTH 

Ln=SQR(Dn*Tn) 

Xn=W*Na/(Na+Nd> {DISTANCE INTO N REGIO 

X P =W-Xn {DISTANCE INTO P REGIO 

U1 = (X j-Xp) /Ln 

U2=Dn/Ln 

U4=D P /Lp 

FI =FNRatio<U1 ,Sn,U2) 

U3=<B-Xj-Xn)/L P 

F2=FNRatxo<U3,S P ,U4) 

Pin=B*Wo*1 000 {POWER INTO CELL B by 1 CM (WAT 

lmax=0 

Kmax=60 

FOR J=1 TO Kmax 

X< J)=l .000+(J-1 >*.03 { ENERGY IN EV 

Xx < J) = 1 . 2402/X( J) {WAVELENGTH IN MICRON 

Phio=Wo*Xx( J>*5 . 033E+21 

Ener=X(J) 

IF- (Ener>l .7) THEN Xkk=Xkko 

IF (Ener<1 . 7 ) THEN Xkk=(Xkk 1*Ener-1 )*1 .OE-3 

Dei=Xkk*(T-TO> 

E=Ener+DeI 

Alpha=FNAbso(Nsi ,Esi(*> ,Asi<>) ,E) {ABSORPTION COEFFICIENT 

{PRINT "ALPHA =”; Alpha 
Z z I =A lpha*A 
Yyy 1 =0 

IF (Zzl <700) THEN Yyy t =EXP(-Zz1 > ! NO NUMERICAL INTEGRATION 

In=-Ln*F 1*Q*Phio*(1-Re)*(1-Yyy1 ) {CURRENT FROM RIGHT SIDE AMPS 
Ip=-Lp*F2*Q*Phio*( 1 -Re>*( 1 - Yyy 1 > {CURRENT FROM LEFT SIDE AMPS 
Isr=< In+Ip)/Pin 
IF < Isr> Imax ) - THEN Imax = Isr 

{PRINT In, Ip {OPTIONAL PRINT STATEMENT 

Aa( J , 1 )=In/Pin {SPECTRAL RESPONSE FROM RIGHT. 

Aa< J,2)=Ip/Pin {LEFT 

Aa( J , 3) =Isr {TOTAL RESPONSE 

NEXT J 

t *•*••**■■*■■*•***-**••*•■*•**■*■■*•«*•»■*■*•»■*■■*•»**•*■*■*•■*•■*■**■■*-*»•*■*•**■*■**■*-*■■**•*•*•*■*■* **■*■*■ *■*■*■*•■*•■*••*■ 


CONSTANT*TEMP 


{MOBILITIES 

{DIFFUSION COEFFICIENTS 

{PHOTONS/CM2-S 
{INTRINSIC CARRIER DENSITY 
{BUILT IN VOLTAGE 
{SQUARE OF DEPLETION WIDTH 
{DEPLETION WIDTH 
{DIFFUSION LENGTH 


{DISTANCE 

{DISTANCE 


INTO N REGION 
INTO P REGION 


(WATTS) 


{WAVELENGTH IN MICRONS 


! NO NUMERICAL INTEGRATION 
{CURRENT FROM RIGHT SIDE AMPS 
{CURRENT FROM LEFT SIDE AMPS 


J 

? PLOT SPECTRAL RESPONSE CURVES 

j 

GCLEAR 

WINDOW -.25.1 .3, -.2. 1 .1 
FRAME 

f *■*•*•■*■*•*■**•*•■*■*••*■)* *■♦■*■*•»*■*■*•**■»•*■*■*■**■*•*■*•■*■*••»■* ■**•*•***••*■*■•* »■*■*■*■■*■»■*■*■* * ■*■*** 
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1810 ! DRAW AND LABEL AXES 

1 820 ! **»*»•**»*■»*»****■»*******»***■*■**■******■**■*•*'***•**•*■»*•****■***■**•**'*■■*'**'*■*•**»•* 
1830 AXES .1 , .1 ,0. .0,5, 5, 6 
1840 CSIZE 4 
1850 LORG 6 

1860 FOR X_lafael= . 5 TO 1.5 STEP .5 

1870 MOVE X_label ,-.0125 

1880 LABEL X_la bel 

1890 NEXT X_labe 1 

1900 LORG 8 

1910 FOR Y_label=0 TO 1. STEP ,1 
1920 MOVE (-.005) ,Y_label 
1930 LABEL Y_label 
1940 NEXT Y_iabe 1 

1 950 * #»*#*#*#»#**##*»»»**#####*»*#»***♦»****»**»«#**»*#*»*»«#***#*##**#** 
1960 ! PUT LABELS ON THE AXES 

1 97 0 * it*^***********************#******^************************^*********^ 

1980 LORG 6 

1990 DEG 

2000 LDIR 0 

2010 MOVE .75,-. 1 

2020 LABEL "WAVELENGTH ( M ICRONS ) " 

2030 LDIR 90 
2040 MOVE -.19, .5 

2050 LABEL "NORMALIZED SPECTRAL RESPONSE" 

2060 ! *•***-*•******■**-***»*■*•*****■*****■**-*****■»*♦•***•***-**•**■*•*•**-***•****•*-*•**■*■*•*** 
2070 ! 

2080 ! DRAW THE SPECTRAL RESPONSE CURVES 

2090 ! 

2100 ’ ***■*■*•*■***■**•**■***■*■*•**■!**•»***'»'*-***'**■»■***•»*•*■*■»*•»***•*** •»*■***<*■**-***■»-■*■•**■■*••*••* 

2110 FOR 1=1 TO 3 

2120 FOR J-1 TO Kmax 

2130 Z< J)=Aa< J , I) /Imax 

2140 NEXT J 

2150 MOVE Xx < 1 ) , Z< 1 > 

2160 FOR K= 1 TO Kmax 
2170 DRAW Xx <K) , Z(K> 

2130 NEXT K 

2190 NEXT I 

2200 LDIR 0 

2210 MOVE -.15, -.05 

2220 LABEL "Si" 

2230 ! ************** ************ #***»«***#»*##«»***#*#*»##**#*»**»»#***»*##* 
2240 ? 

2250 ! PAUSE 

2260 ! TYPE DUMP GRAPHICS FOR HARD COPY . 

2270 1 ********************************************************************* 
2280 ! 

2290 PAUSE 
2300 ! 

231 0 ! ******************************************************************** 
2320 LDIR 0 

2330 GCLEAR !PUT A LABEL ON THE SCREEN 

2340 WINDOW 0,100,0,100 

2350 LORG 5 

2360 MOVE 50,25 

2370 LABEL "CURRENT-VOLTAGE" 

2380 MOVE 50,20 

2390 LABEL "CURVE IS NEXT" 

2400 MOVE 50,15 

2410 LABEL "CALCULATIONS TAKE 5 MINUTES " 
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2420! 

2430! 

2440 

2450! 

2460! 

2470! 

2480 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

2560 

2570 

2580 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

2670 

2680 

2690 

2700 

2710 

2720 

2730 

2740 

2750 

2760 

2770 

2780 

2790 

2800 

2810 

2820 

2830 

2840 

2850 

2860 

2870 

2880 

2890 

2900 

2910 

2920 

2930 

2940 


*** *-**•* ■*•*-* »■**-* •»■*■-*•* ***■**■*•*•*■»**•»*•*■*•*****■*•!*•*■•*••*■■*■**•*••*■•**■■*■■***■*■*■***•*■*•*• *•*■***•*■*■* *■•*■*■** 


PRINT " Isc 


Io 


Irec 


Volts" ! OPTIONAL 


* *-*-*>•*- * *■* * ■*■**--* -* 4* -*-*-*• •**-*•*-* ■*•*■■*■■*■*■■*■•**••* «*****•*•***-<*'***<*•* ■*•*■**■»■*• •**•*•**•■* •**■■*’■*■■*> -*-*■*•■*■*• 

TEMPERATURE DISTRIBUTION FROM TOP TO BOTTOM AS FUNCTION OF Z IS CONSTANT 
Mun=FNMun(Na,T,AbO.Ab1 ,Ab2,Ab3> MOBILITIES AS FUNCTION OF T 
Mup=FNMup < Nd , T , BaO . Bal , Ba2 , Ba3 > 

Kbtg=Kb*T {BOLTZMANN'S CONSTANT TIMES TEMPERATURE 

Egt=Eggo- < A Iphao )*T*T/ < T+Be tao ) !BANDGAP IN EV 

Dn=Kbtg*Mun {DIFFUSION COEFFICIENT <CK**2)/S 
Dp=Kbtg*Mup 

Ni=<Nia0>*<Nia1>*<T''1 .5)*EXP<-Egt/<Kbtg*2)> {INTRINSIC DENSITY 
Vbi=Kb*T*LOG(Na*Nd/(Ni*Ni ) ) !BUILT IN VOLTAGE (VOLTS) 

Lp=SQR<Dp*Taup) {DIFFUSION LENGTHS (CM) 


Ln=SQR(Dn 
U2=Dn/Ln 
U4=Dp/Lp 
Ppo=Na 
Nno=Nd ! 

Npo=Ni*Ni/Ppo 
Pno=Ni*N i /Nno 
U7=S P /U4 
U9=Sn/U2 


*Taun) 

!CM/S 
! CM/S 

! DOPING DENSITY 


CM-3 

CM-3 

CM-3 

CM-? 


NO DIMENSIONS 
NO DIMENSIONS 


f *»***#•*>***##»*##****•*#*#*«#■»**»***#****•»*»*#****#***»**!**»*#*»■»********* 

! CALCULATE CURRENT VOLTAGE RELATIONSHIP 

f 

f *****************************************************-****-i**-*‘****-**-*'*-*** > *'*fr 

Pmax=0 

Maxy=0 

Phi0=Wo*Lambda*5. 03306E+21 ! PHOT0NS/CM2-S 

Engo=1 . 2402/Lambda ! ENERGY IN EV 

IF ( Engo<- 1.7) THEN Xkk =( 1 . 3-Engo- 1 . >* 1 . OE-3 ! SEE ABOVE FOR Xkk 
Dele=Xkk*(T-T0) ! SHIFT IN BANDGAP DUE TO TEMP 

E=Engo+Dele ! EV 

AIpha=FNAbso(Nsi ,Esi(*> ,Asi(*> ,E) ! ABSORPTION COEFFICIENT CM- I 
Start: V=V+Delvo ! (VOLTS) 

Icount = IcouLnt + 1 
Vol ts( Icount)=V-*1 .QE+3 
IF ( Vo 1 ts< Icoan t ) >Max ) 

1 i = ° 

Again: Ii=Ii+1 

Z2=< Ii-1 )*Del z 
Z1 =Z2+De 1 z/2 

IF (Z1>A) THEM GOTO Over 
W2-(2*Epso*(8.85E-14)/Q)*<Vbi-V)»< 1 
W=SQR(W2) 

Xn=W*Na/(Na+Nd) 

X P =W-Xn 
U1 =(X j-Xp)/Ln 
FI =FNRat io(UV,Sn,U2) 

U3=<B-X j-Xn)/Lp 


! MILLIVOLTS 


THEN Max=Vo 1 ts ( Icount ) 


CM 

CM 


/Nd+1 ./Na) 


CM**2 

CM 


! CM 
! CM 

! NO . DIMENSIONS 
! NO DIMENSIONS 
! NO DIMENSIONS 


2950 

F2-FNRatio(U3.Sp,U4> 

! NO DIMENSIONS 

2960 

In=-Ln*F1 

! CM 

2970 

Ip=-Lp*F2 

!CM 

2980 

! SHORT CIRCUIT CURRENT 


2990 

Yy 1 =Aipha*Z1 

!NO DIMENSIONS 

3000 

Yy2 = 0 


3010 

IF ( Yy 1 <700 ) THEN Yy2=EXP< -Yy 1 ) 
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! AMPS/CM2 


3020 

3030 

3040 

3050 

3060 

3070 

3080 

3090 

3110 

3120 

3129 

3130 
3140 
3150 
3160 
3170 
3180 
3190 
3200 
3210 
3220 
3230 
3240 
3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
'3420 
3430 
3440 
3450 
3460 
3470 
3480 
3490 
3500 
35 1 0 
3520 
3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
3610 


Isc=(In+Ip)*-Phi0*(1-Re)*Alpha*Yy2*Q 
Y8=FNA < V , Ni , W , Kb , T , Vb i , T au.n , T aup ) 

Irec=Y8 ! RECOMBINATION CURRENT AMPS/CM2 

t *******#**»******«*«*»****«*****«*•*#*•**«'***'»*«*»«********#«****** 
f 

IPRINT Y8;V;Ni;W;Kb;T;Vbi;Taun;Taup ! OPTIONAL PRINT STATEMENT 
| ******************************************************************** 


R32=FNCur(U3,U7) 
Ipd=R32*Pno*Q*U4 
U8=<Xj-Xp)/Ln 
R33=FNCur(U8 ,U9) 
Ind=R33*Npo*Q*U2 
Io=? ( Ipd+Ind) 

Y isc< I i ) =Isc 
Yidark <Ii) = Io 
Y rec ( I i ) = 1 rec 
GOTO Again 


! NO DIMENSIONS 
! AMPS/CM**2 
! NO DIMENSIONS 
! NO DIMENSIONS 
! AMPS/CM**2 
! AMPS/CM2 
! AMPS/CM2 
? AMPS/CM2 
’ AMPS/CM2 


| *»********«*-»#********»»***iHt*************iH(*iHHH»»**'»****## + ***'* 

f 


! NUMERICAL INTEGRATION OF PREVIOUS RESULTS 


! 


Over: lrec=0 1AMPS/CM2 

lsc=0 ’AMPS/CM2 

I o=0 ! AMPS/CM2 

Jki-Ii-2 

FOR Ii-1 TO Jki STEP 2 

Isc=Isc+(Del z/3 .)*<Yisc(I i )+4*Y isc( Ii+1)+Yisc(Ii+2>) 
Io=Io+(Delz/3.)*(Y idark < I i )+4*Y idark < I i + 1 > + Y idark < I i +2) > 
Irec=Irec+<Del z/3.)*(Yrec(Ii) +4*Yrec< Ii + 1 )+Yrec(Ii+2)) 

NEXT Ii 

’ WE INTEGRATE OVER THE HEIGHT OF CELL— UNITS ARE NOW AMPS FOR ALL I'S 


U6=V/(Kb*T) ! NO DIMENSIONS 

j *•■*-**■■#■■*•*■*•*-*•**-*• *•*■*•*■*■•*■■*■»•**-*■*■**■■*■■*•*■*■-»■*■*•*•*■■**■*■**■ *-#•**■**■*■*■•**•*■«■*•■* «-**•*■■«•■*•■*■***■■■«■*■** 

t 

'PRINT Isc ; I o ; Irec : V ! OPTIONAL PRINT STATEMENT 

» 

t *♦*•■**■*■■*••*****#■*■**•**•■*■*•■»•»■■*■*■*■*•■*■*■*■*■*■*•■*•*■*■*•**•*«■*■■»■*■»«••*■•»•*■ #■*■■*■***■#■*■****#**■*****■*■**** 

IF (Rs=0) THEN GOTO Skip ! TEST TO SEE IF SERIES RESISTANCE IS ZERO 
I=FNNI < Rs , V , T , Io , Isc , Irec ,Kb ,Rsh) ! NONLINEAR CALCULATION OF CURRENT 
GOTO Ready 

Skip: I=Isc-Io*<EXP(U6)-1)-Irec-V/<Rsh) ’ TOTAL CURRENT AMPS 

Ready : ! 

Amps ( Icount > = 1* 1000 ’SAVE FOR PLOTTING CONVERT TO MILLIAMPS 
P=I*V ’POWER (WATTS) 

f i(*)(********iW*******iH(»****************)(***Jtif*i(******it**iti(***iHt*»*** 

’ FIND MAXIMUM POWER IN WATTS 

IF <P>Pmax> THEN Pmax=P 

IF < Amps ( Icount ) >Maxy ) THEN Maxy=Amps ( Icount ) 

IF <I<0) OR <Vb i<V+Delvo) THEM GOTO Finish 
IF < V>= . 45) THEN Delvo=Xstep 
GOTO Start 

Finish: Amps(Icount)=0 

f *■* ■*■*■*■*■»■■*■>*****■**'*•*•*•*■»■*■*•-*•■*•»•*■*■*■***■*•*••*■*■>*■*■>«■*■*■■*•*•*•«■*•*•**■*•*•*•*•*■»•*■» It****-*-*****-** 

t 

’NOW DO SOME GRAPHICS 

f 

J >*■*•*■*»*■*■!*■*•■»■■*■*■**■■**■*•*••»■*■**•*■**•* **-»■*■**•**■**•*■***■*■**■*••■*•*•*!•■■*■•*• t**#*#*** •»*■*•**■*»■*■■*•*■•* 
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3620 

3630 

3640 

3650 

3660 

3670 

3680 

3690 

3700 

3710 

3720 

3730 

3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3860 

3370 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

4030 

4040 

4050 

4060 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 

4210 


GCLEAR 

DUMP DEVICE IS 701 ! BE SURE TO TURN THE PRINTING DEVICE ON 

Yl=.2*Maxy 

Ddy =Maxy / 1 0 

WINDOW -< 150) ,(1200) ,-<Yl> .Maxy+Yl 

FRAME 

DEG 

AXES 10, (Ddy) ,0.0,10,2,6 

t ********•**■**•****♦■**•*•»*******■*•**■*■»■»*•»•***•**■**■*■*■*•***■*•*■***•*■*••*••*■*■*••**•*•*•**'*•*•*■ 

1 LABEL AXES X-AXIS FIRST 

f *•**■*■**•*■**•*■*•■*■***•*••***-»■-*****•*■■**•***■■*■**■*■****■**•*•* 

LORG C 
LDIR -90 

FOR X_labei=0 TO 1200 STEP 100 
MOVE X_label , -7*Y1/1 6 
LABEL X_label 
NEXT X_label 

f **■*•**•*•*•*•*•**•■**•*•* *■***■**•**•**»**■***■***-***■***********•** **■*■*■***•* **■*■»* *•*•**•* 

! NOW LABEL Y-AXIS 

1 *•*•**•*•*•*■*•*■**•■*■»■ *****-*•)«••*■*•'*■*■**■*•*■*•**■*■*■*■*••*■*■*****■*■*•*•*■•*■*■*■*■*****•■***• **■■* **■* ■»•»■»* 
LORG 8 
LDIR 0 

FOR Y_label=0 TO Maxy STEP Ddy 

MOVE 0 , Y_label 

LABEL DROUND < Y_1 abe 1 , 2 ) 

NEXT Y_label 
LORG 5 

t *****************»*********•*•*•*•*•*•*•*•*•*•*■*•*•*•*■***■*■*•***•**■•*•*■**■*■*■*•»■■»■*•***•**■*** 


♦LABEL GRAPH 
| 

MOVE (350) , -< Y 1*1 3/1 6) 

LDIR 0 

LABEL "MILLIVOLTS" 

MOVE -< 130) , (Maxy/2) 

LDIR 90 

LABEL "MILLIAMPS" 

CSIZE 4 
LDIR 0 

MOVE <-75>,-< 13*Y1/16) 

LABEL "Wo=" ; Wo 

f *********■■*■*•**■*•*****-*•-*■*■■***•*■**•***•*-***■■*■•*■-***■-**■**■*■*■#•*■*■**•*••*■-*-*•■*•**■****•#■ 


DRAW THE CURVE 


f * *-*■*■* -*•*•*-*** *>-*-**•*--*•* ■* ■#■■#*■-*-*■*■■*-**■-*■*■■*■*■-*■*■ -*■*-*■**■ •*■•*•* ■*■**'** -*•*■■*■*■’■*•*-*■'*■ 

MOVE Vo 1 ts < 1 ) , Amps < 1 ) 

FOR I j = 1 TO Icount 
DRAW Vo 1 ts< I j) , Amps< I j ) 

NEXT Ij 

MOVE 600 , 1 . 05*Maxy 

LABEL "SILICON VERTICAL JUNCTION CELL" 

MOVE 600, , 9*Maxy 
LABEL "Rs=" : Rs 
MOVE 600 , . 85*Maxy 
LABEL "T ="; DROUND <T, 4) 

MOVE 600 . . 8*Maxy 
LABEL "Sn=" ; Sn 
MOVE 600, . 75*Maxy 
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4220 LABEL "Sp=";Sp 
4230 MOVE 500, .7*Maxy 
4240 LABEL "Pout"" ; Pmax 
4250 MOVE 500, .65»Maxy 

4260 LABEL "Eff = " ;DROUND<Pmax*100/Pin,4> 

4270 MOVE 500 , . 6*Maxy 
4280 LABEL ”L£unbda=" ; Lambda 
4290 MOVE 500 , . 55*Maxy 

4300 LABEL "Jsc=" ; DROUND< Isc/8 ,2) AMP/CM2” 

4310 MOVE 500 , . 5*Maxy 

4320 LABEL "Isc=";Isc 

4330 MOVE 400 , . 45*Maxy 

4340 LABEL "WIDTH=" ; B ; " THICKNESS=" ; A 

4350 MOVE 500 , . 4*Maxy 

4360 LABEL "Xj=";Xj;" Voc=" ;Volts<Icount> 

4370 MOVE 500 , . 35*Maxy 

4380. LABEL "Ln=" ; DROUND<Ln , 5> Lp='' ; DROUND < Lp , 5 > 

4390 MOVE 500, .3*Maxy 

4400 LABEL "Mun=" ;DROUND<Mu.n ,5) ; '* Miip=” ; DROUND (Mup, 5) 

4410 MOVE 500 , . 25*Maxy 
4420 LABEL "Na=";Na:" Nd=";Nd 

4430 ! *** -****■ ■#■**■*■ **■*-*■* *■**■*■ *■**■*■-* *********** ****** *■****-*■**-**■ ******-* -*■***-*-**-* 

4440 I 

4450 PAUSE ! TYPE ' DUMP GRAPHICS' IF YOU LIKE WHAT YOU SEE 

4460 ! 

447 0 ! *■*■**■!*-*■***** **-*******■*■* * : *- j ***+** : ****-**-**-* j *****-**-*-*-*****-*i t -*■**■"+**■*■***'****** 

4480 GCLEAR 

4490 WINDOW 0,100,0,100 
4500 LORG 5 
4510 MOVE 50,50 

4520 LABEL "FINISHED WITH PROGRAM" 

4530 FOR K= 1 TO 5000 1 PAUSE AND THEN CLEAR THE SCREEN 

4540 NEXT K 
4550 GCLEAR 
4560 end 

457 0 * *-»)Ht**-»****i(*********«****»************** ■***•*■*■**■****■*•*•»•■**■***■*■*■*■*■*■*■*■*■■*•■*■* 

4580 ! 

4590 ! FUNCTIONS USED IN THE PROGRAM 

4600 ? 

461 0 ! 

4620 def FNMim<N3 , T , AbO , Abl , Ab2 , Ab3 ) ! MOBILITY FUNCTION 

4630 B1 =AbJ *T*T/<N3'' .66667) 

4640 D=N3*B1 *( 1 . - . 5*B1 ) 

4650 XI =Ab0*<T" 1 ,5>/D 
4660 X2=Ab2*EXP<-Ab3*L0G<T) ) 

4670 R=1 ./X2+1 ./XI 

4680 RETURN 1./R 
4690 FNEND 

4700 1 *****-***»****-**************-»'»********* 4 *i***-»^********** i, ******** < * <t ******* 

4710 def FNMup< N4 , T , BaO , Bal ,Ba2,Ba3) ! MOBILITY FUNCTION 
4720 X I =Ba2*EXP<-Ba3*L0G(T> > 

4730 B1 =Ba1 *T*T/<N4' .66667) 

4740 D=N4*B1 *< 1 . - .5*B1 ) 

4750 X2=BaO*<T'1 .5>/D 

4760 R=1 ./X2+1 ./X! 

4770 RETURN 1 . /R 
4780 FNEND 

4790 f *•*■*•»•**-*■*■*■*••*•»■***•*■*•*•*••*■«■■**■*■■*•**•»■»*•■*•#■» *•*•*■*•*■■* i****-*******************-*** 

4800 def FNSinh(X) JHYPERBOLIC SINE 

4810 RETURN <EXP<X)-EXP< -X> >/2 

4820 FNEND „ Q 



1830 * **■*•-**•■*■*■*»•*•**■*•*•*'»*•**•*■*•*■**■**•»■*•*•■*■•*•»■*■*■*•*•*•*■*••»*•*'*■***•*•*•*•**••*•*■**• *■»***■>* ******** ******* 

1840 def FNCosh(X) !HYPERB0L0IC COSINE 

1850 RETURN <EXP<X) +EXP< -X) )/2 
1860 FNEND 

187 0 * *•»•»■»****•*»■*»**•**■**•*»»*•»****•)»•»***»»'*******•»***»******»*•**»»*■*'**»*•****•»■***** 

1880 def FNRatio<X1 ,S,X2) • 'SPECIAL FUNCTION 

1890 IF <X1 >230) THEN GOTO Alternative 

1900 F1n=(S/FNCosh(X1 >>-X2*FNTanh<X1 >-S 

1910 F 1 d=X2+S*FNTanh(X1 ) 

1920 FI =F1 n/F 1 d 
1930 GOTO Next 
1940 Alternative: F 1 = — 1 
1950 Next:. ! 

1960 RETURN FI 
1970 FNEND 

1980 ' -t-******-*****-**********-**********************-*-*-**-***-*********"*-*-*-***-*** 

1990 def FNTanh<X) 1HYPERB0LIC TANGENT 

>000 IF <X1 >230) THEN GOTO Next 

>010 FI =( 1 . -EXP( -2 . *X) > /< 1 . +EXP<-2 . *X) ) 

>020 GOTO Ret 

>030 Next: FI = 1 

>040 Ret: RETURN FI 

>050 FNEND 

>060 ' ** **************************************** it*************************** 

>070 def FNCur(A.B) ! ANOTHER RATIO TYPE FUNCTION 

>080 FI =<B+FNTanh<A> )/ <B*FNTanh<A)+1 . ) ' 

>090 RETURN FI 
>100 FNEND 

>110 ! I****--*-*-****-*-*-*-*-**-****-*-***-*-***-*-***-****-*-*-*-*--*-**-***-*--*-*-*-**-*-***-*-*-*-*-*-*-*-******* 

>120 def FNN1 <Rs , V , T , Io , Isc , Irec , Kb , Rsh) ? NONLINEAR EQUATION SOLVER 
>130 Ib=Isc-Irec-V/Rsh+Io 'INITIAL GUESS AT CURRENT 


>1 40 Kbtg=Kb*T 
>150 \J ~ . 9999 

>160 Startl : U1=1+Rs/Rsh 


>170 

>180 

>190 

>200 

>210 

>220 

>230 

>240 


F=Ib*U*U1 -Ib+Io*EXP< <V+U*Ib*Rs)/Kbtg) 

Fp=+Ib*U1 +Io*EXP < < V+Ib*U*Rs) /Kbtg)*(Rs*Ib/Kbtg> 

U1=U-F/F P ! IMPROVED GUESS BY NEWTONS METHOD 

Error=ABS(U1-U) ! ITERATE UNTIL ERROR IS LESS THAN E-S 

IF (Error<1 .0E-6) THEN GOTO End 
U=U1 ’UPDATE OLD GUESS 

GOTO Startl ! DO EVERYTHING AGAIN 

End: RETURN U1*Ib 


>250 FNEND 


j 26 0 ' *’*****■*-******■**************•■*-*************************************** 
>270 'def FNA(X , Ni , W , Kb , T , Vbi , Taun , Taup) ! RECOMBINATION CURRENT 

>280 Y 1 = < 1 .6E-19)*PI*Ni*W*FNSinh<X/<2*Kb*T>> 

>290 Y2=SQR<Taup*Tau.n)*(Vbi-X)/ <Kb*T) 

>300 Y=Y1/Y2 

>31 0 RETURN Y 
>320 FNEND 


>330 ! ***•»•*••*•*•*•*'**$*'*■*■»•**'*■*■**■■*•■»■»■»■*•*•*■*■*■■*•*■*•*■*****♦•*"*•■*■»■*■* *•*•*■**■ ******* *■*■»*■*■**•*■* 

>340 def FNB (X,Ni , W , Kb , T , Vb i , T au.n , T aup ) ! DERIVATIVE OF ABOVE FUNCTION 

>350 Y 1 = ( 1 .6E-19)*PI*Ni*W*SQR<Taup*Taun>/(Kb*T) 

>360 Y2=C< Vbi-X)/<Kb*T) )*FNCosh<X/(Kb*T> ) 

>370 Y3=FNSinh(X/<Kb*T> ) 

>380 Y4 = Tau.p* Taun* < <Vbi-X)/(Kb*T) >"2 
390 Y=Y1 *< Y2 + Y3) /Y4 
>400 RETURN Y*1000 
410 FNEND 

420 def FNAbso (Nsi , Esi <*) , Asi < *) , Energ) 
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5430 

5440 

5450 

5460 

5470 

5480 

5490 

5500 

5510 

5520 

5530 

5540 

5550 


1 = 0 

IF (Energ<Es i ( 1 ) ) THEN Aipha=0 
IF (Energ<Esi < 1 ) ) THEN Stop3 
St3 : 1=1+1 

IF ( I>Ns i ) THEN Last3 

IF <Es i ( I ) <=Energ) AND (Esi (1+1 >>=Energ) THEN Ratio3 
GOTO St3 

Last3: Alpha=Asi(Nsi> 

GOTO Stop3 

Rat io3 : Alpha=Asi <I) +< Asi ( 1+1 )-Asi <I))*(Energ-Esi(I>>/(Esi(I+1)-Esi(I)) 
Stop3: ! 

RETURN Alpha 
FNEND 
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FIGURE Bl. Representative output from the program "VJSCP" 

illustrating spectral response and current-voltage 
relationship for vertical junction photoconverter. 
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APPENDIX C 

PROGRAM "SC3D" FOR SOLAR CELL 
3-DIMENSIONAL 
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PROGRAM SC3D(INPUT,OUTPUT,TAPE5=INPUT,TAPE6=OUTPUT, TAPES) 

C 

C SC3D=SOLAR CELL 3 DIMENSIONAL (ILLUMINATED FROM TOP Z=0) 

C 

C GREEN’S FUNCTION SOLUTION - SYMMETRIC IN X AND Y 
C 

COMMON / BLKl / NDATA,ENERG(25),ALP(25) 

COMMON/BLK2/POW.L AMBD A,N A,ND,B' < !,SN,SP 1 SO,Sl,RE 1 RS,RSH,A,B,H 

COMMON/BLK3/EPSO.A0 1 A1,A2,A3,BO,B1,B2,B3,EGGO,ALPHAO,BETAO 

COMMON/BLK'J/XKKO.MUN.MUP.TAUN.TAUP.DN.DP.PHIO.NI.LN.LP.VBI 

C0MM0N/BLK5/ETAN(200),XIN(200),XIP(200),ETAP(200) 

COMMON/BLKe/PI.PIS.Q.ZJ.T.TO.NIAl.NIAS.TNO.TPO.W, ALPHA 

COMMON/BLK7/NEIGV,KB,KBTG,XGRAPH 

COMMON/BLK8/VOC.PMAX.EFF.ISC 

DIMENSION ENERGS(25),ENERGGA(25),ENERGGE(25) 

DIMENSION ALPS(25), ALPG A(25),ALPGE(25) 


REAL LAMBDA,NA,ND,NIA1,NIA2,MUN,MUP,KB.KBTG,LN.LP,NI,ISC 
REAL NIAGE1,NIAGE2,NIAS1.NIAS2,NIAGA1,NIAGA2 

’ CALL PSEUDO 

NAMELIST VARIABLES 

NAMELIST /P AR AM/PO W,N A,ND,L AMBD A, A,B,H,Z J,T,RS,RSH,SN,SP, 
1 SO,Sl,NEIGV,ITYPE .ISWITCH 
C 

C POW = LASER POWER (KW/CM2) 

C NA=DOPING DENSITY (CM-3) 

C ND = DOPING DENSITY (CM-3) 

C LAMBDA=WAVELENGTH (MICRONS) 

C A = H ALF WIDTH OF X-DIRECTION C -A .LE. X .LE. A) 

C B = HALF WIDTH OF Y-DIRECTION C -B .LE. Y .LE. B) 

C H = DEPTH OF CELL IN Z-DIRECTION C 0 .LE. Z .LE. H) 

C ZJ=JUNCTION DEPTH FROM Z = 0 (CM) 

C T = TEMPERATURE (DEG K) 

C RS=SERIES RESISTANCE OHMS 
C RSH = SHUNT RESISTANCE OHMS 

C SNrSURFACE RECOMBINATION VELOCITY TOP (CM/S) 

C SP = SURFACE RECOMBINATION VELOCITY BOTTOM (CM/S) 

C SO = SURFACE RECOMBINATION VELOCITY TOP SIDES 
C S1 = SURF ACE RECOMBINATION VELOCITY BOTTOM SIDES 
C NEIGVzNUMBER OF EIGENVALUES 

C ITYPE = TYPE OF MATERIAL 

C ITYPE = 1 FOR SILICON (SI) 

C ITYPE = 2 FOR GALLIUM ARSENIDE (GAAS) 

C ITYPE = 3 FOR GERMANIUM (GE) 

C ISWITCH TURNS CONTOUR PLOT ON AND OFF 
C ISWITCH = 0 CONTOUR PLOT IS OFF 

C ISWITCH = 1 CONTOUR PLOT IS ON 

C 
C 
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NOMINAL VALUES FOR NAMELIST VARIABLES 


POWr.OOl 

LAMBDAsl.06 

T0=300. 

RS=0.0 

SN=1.0 

S0=1.0 

NA = 1.25E17 
A=.5 

H=5.00E-4 

NEIGVrlO 

ISWITCH=0 


$ RE = .05 

$ ITYPE = 1 
$ T = 300. 

$ RSH = 1.0 E 6 
$ SP = 1.0 
$ Sl = 1.0 
$ ND = 5.0E19 
$ B = .5 

$ ZJ = 2.5E-4 


SILICON PARAMETERS 
NDATAS=20 

DATA (ENERGSCK).Ks 1,8X1.1.1.15, 1.24, 1.378, l'.46,1.77,2.067, 2.48/ 

DATA CENERGSCK),K:9, 15X2.76,3. 1,3. 26.3.5,3.6.4.0. 4. 2/ 

DATA (ENERGSC K),K = 16, 20X5. 0,5.2, 6. 0,8. 0,10,0/ 

DATA CA LPSCKD.Ksl, 7X1. 0,10. 0,100. 0,350. O', 950.0, 2. E3.4.5E3/ 

DATA CALPSCK),Ks8,14X1.0E4,2.0E4,5.0E4,1.0E5,1.0E6,l.lE6,1.8E6/ 
DATA CALFSCK),K = 15,20X2.5E6,1.9E6.2.0E6,1.5E6,1.3E6,1.1E6/ 
CONSTANTS FOR MOBILITIES 

DATA ASO.AS1.AS2.AS3/65.02.5.72E-9.2.4 9E9.2.5/ 

DATA BSO,BSl,BS2,BS3/54.4 6,6.7 6 E - 9 ,2 . 93 E 9,2.7/ 

CONSTANTS FOR BANDGAP AS FUNCTION. OF TEMPERATURE 
DATA EGGOS.ALPH AOS, BET A0S/1.16.7.02E-4, 1108.0/ 

CONSTANTS FOR INTRINSIC CARRIER DENSITY 
DATA NIAS1.NIAS2/4.9E15.1.61/ 

CONSTANTS FOR LIFETIMES 

DATA TNOS,TPOS,B4/2.79E-4,3.52E-3,7.1E15/ 

EPSOS-11.6 

GALLIUM ARSENIDE PARAMETERS 
ABSORPTION COEFFICIENT CURVE FOR GAAS 
NDATAGA=16 

DATA (ENERGG AC K),K = 1,8X1.37, 1.38, 1.4,1.41,1.42,1.43,1.4 4,1.5/ 
DATA CENERGGACK),Ks9,16D / 1.6,1.8,2.0.2.4,3.0.3.5,4.0.4.5/ 

DATA C ALPGACKD,K = 1,8)/18..32.,119.,175.,2.5E3,4.5E3,6.E3,1.4E4/ 
DATA CALFGACK).K:9.16)/1.8E4.3.E4,4.E4.1.E5.2.E5,3.£5,4.£5.5.E5/ 

CONSTANTS FOR MOBILITIES 

DATA AGA0,AGA1,AGA2,AGA3/5076.11,5.24E-9,6.43E5.1.0/ 

DATA BGA0,BGA1,BGA2.BGA3/1477.3,8.59E-8,6.05E7,2.1/ 

CONSTANTS FOR BANDGAP 

DATA EGG0GA.ALPHGA.BETGA/1.519.5.4E-4.204./ 

CONSTANTS FOR INTRINSIC CARRIER DENSITY 
DATA NI AG Al.NI AG A2/4.9E15,. 0488844/ 

CONSTANTS FOR LIFETIMES 

DATA TNOG A.TPOGA.B 4/8. 8179 E -7,1.404 E-7.7.1E15/ 

EPS0GA = 13.1 
GERMANIUM PARAMETERS 
GERMANIUM SOLAR CELL PARAMETERS 
NDATAGE=22 
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DATA CENERGGECK),K = 1.8D/.66 I .6a,.7,.73,.76 1 .79,.a,.61/ 

DATA CENERGGECKD,K = 9.16)/.a3,.a5,.875,.9,.9 7,1.7,2.6,3./ 

DATA CENERGGECKD,K = 17,22D / 3.6,4.,4.4,5.0,6.0,8./ 

DATA CALPGECKD,K = l,8)/l.,10.,20.,6O.,10O.,6OO.,l.E3,^.E3/ 

DATA CALPGECKD,K = 9,16D / 5.E3,8.E3,8.5E3,9.E3,1.E4,1.E5 1 7.E5,8.E5/ 
DATA CALPGE(K),K = 17,22) / 1.E6,1.'4E6,2.0E6,1.3E6,1.2E6,1.0E6 / 

C 

DATA AGEO,AGEl t AGE2,AGE3/458.68,7.92E-9 1 6.08E7,1.6 6'' 

DATA BGE0,BGE1,BGE2,BGE3 / 3140.27,1.79E-7,1.18E9 1 2.33/ 

DATA EGGOGE,ALPHGE,BETGE'\741,4.77E-4,235./ 

DATA NIAGE1.NIAGE2/4.9E15,. 33313' 

DATA TN0GE,TP0GE,B^/2.1022E-5,2.1007E-5,7.1E15/ 

EPS0GE=16. 

C 

C READ IN ANY CHANGES TO THE ABOVE PARAMETERS 
C USE NAMELIST VARIABLES 

C 

C LABEL GRAPHICS OUTPUT AND PRINTED OUTPUT 

XGRAPHrO. 
lOO CONTINUE 

READ(5,PARAM} 

IFCE0F(5))600,601 

600 WRITEC6.603D 

603 FORMATC1X.28HEND OF FILE ENCOUNTERED-STOP ) 

STOP 1313 

601 CONTINUE 
C 

WRITEC6.H0) 

110 FORMAT(//) 

Q=1.6E-19 
XKKO=1.2E-3 
PI = 3. 1415926536 
PI2=PI*.5 
KB=8.6172E-5 
KBTG=KB*T 

C EPSO IS THE PERMITIVITY 

XGRAPH=XGRAPH+ 1 
C 

C SETUP PARAMETERS FOR CORRECT TYPE OF PHOTOCONVERTER MATERIALS 

C ITYPE = 1 CSI3 

C ITYPE = 2 CGAASD 

C ITYPE = 3 CUED 

IFCITYPE .EQ. ID WRITEC6.111D 
IFCITYPE .EQ. 2D WRITEC6.112D 
IFCITYPE .EQ. 3D WRITEC6.113D 

111 FORMATC1X.25HSILICON PHOTOCONVERTER , ' D 

112 FORMATC1X.35HG ALLIUM ARSENIDE PHOTOCONERTER .^D 

113 FORMATC1X.25HGERMANIUM PHOTOCONVERTER . ' D 

IFCITYPE .EQ. ID GO TO 10 
IFCITYPE .EQ. 2D GO TO 20 
C GERMANIUM DATA GOES HERE 

NDATA=ND AT AGE 
DO 40 I=1,NDATA 
ENERGCID=ENERGGE(I) 

40 ALPCID = ALPGECID 
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10 

c 


41 


20 

C 


42 


30 

C 

C 

C 

199 

200 

201 

202 

203 

204 


A0=AGE0 $ CzAGEl $ A2=AGE2 $ A3=AGE3 
B0 = BGE0 $ B1 = BGE1 $ B2 = BGE2 $ B3=BGE3 
EGG0 = EGG0GE $ ALPHA0=ALPHGE $ BETA0=BETGE 
NIA1 = NIAGE1 $ NIA2=NIAGE2 
TNOrTNOGE $ TPO=TPOG£ 

EPSO=EPSOGE 
GO TO 30 
CONTINUE 

SILICON PARAMETERS 
NDATA=NDATAS 
DO 41 Irl.NDATA 
ENERGCI) = ENERGSCI) 

ALPCO = ALPS(I) 

A0 = ASO $ A 1 = A SI $ A 2 = A S 2 $ A3 = AS3 

BO = BS0 $ B1 = BS1 $ B 2 = B S 2 $ B3 = BS3 

EGGO=EGGOS $ ALPHAO=ALPHAOS $ BETAO=BETAOS 
NIA1=NIAS1 $ NIA2=NIAS2 
TNOrTNOS $ TPO=TPOS 
EPSO=EPSOS 
GO TO 30 
CONTINUE 

GALLIUM ARSENIDE PARAMETERS 
NDATArNDATAGA 
DO 42 I=1,NDATA 
ENERGCI)=ENERGGA(I) 

ALPCI) = ALPGACn 

AO = AGAO $ A1 = AGA1 $ A2 = AGA 2 $ A 3 = A G A 3 

BOrBGAO $ B1 = BGA1 $ B2 = BG A2 $ B3 = BGA3 

EGGOsEGGOGA $ ALPHAO=ALPHGA $ BETAO=BETGA 
NIA1 = NIAGA1 $ NIA2 = NIAGA2 
TNOrTNOGA $ TPOrTPOGA 
EPSO=EPSOGA 
CONTINUE 

CALCULATE CONSTANTS 
CALL CONST 

CALCULATE SPECTRAL RESPONSE 
CALL SPEC 

CALCULATE I-V CURVE 
CALL IV(LAMBDA) 

WRITEC6.199DXGRAPH 
F0RMATCT20.6HGRAPH ,F3.0 ) 

WRITEC6,200)ISC,VOC,EFF,RSH 1 DN,NEIGV 
FORMAT( / ,T2,4HISC=,E12.5,T20 i 4HV0C = ,E12.5,T4O,4HEFF = , 

1 E12.5,T60,4HRSH = ,E12.5,T80,3HDN = ,E12.5,T99,8HNEIGV = ,15) 

WRITEC6,201)NA,ND,2J,Sl f DP.NIAl 
FORMAT (T 2,3HNA=,E12.5,T20 I 3HND = ,E12.5,T40 I 3HZJ = ,E12.5, 
1T60,3HS1=,E12.5,T80,3HDP = ,E12.5,T99,7HNIA1 = ,E12.5) 

WRITEC6,202)TAUN,TAUP,LN,LP,TN0,NIA2 
FORM ATCT2.5HT AUN = ,E12.5,T20,5HTAUP = ,E12.5,T40,3HLN = ,E12.5, 

1 T60,3HLP = ,E12.5,T80,6HTN0 = ,E12.5,T99,7HNI A2 = ,E12.5) 

WRITEC6,203)A,B,H,S0,TP0 

F0RMATCT2,2HA = I E12.5,T20,2HB:,E12.5,T40,2HH = ,£12.5,T6 0,3HS0 = , 
1 E12.5.T80.6HTP0 = ,E12.5) 

WRITEC6,204)MUN,MUP,LAMBDA,P0W 

F0RMATCT2,4HMUN=.E12.5,T20,4HMUP = ,E12.5,T40,7HLAMBD A-.F10.4, 
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lT60,iHP0W = ,E12.5 } 

WRITE(6,205DT,RS,SN,SP 

205 F0RMAT(T2,2HT = ,E12.5,T20,3HRSr,E12.5,T4 0,3HSN = ,E12.5, 

1 T60,3HSP = ,E12.5) 

WRITE(6,206)W,PMAX.NI,VBI, ALPHA 

206 F0RMATCT2,2HW=,E12.5 1 T20.5HPMAX = ,E12.5,T4 0,3HNI = ,E12.5, 

1 T60,4HVBI = ,E12.5,T80,6HALPHA = ,E12.5,///) 

C TEST TO SEE IF CONTOUR PLOT IS WANTED 

IF(ISWITCH.EQ.l) CALL CONTOUR 
GO TO 100 
END 

SUBROUTINE CONTOUR 
CONTOUR PLOTS OF CURRENT DENSITY 
VALUES STORED IN AAA AND BBB ARRAYS 
DIMENSION A A A(2Q,20),BBB(20,20) 

AAA IS 20X20 ARRAY OF CURRENT DENSITY VALUES AT X,Y AMD Z=ZJ 
BBB IS 20X20 ARRAY OF CURRENT DENSITY VALUES AT X,Y AND Z=ZJ+W 


N=20 

CALCULATE JP AND STORE IN AAA. FIND AMAX AND AMIN VALUES FOR ARRAY 
CALL cdenabcaaa.bbb.amax.amin.bmax.bmin.nd 
FINC rO $ NSET=0 $ NHI=-1 $ ND0T=0 

CALCULATE JN AND STORE IN BBB. FIND BMAX AND BMIN VALURS FOR ARRAY 


PLOT CONTOUR GRAPH OF AAA 

CALL CONRECCAAA.N.N.N.AMIN.AMAX.FINC.NSET.NHI.NDOT} 

CALL EZTRN 
CALL NFRAME 

PLOT CONTOUR GRAPH OF BBB ' 

CALL CONRECCBBB.N.N.N.BMIN.BMAX.FINC.NSET.NHI.NDOT) 

CALL EZTRN 
CALL NFRAME 
RETURN 
END 

SUBROUTINE CDENAB(AAA,BBB,AMAX,AMIN,BMAX,BMIN,N) 

CALCULATE CURRENT DENSITY AT X,Y COORDINATES SPECIFIED 
PUT RESULTS IN AAA AND BBB ARRAYS 
ALSO FIND MAX AND MIN VALUES OF ARRAY ELEMENTS 
DIMENSION AAAC20.20D.BBBC20.20) 

COMMON/BLK2/PO W,L AMBD A.N A.ND.B'I.SN.SP.SO.Sl.RE.RS.RSH.A.B.H 
COMMON / BLK4 / XKKO,MUN,MUP,T AUN.TAUP.DN.DP.PHIO.NI.LN.LP.VBI 
COMMON / BLK6 / PI,PI2,Q,Z J,T,TO,NI Al.NI A2,TN0,TP0,W,ALPHA 
C CALCULATE JPCX.Y) AND STORE IN AAA 

REAL NA.ND 

C CALCULATE JN(X,YD AND STORE IN BBB 

DELX = 2.a A^FLOATCN+l) 

AMIN=5.E5 
YWQ=ZJ +W 
AMAX=0.0 

DEL Y = 2.*B / FL0AT(N+1) 

BMIN=5. 0 E5 
BMAX=0.0 
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DO 10 J=1.N 
WRITEC6.99) 

99 F0RMATCT5,1HX,T15,1HY,T25,2HZ J,T43,2HJP.T55,1HX, 

1 T65,1HY,T75,4HZJ + W,T93,2HJN ) 

DO 20 I=1,N 

X=-A+IaDELX 
Y=-B+JaDELY 

CALL SUMIJCX.Y.XJP.XJN) 

WRITEC6,100)X,Y,ZJ,XJP,X,Y,YWQ,X JN 

100 F0RMATC3F10.4,E15.7,3X,3F10.6,E15.7 ) 

AAACI.J) = 4.0xuxALPHAaPHI0aC1-RE)aXJP 
BBBCI,JD = ^.0aQaALPHAaPHI0aC1-RE)aXJN 
IFCAAACLJ).GT.AMAX) AMAXrAAACLJD 
IFCAAACLJD.LT. AMIN) AMIN=AA ACL J) 

IFCBBBCLJD-GT.BMAXD BMAX = BBBCI. J) 

ifcbbbcljd.lt.bmind BMINsBBBCI.JD 

20 CONTINUE 
10 CONTINUE 
RETURN 
END 

SUBROUTINE SUMIJCX.Y.X JP,X JN) 

C SUM TERMS OF DOUBLE FOURIER SERIES FOR CURRENT 

C DENSITIES JN AND JP 

C 

COMMON / BLK2 / POW,LAMBDA,NA,ND,B'‘LSN,SP,SO,Sl,RE,RS,RSH,A,B,H 
COMMON/BLK^/XKKO.MUN.MUP.TAUN.TAUP.DN.DP.PHIO.NI.LN.LP.VBI 
C0MMON/BLK5/ETANC2OO),XINC2OO),XIPC2OO).ETAPC20O) 
COMMON/BLK6/PI,PI2, Q.ZJ.T, TO, NI Al.NI A2.TNO.TPO, W, ALPHA 
C0MM0N/BLK7/NEIGV,KB,KBTG,XGRAPH 
C INSERT COMMON BLKS 2.4, 5, 6,7 

REAL NA.ND,MUN,MUP,NI,LN,LP,NIA1,NIA2KB.KBTG 
C 

XJP=0.0 
X JN=0.0 

DO 10 I=1,NEIGV 
DO 20 J=1,NEIGV 

CALL TERNIJCX.Y.I.J.FPIJ.FNIJ) 

X JP=X JP+FFI J 
XJN=X JN+FNIJ 
20 CONTINUE 
10 CONTINUE 
RETURN 
END 

SUBROUTINE TERNIJCX.Y.I.J.FPIJ.FNIJ) 

C I.J TERM OF DOUBLE FOURIER SERIES 

C 

COMMON/BLK2/POW,LAMBDA,NA,ND,B4,SN,SP,SO,S1,RE,RS,RSH.A,B,H 
COMMON/BLK^/XKKO.MUN.MUP.TAUN.TAUP.DN.DP.PHIO.NI.LN.LP.VBI 
COMMON/BLK5/ETANC20O),XINC2OO).XIPC2OO),ETAPC20O) 
COMMON/BLK6/PI.PI2.Q.Z J.T.TO.NI Al.NI A2.TN0.TP0, W, ALPHA 
C0MM0N/BLK7/NEIGV.KB.KBTG.XGRAPH 
REAL NA,ND.LP,LN,NP,NN,MUN,MUP,NLNIA1,NIA2,KB,KBTG 
GCZ) = FUNC-ALPHAaZ) 

HHCAA,X)=CALPHAaCOSCAAaX)+AAaSINCAAaX))/CALPHAaALPHA+AAaAA) 
SCAA,X) = HHCAA.X)-GC2aX)aHHCAA,-X) 
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Y1CAA.X) = AAaSINHCX) + C0SHCXD 
YHCAA.B,X) = AAaCBaCOSHCXD + SINHCX)-BD 
C FPIJ (TOP REGION) SOLVE FOR JP 

UJPNOR = A+(SOaAaA/DP)aCGOS(XIP(J))/XIPCJ)Dxa2 

XMUI = ETAP(I)/B 

XL AM J = (XIP( J))/ A 

DEN2 = (XMUI a a2)+(XLAM Jaa2) 

GIJ.2=LPaLP/(1. + LPaLPaDEN2) 

GIJ=SQRT(GIJ2) 

VIPNORrB+CSOABAB/DP)ACCOS(ETAPCI))/ETAP(I')'lAA2 

NP=SPaGIJ/DP 

XXT=ZJ/GIJ 

SPI1 = BaSINCETAP(I))/ETAP(I) 

SPI2 = AaSINCXIP(J))/XIP( J) 

SPI3=COS(XLAMJaX)/UJPNOR 
SPH = COS(XMUIaY)/VIPNOR 
SPI6 = NPaSINHCXXT>COSH(XXTD 
ZO=ALPHAaGIJ-1 
XO=C1-EXP(-XXTaZO))/ZO 
X1=(1+(NP + 1)aXO) / SPI6 
X2=X1-EXP(-ALPHAaZJ) 

X3 = 1./(1 + ALPHAaGIJ) 

GMN=X3aX2 

FPIJ=GIJaSPHaSPI3aSPI2aSPI1aGMN 
C CALCULATE FNIJ 

HP = H-(ZJ+W) 

U JNNOR = A + (S1a A a A / DN)a(COS(XIN( J))/XIN( J))aa 2 
VINNOR = B + (Sl aBaB/DN)a(COS(ET AN(I))/ETAN(I))aa2 
XMUI=ETAN(I)/B 
XLAMJ=XIN(J)/A 
DEN3 = (XL AM J aa2)+(XMUIaa2) 

SIJ2sLNaLN^(1.+LNaLNaDEN3J " 

SIJ=SQRT(SIJ2) 

NN=SNaSIJ/DN 
XXB=HP/SIJ 

SNI1 = BaSIN(ETAN(I))/ETAN(I) 

SNI2 = AaSIN(XIN(J))/XIN( J) 

SNI3=COS(XLAMJaX)/UJNNOR 
SNHsCOSCXMUIaYKVINNOR 
SNI5 = YI4(SI J.NN.XXB) 

Z1=EXP(-ALPHAaH) 

Z2 = EXP(- ALPHAa(ZJ+W)) 

Z3 = ( ALPHA aSIJ)aa2 
Z4=1./CZ3-1) 

Z5=NN-ALPHAaSIJ 
Z6=NNaALPHAaSIJ-1 
SNI6=COSH(XXB>NNaSINH(XXB) 

HMN = ZiA(ZlAZ5 + Z2 aZ6aSINH(XXB)+Z2a(-Z5)aCOSH(XXB)) 
FNIJ=SIJaSNI1aSNI2aSNI3aSNI4aHMN/SNI6 
RETURN 
. END 

SUBROUTINE CONST 

CALCULATE CONSTANTS 

CALCULATES OTHER CONSTANTS AND PARAMETERS FOR PROBLEM 
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C0MM0N/BLK2/P0 W,L AMBD A.N A.ND.B-I.SN.SP.SO.Sl.RE.RS.RSH.A.B.H 

C0MM0N''BLK3''EPS0 1 A0 1 Ai 1 A2.A3,B0 1 Bl,B2 1 B3,EGG0 1 ALPHA0 1 BETA0 

COMMON/BLK^/XKKO.MUN.MUP.TAUN.TAUP.DN.DP.PHIO.NI.LN.LP.VBI 

COMMON/BLK6/PI,PI2,Q,2J.T,TO,NIA1,NIA2,TNO,TPO,W,ALPHA 

C0MH0N/BLK7/NEIGV,KB,KBTG,XGRAPH 


REAL MUN.MUP.KBTG.KB.LAMBD A,N A,ND 1 NIA1,NIA2,NI,LN,LP 


CONVERT WAVELENGTH (MICRONS) TO ENERGY (EV) 
ENER=1.2^02/LAMBDA 

CALCULATE BANDGAP AS FUNCTION OF TEMPERATURE 
EGT=EGGO- ALPHAO* TaT/(T>BETAO) 

CALCULATE ENERGY SHIFT DUE TO TEMPERATURE 
DELE=XKKO*(T-TO) 

E=ENER+DELE 

CALCUALTE ABSORPTION COEFFICIENT 
CALL ABSORP(E.ALPHAl) 

ALPHA=ALPHA1 
CALCULATE LIFETIMES 
TAUN:TNO / (l.+NA / B , 1) 

TAUP = TPO/(l.+ND/B^) 

CALCULATE MOBILITIES 
CALL MOBILCY1.Y2) 

MUN=Y1 $ MUP = Y2 
CALCULATE DIFFUSION COEFFICIENTS 
DNiKBTGaMUN $ DP = KBTGaMUP 
CALCUALTE DENSITY OF PHOTONS IN PH0T0N’S/CM2-S 
PHI0=P0WaLAMBDAa5.03306E21 ' 

CALCULATE INTRINSIC CARRIER DENSITY 
POWER = -EGT / (2.aKBTG) 

CALL EXPON(POWER.ANS) 

NI = NIA1aNIA2a(Taa1.5)aANS 
CALCULATE BUILT IN VOLTAGE 
VBI = KBTGaALOGCNAaND/CNIaNI)) 

VOLTSrO.O 

CALL WIDTHC VOLTS) 

CALCULATE DIFFUSION LENGTHS 

LP=SQRTCDP*TAUP) $ LN=SQRTCDNaTAUN) 

CALCULATE ALL NECESSARY EIGENVALUES FOR PROBLEM 
CALL ESET 
RETURN 
END 

SUBROUTINE ABSORP(E, ALPHA) 

CALCULATE THE ABSORPTION COEFFICIENT FOR GIVEN ENERGY E 


COMMON/BLKl/NDATA,ENERG(25),ALP(25) 


1=0 
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20 


IFCE.LT. ENERGC1DDG0 TO 10 
CONTINUE 
1 = 1+1 

IFCI.GT.NDATA)GO TO 30 

IFCCENERGCI).LE.ED.AND.(ENERGCI+1) .ge. eddgo to 40 
GO TO 20 

40 ALPHA=ALPCI)+CALPCI+1)-ALPCI))aCE-ENERG(I))/CENERGCI + 1)-ENERGCI^ 
RETURN 

10 A L P H A = 0 . 

RETURN 

30 ALPHA: ALPCNDATA) 

RETURN 

END 

SUBROUTINE EXPONCX.ANS) 

C EXPONENTIAL TERM IN CASE ARGUMENTS GET LARGE 

IF(X.LT.-500) GO TO 10 
IFCX .GT. 700)GO TO 20 
ANS=EXP(X) 

RETURN 
10 A N S = O . O 
RETURN 

20 WRITEC6.30) 

30 F0RMATC1X.2SH ARG TOO LARGE, EXP FUNCTION ) 

STOP 444 
END 

SUBROUTINE ESET 

CALCULATE EIGENVALUES FOR ORTHOGONAL FUNCTIONS 


COMMON/BLK2/POW l LAMBDA.NA,ND,B4,SN,SP,SO,Sl,RE,RS,RSH,A,B,H 
COMMON/BLK4/XKKO,MUN 1 MUP,TAUN,TAUP,DN,DP,PHIO,NI,LN.LP,VBI 
COMMON/BLK5/ETANC200D,XINC200),XIPC200),ETAPC200) 

COMMON/BLK6/PI, PI2, Q.ZJ.T, TO, NI Al.NI A2.TNO.TPO, W. ALPHA 
C0MM0N/BLK7/NEIGV,KB,KBTG.XGRAPH 

REAL MUN,MUP,NI,LN,LP,NIA1.NIA2,KB 1 KBTG,NA.ND, LAMBDA 

CALCULATE EIGENVALUES FOR TOP REGION X-DIRECTION 
R1 = S0 $ R2 = DP $ R3 = A $ N4 = NEIGV 
CALL EIGEN(XIF,R1,R2,R3,N4) 

CALCULATE EIGENVALUES FOR TOP REGION Y-DIRECTION 
R1=S0 $ R3=B 

CALL EIGENCETAP,Rl,R2,R3,N423nP CALCULATE EIGENVALUES FOR BOTTOM R 
R1 = S1 $ R2 = DN $ R3 = A 
CALL EIGEN(XIN,R1,R2,R3,N4) 

CALCULATE EIGENVALUES FOR BOJOM REGION 
R 1 = S 1 $ R 3 = B 

CALL EIGENCET AN,R1,R2,R3,N4) 

RETURN 

END 

SUBROUTINE EIGENCXX, S.D.H.N} 

FINDS EIGENVALUES FOR EQUATION TANCX) = C^X, C=CONSTANT 

C0MM0N/BLK6/PI,PI2,Q,ZJ,T,T0 1 NIA1,NIA2,TN0,TP0,W, ALPHA 
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REAL • NIA1.NIA2 


ORIGINAL 'PAGE- S3 
OF POOR QUALITY 


DIMENSION XXC200) 

F(2) = TAN([Z)-C / 2 

FPCZ) = C/CZxZ) + i./CCOS(2)AA2) 

C SOLVES EQUATION X*TANX=C 
C N MUST BE .LE. 200 
C=SaH^D 

IF(C.EQ.O) GO TO 1001 
C GUESS AT FIRST ROOT 
X1 = PI'M. 

M = 1 $ ICsl 

200 Y1=F(X1) 

IF(Y1.GT.0)G0 TO 10 
C Y 1 IS NEGATIVE 
X2=X1 
Y2= Y 1 

5 X1:.5a(X2 + (2aM-1)aPI2) 

Yl=F(Xl) 

IFC Y1.GT.0)G0 TO 20 
C Y 1 IS NEGATIVE 

X 2 - X 1 
Y2 = Y1 
GO TO 5 

10 X 2 = C M - 1 ) a P I * . 8 a C X 1 - C M - 1 ) a P I) 

Y2=FCX2) 

IFCY2.LT. 0)G0 TO 20 

XI = X2 
Y1=Y2 

GO TO 10 

2 O X -- C X 1 * X 2 ) / 2 . 

Y = F ( X ) 

I F (' A B 3 ( Y ) .LT. 1.0E-5)G0 TO LOO 
IF( Y .LT. O) GO TO 30 
X 1 T X 
Yl-Y 

GO TO 20 

3 0 X 2 = X 

Y2 = Y 

GO TO 20 
i 0 O X X ( I C ) - X 

M = M + 1 $ I C - I C 1 

JC1 = X*PI 

IFCIC.GT.N) GO TO 300 
GO TO 2 O 0 
300 CONTINUE 
RETURN 

i O C 1 WRITE (6,234 ) 

234 FORM ATCH5HREGOMBINA.TION VELOCITY OF ZERO NOT ALLOWED. 
END 

S TJBRO uti n e mobilcx M U N , X M TJ p ) 

C CALCULATE MOBILITIES IN N AND P REGIONS 
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COMMON/BLK2/POW,LAMBDA,NA,ND,B4,SN,SP.SO,Sl,RE.RS,RSH,A,B,H 
COMMON/BLK3/EPSO, AO, Al, A2, A3,BO,Bl,B2,B3,EGGO,ALPHAO,BETAO 
COMMON''BLK4''XKK0,MUN,MUP,TAUN,TAUP,DN,DP,PH2O,NI.LN,LP,VBI 
C0MM0N/BLK6/PI, PI2, Q,ZJ,T, TO, NIA1,NIA2,TN0,TP0,W, ALPHA 


REAL LAMBDA, NA,ND,MUN,MUP,NI.LN,LP.NIA1,NIA2 

MOBILITY OF ELECTRONS 
RQ=2./3. 

BB1=A1xTaT / (NAxxRQ) 

D=NAaBB1aC1.-.5aBB1) 

X1 = A0a(Txx1.5KD 
ZZ = - A3 a ALOGCTD 
CALL EXPONCZZ.ANSD 
X2=A2aANS 
R = l./X2 + iyxi 
XMUN=1./R 

MOBILITY OF HOLES 
ZZZ=-B3aAL0G(TD 
CALL EXPONCZZZ.ANS) 

X1=B2aANS 

BB1=B1aTaT^CNDaaRQ) 

D = NDaBB1aC1.-.5aBB1) 

X2 = BOx(Tax1.5} / D 

R = l./X2 tl./Xl 

XMUP=1./R 

RETURN 

END 

SUBROUTINE WIDTHCVOLTS) 

CALCULATE WIDTH OF DEPLETION REGION 


COMMON/BLK2/PO W.LAMBD A,N A.ND.B-I.SN.SP.SO.Sl.RE.RS.RSH.A.B.H 
COMMON/BLK3/EPSO,AO,A1,A2,A3,BO,B1,B2,B3,EGGO,ALPHAO,BETAO 
COMMON/BLKi/XKKO,MUN,MUP,TAUN,TAUP,DN,DP,PHIO,NI,LN,LP,VBI 
COMMON / BLK6 / PI,PI2,Q,Z J,T,TO,NI Al.NI A2.TN0.TP0, W, ALPHA 
COMMON/BLK7/NEIGV.KB.KBTG.XGRAPH 

REAL KB, KBTG, LAMBDA, LN,LP,MUN,MUP,NA.ND,NI,NIA1,NIA2 
ABC=C1.^ND ♦ l./NA) 

W2=(2. aEPSOa(8.85E-1'4) / Q)aC VBI- VOLTS)a ABC 
W = SQRT(W2) 

RETURN 

END 

SUBROUTINE SPEC 

CALCULATE SPECTRAL RESPONSE 

COMMON / BLK3 / EPSO,AO, Al,A2,A3,BO,Bl,B2,B3,EGGO,ALPHAO,BETAO 
CALCULATE SPECTRAL RESPONSE 


DIMENSION WAVC200),SPC200),SNC200),SDRC200),SRC200) 
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1=0 

9 0 1 = 1 + 1 

E=EGG0-.05 + .025*(I-l) 

IFCCE.GT.5.0D.0R.CI.GT.200)DG0 TO 100 

WAV(ID = 1. 2402/E 

CALL SH0RTCS1P.S1N.S1DR.E) 

SNCID=S1N 
SPCI)=S1F 
SDRCXD = S1DR 
3RCI)=S1N + S1P + S1DR 
WRITEC6,aO)I.E.WA VCID.S1P.S1N.S1DR 
80 FORM ATCI5.5E15.5) 

GO TO 90 
100 CONTINUE 
MAX=I-1 
1 = 1-1 

C NORMALIZE SPECTRAL RESPONSE 

XMAX=0. 

DO 110 K=1,MAX 

ifcsrck).gt.xmaxdxmax=srckd • 

110 CONTINUE 

DO 120 K=1,MAX 

SN(K)=SNCK)/XMAX 

SPCK)=SP(KXXMAX 

sdrckd=sdrckd/xmax 

SRCK)=SR(K)/XMAX 
120 CONTINUE 
C PLOT GRAPH 

NXL=20 $ SX=10. ' $ XOFF = . 5 

NYL=28 $ SY=10. $ Y0FF = .5 

I S Y M = O $ XMINrO. $ ' X M A X = 1 . 2 
IEC=0 

CALL GR APH(W A V,SN,SF,SDR,SR,I) 

RETURN 

END 

SUBROUTINE GRAFHCWA V,SN,SP,SDR,SR,I) 

C GRAPHICS OUTPUT OF PREVIOUS RESULTS FOR SPECTRAL RESPONSE 

DIMENSION W A V(200),SN(200),SP(200),SDR(200),SR(200) 

CALL CALPLTC1. 5,1. 0,-3) 

XPG=5.0 

XDV=10.0 

XTIC=.5 

YPG=5.0 

YDVrlO.O 

YTIC=.5 

CALL ASCALECWAV.XPG.I.l.XDV) 

CALL ASC ALECSR, YPG,I,1,YDV) 

CALL AXES(0.,0.,0.,XPG,0.0,WAV(I + 2),-XTIC,XDV, 

120HWA VELENGTH CMICRONS) ,.1-4, -20) 

CALL AXES(0.,0.,90.,YPG,0.0,.2,-YTIC,YDV, 

1 28HNORMALIZED SPECTRAL RESPONSE ,.1-4,28) 

SFX = 1./WA VCI + 2) 

SFY=5.0 

CALL CALPLT(WAV(1)*SFX,SNC1)*SFY,3) 
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DO 100 J=1.I 

CALL CALPLTCWAVCJ>SFX,SNCJ>SFY.2) 

100 CONTINUE 

CALL CALPLTCWAVC1)^SFX,SPC1D*SFY,3) 
DO 101 J = 1.I 

CALL CALPLTCWAVCJ)aSFX,SPCJ>SFY,2D 

101 CONTINUE 

CALL CALPLTCWAV(i;)*SFX.SDRCl)*SFY,3) 
DO 102 J = 1,I 

CALL CALPLTCWAVCJDaSFX,SDR(J)aSFY,2) 

102 CONTINUE 

CALL C ALPLT(W A V(1)*SFX,SR(1)*SFY,3) 
DO 103 J=1,I 

CALL CALPLTCWAVCJ>SFX,SRCJ>SFY,2]f 

103 CONTINUE 

CALL CALPLTCO.,6.,-3) 

RETURN 

END 

SUBROUTINE SHORTCSIP.SIN.SIDR.ED 
SHORT CIRCUIT CURRENT 


COMMON/BLK2/POW,LAMBDA,NA,ND,B4,SN,SP,SO,Sl.RE,RS,RSH,A.B,H 

common/blk^/xkko,mun,mup,taun,taup,dn,dp,phio,nlln,lp,vbi 
COMMON- / 'BLK 6 / PI,PI2,Q,2J,T 1 TO,NIAl,NIA2,TNO,TPO,W, ALPHA 
C0MM0N/BLK7/NEIGV,KB,KBTG,XGRAPH 

REAL L AMBD A,NA,ND,MUN,MUP,NI,LN,LP,KB,KBTG,NIA1,NIA2 

CALCULATES THE SPECTRAL RESPONSE COMPONENTS 
CALCULATES DOUBLE SUMMATION" OVER I.J 
CALL ABSORPCE.ALPl) 

ALPHA=ALP1 

S1P = 0 $ S1N = 0 

DO 10 Isl.NEIGV 

DO 20 J=1,NEIGV 

CALL NUMBSCI, J,SNI J.SPI J) 

S1P = S1P + SPI J 
S1N=S1N+SNIJ 
O CONTINUE 
O CONTINUE 

. CALL EXP0NC-ALPHA*ZJ,Q1) 

CALL EXPONC-ALPHAaW.02) 

S1DR-Q1a(1.-Q2)a‘5.a AaB 

RETURN 

END 

SUBROUTINE DCURENT(IO) 

DARK CURRENT 


C0MM0N/BLK2/P0W,LAMBDA,NA,ND,B^ 1 SN,SP,S0 1 S1,RE,RS,RSH,A,B,H 
COMMON / BLK^ / XKKO,MUN,MUP,T AUN.TAUP.DN.DP.PHIO.NI.LN.LP.VBI 
C0MMON/BLK5/ETANC20O),XINC20O),XIPC2OO).ETAP(20O) 
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C0MM0N / BLK6 / PI,PI2,Q,ZJ,T,T0,NIA1,NIA2,TN0 1 TP0,W, ALPHA 
C0MM0N/BLK7/NEIGV,KB,KBTG,XGRAPH 
C 

REAL NNO.NPO, LAMBDA, NA.ND'.MUN.MUP.NI.LN.LP.NIAl.NIAB 
REAL KB.KBTG.IO.NN.NP 
C 

Y1CX,ZD = ZaCOSHCX) + SINHCXD 
Y2CX.Z) = ZaSINHCX>COSHCX) 

Y3CX.Z) = CZ + TANHCX)KC1 + ZaTANH(X)) 

C DOPING DENSITIES 
PPO=NA 
NNO=ND 

NPO=NIaNI/PPO 
PNO=NI*NI/NNO 
SUMlrO $ SUM 2 = 0 
DO 10 Isl.NEIGV 
DO 20 J=1,NEIGV 

DEN2 = CETAPCJVB)aa2 +(XIPCID/A>a2 
GIJ2 = LPaLP/C1.aLPaLPaDEN2) 

GIJ=SQRTCGIJ2D 

DEN3 = CXINCID y/ ADAA2 + (ETAN(J)/B)aa2 
SI J2 = LNaLN / C1* + LNaLNaDEN3) 

SIJ=SQRTCSI J2D 

NP = SPaGIJ/DP $ NN = SNaSIJ/DN 

UJNNORM = A-*-CSOaAaA/DN>(COS(XINCJ))/XINCJ))aa2 

VINN0RM=B+CS1aBaB/DN)aCC0SCETANCI)KETANCIDDaa2 

F1 = SINCXINCJD) $ F2 = SIN(ETAN(I)) 

HP=H-ZJ-W 

F3 = 16aAaAaBaBaQaDNaNIaNI/CNAaSIJ) 

F^ = CF1/XINC J))aa 2 $ F5 = (F2/ETANCn)AA2 

F6=1./CUJNN0RMaVINN0RM) 

TERMN=F3aF^aF5aF6aY3CCHP/SIJ),NN) 

UJPN0RM = A + CS0aAaA/DP)aCC0SCXIPCJD) / XIPCJ)Daa2 
VIPN0RM = B + CS1aBaB/DPDaCC0S(ETAPCI)KETAPCI^aa2 
F1 = SINCXIPCJ)D $ F2 = SINCETAP(I)) 

G3 = 16aAaAaBaBaQaDPaNIaNI/CNDaGIJ) 

G-i = (Fl/XIP(J))AA2 $ G5=(F2 / ETAPCI))aa2 
G6 = 1./CUJPN0RMaVIPN0RM) 

TERMP = G3aG^aG5aG6aY3CCZJ/GIJ),NP) 

SUM1=SUM 1 + T E R M N 
SUM2=SUM2+TERMP 
20 CONTINUE 
10 CONTINUE 

I0 = SUM1+SUM2 
RETURN 
END 

SUBROUTINE NUMBSCI, J.SNI J.SFI J) 

I,J TH TERM OF DOUBLE FOURIER SERIES 


C0MM0N/BLK2/P0W.LAMBDA,NA,ND,B^ t SN,SP,S0,Sl,RE,RS,RSH,A,B t H 
COMMON/BLK^/XKKO,MUN,MUP,TAUN,TAUP,DN,DP,PHIO,NI,LN,LP,VBI 
COMMON/BLK5/ETANC2OO),XINC2OO),XIPC2O0),ETAPC20O) 
COMMON/BLK6/PI,PI2 ,Q,ZJ,T i TO,NIA1,NIA2,TNO,TPO,W, ALPHA 
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C0MM0N/BLK7/NEIGV.KB,KBTG,XGRAPH 

C 

REAL LAMBDA.N A,ND,MUN,MUP,NI,LN,LP,NIA1,NIA2,KB,KBTG,NN,NP 
EXTERNAL FUN 
C 

G(Z}=FUN(-ALPHAaZ3 

Y1CX,Z)=ZaC0SHCX)+SINHCXD 

YECX,Z}=Z*SINH(X>COSHCX) 


UJPNOR = A+CSOaAxA/DPDaCCOSCXIPCJ)KXIPCJ)Daa2 
DEN2=CETAPCID x B)aa2 + (XIPC JXADaaP 
GIJ2 = LPaLP''(1. + LPaLPaDEN2} 

GIJ = SQRTCGIJ2) 

VIPNOR=B+(SOaBaB / DP)a(COS(ET APCI)) / ET AP(I})aa2 
UJNNOR^A+CSIaAaA/DNDaCCOSCXINCJDD/XINCJDDaaB 

VINN0R = B-KS1*B*B/DND*CC0S(ETANCI)) / ETANCID)**2 

DEN3 = CXINCJ)/A)aa2 +CET AN(IKB)aa2 
SIJ2=LNaLN/C1.+LNaLNaDEN3) 

SIJ = SQRTCSI J2) 

NP=SPaGIJ/DP $ NN = SNaSIJ/DN 

TEST = CALPHAaGIJ-1) 

IFCABSCTEST3.LE..001DG0 TO 10 

SPI1 = 16a ALPHAaAaAaBaBaGI J/(U JPNORa VIPNORaC(ALPHAaGI J)aa2-1)) 
F4=SINCXIPCJ)) $ F5=SIN(ETAP(I)) 

SPI2 = CCF5/ETAPCI)]aa2)aCCF4/XIPC J))*»2) 

SPI3 = NP + ALPHAaGIJ-GCZJ)aY1CCZJ/GIJD,NP) 
SPM=CSPI3/Y2CCZJ/GIJ),NP))-ALPHAaGIJaGCZJ) 
SPIJ=SPI1aSPI2aSPM 
GO TO 40 

10 SPI1 = 16a ALPH Aa A a A aB a B aGI J^CU JPNORa VIPNORaCALPHAaGI J+1)) 

F4 = SINCXIPCJDD $ F5 = SIN(ETAP(I)) 

SPI2 = CCF5/ETAPCI))aa2)a(CF4/XIPCJ))aa2) 

BX=ZJ/GIJ 

SUM=BX $ XL=BX 

20 RATI0=(1 + CNP + 1)aSUM)/(NPaSINHCBX) + C0SHCBX)) 

SPI3 = RATI0-GCZJ) 

SPIJ=SPI1aSPI2aSPI3 
40 CONTINUE 
HP=H-ZJ-W 

G5 = SINCXINC JD) $ G6 = SINCET AN(I)) 

SNI1 = 16 a ALPH Aa A aAaBaBaSI JaGCCZ J+W)) 

SNI2=(G5/XINC J))aa2 
SNI3=CG6/ETANCI))aa2 
SNI4 = 1./'(UJNN0RaVINN0R) 

TEST=ALPHAaSIJ-1 
IFCABSCTEST).LE..001)G0 TO 50 
SNI5=1./CCALPHAaSIJ)aa2 -1.) 

XA = NNaCCOSHCHP/SIJ)-GCHPDD^SINHCHP/SIJD 
XB=ALPHAaSIJaGCHP) 

SNI6 = ALPHA aSIJ-(XA + XB) / Y2((HP''SIJ),NN3 

SNIJ=SNI1aSNI2aSNI3aSNI4aSNI5aSNI6 

RETURN 

50 CONTINUE 
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BY=HP/SIJ 

SUM=BY 

60 RATI0=C(NN-1.)*SUM - 1)xG(HP)''CNN*SINHCBY) + C0SH(BY)) 
SNI5=(1.+RATI0KCALPHAaSIJ + 1.) 
SNIJ=SNI1*SNI2*SNI3*SNI4*SNI5 
RETURN 
END 

FUNCTION FUN(X) 

IFCX.LT.-230)G0 TO 10 
FUN=EXP(X) 

RETURN 
10 F U N : O 

RETURN 

END 

SUBROUTINE IV(WAVE) 

C CURRENT- VOLTAGE CURVE 

C 

C CALCULATE CURRENT- VOLTAGE CURVE 

C 

C 

C0MM0N / BLK2 / P0 W,L AMBD A,N A,ND 1 B-‘!,SN,SP,S0,S1,RE,RS,RSH,A,B,H 

COMMON/BLK4/XKKO.MUN,MUP 1 TAUN,TAUP,DN,DP,PHIO,NI,LN,LP,VBI 

COMMON/BLK5/ETANC200D,XINC200D,XIPC200).ETAPC200) 

COMMON/BLK6/PI,PI2,Q,ZJ,T, TO, NIA1.NI A2.TN0.TP0.W. ALPHA 

C0MM0N/BLK7/NEIGV.KB.KBTG.XGRAPH 

COMMON/BLK8/VOC.PMAX.EFF.ISC 


DIMENSION AMPSCIOO).VOLTS(IOO) 

REAL IREC.IB.IO.ISC.I 

REAL L AMBD A,N A,ND,LN,LP,MUP,MUN p NI A1,NIA2,KB,KBTG,NI 
C CALCULATE IV CURVE 
FFF=Q*PHI0*(1-RE) 

E = l. 2402/WAVE 

ICOUNTrO $ PM AX=-.5 

DEL V0 = . 0125 

DELV=.025 

C CALCULATE ISC (LIGHT CURRENT) 

CALL SH0RT(S1P,S1N,S1DR,E) 

ISC=FFF*CS1P+S1N+S1DR) 

V=-DELV 
100 CONTINUE 

IC0UNT=IC0UNT+1 
V= V+DELV 
CALL WIDTH(V) 

C CALCULATE DARK CURRENT 
CALL DCURENTCIO) 

IREC = 4*A*B*Q*NIaW*PI2*2*SINHCV/C2*KBTG)) 
IREC=IREC/(CVBI-V)/KBTG) 

IREC=IREC/SQRT(TAUN*TAUP) 

C SOLVE FOR I 

U1=1+RS/RSH 

U=,999 

IB = ISC-IREC-V /RSH+IO 
V0LTS(IC0UNT)=V*1000 
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110 CONTINUE 

Z=( V + U*IB*RS} / KBTG 

F = IBaUaU1-IB + I0aEXPCZD 

FP = IBaU1 + I0aCRSaIB/KBTGDaEXPC2D 

U2=U-F/FP 

ERROR: ABSCU2-U) 

IFCERROR.LT.l.OE-6) GO TO 120 
U=U2 

GO TO 110 
120 I = U 2 a I B 

AMPSCIC0UNT)=Ia1000 

P=IaV 

IFCP.GT.PMAX)PMAX=P 

IFCCI.LT. 0).OR.CVBI.LT.(V+DELVD)DGO TO 130 
IFC V.GT..35) DELV=DELVO 
GO TO 100 

130 AMPSCIC0UNT) = 0. 

PIN=POWa1000a1aAaB 

EFF=PMAXa100/PIN 

V0LTSCIC0UNT)=V0LTSCIC0UNT-1>DELV/100. 

VOC=VOLTSCICOUNT) 

N = ICOUNT $ YMIN = 0 $ XMIN = 0 

YM AX = 1.2a AMPSC1D 

XMAX=1.2aVOLTSCICOUNTD 

NXL=20 

NYL=19 

ISYM = 0 $ S X = 5 . 0 $ SY = 5.0 

XOFF = 0.75 $ YOFF = 0.75 

CALL GRAPHI VC VOLTS, AMPS, ICOUNT) 

CALL NFRAME 

RETURN 

END 

SUBROUTINE GRAPHIVCX.Y.I) 

C GRAPHICS OF CURRENT-VOLTAGE RESULTS 

DIMENSION XCIOOD.YCIOO) 
C0MM0N/BLK7/NEIGV,KB,KBTG,XGRAPH 
XPG=5.0 
XDV=10. 

XTIC=.5 

YPG=5.0 

YDVrlO.O 

YTIC=.5 

CALL ASC ALEC X.XPG, 1,1, XDV} 

CALL ASCALECY.YPG.I.l.YDVD 

CALL AXESC0.,0.,0.,XPG,0.0,XCI + 2),-XTIC,XD V, 
1 20HV0LT AGE CMILLI VOLTS) ,.11,-20) 

CALL AXESCO.,0.,9 0., Y P G , 0 . 0 , Y C I + 2 ) , - YTIC, YD V 
1 19HCURRENT CMILLIAMPS) ,.11,19) 

SFX = l.^XCI + 2) 

SFYsl.^YCI + 2) 

CALL CALPLTCXC1)aSFX,YC1)aSFY,3) 

DO 100 J=1,I 

CALL CALPLTCXCJ)aSFX,YCJ)aSFY, 2) 
lOO CONTINUE 

CALL CALPLTCO.,6.,-3) 
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CALL CHARACTC0.0,0.0,.1.5HGRAPH,0.0,5) 
CALL WHERE(XX,YY,IXX} 

CALL NUMBER(XX + .07,0.,.l,XGRAPH,0.0 i -l) 

CALL NFRAME 

RETURN 

END 




TOTAL 



FIGURE Cl. Representative output from program "SC3D" illustrating the spectral 
response, the current-voltage relation and contour plot of current 
density in the planes Z=Z- and Z=Z-+W. 

J . J 
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APPENDIX D 

PROGRAM "VSC3D" FOR VERTICAL JUNCTION SOLAR CELL 
3-DIMENSIONAL 
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PROGRAM VSC3D(INPUT, OUTPUT, TAPE5=INPUT,TAPE6=OUTPUT,TAPE8D 
C 

C SC3D=SOLAR CELL 3 DIMENSIONAL (ILLUMINATED FROM SIDE X=A) 

C 

C GREEN’S FUNCTION SOLUTION 
C 

COMMON/BLKl/NDATA,ENERGC25D,ALP(25D 

COMMON/BLK2/POW,LAMBDA,NA,ND,B4,SN,SP,SO,Sl,RE,RS,RSH,A,B,H 

COMMON/BLK3/EPSO, AO, Al, A2, A3,BO,Bl,B2,B3,EGGO,ALPHAO,BETAO 

common/blk-^^xkko.mun.mup.taun.taup.dn.dp.phio.ni.ln.lp.vbi 

C0MM0N/BLK5/ETANC50D,XINC50),XIP(50),ETAP(50),XNC50),XPC50) 

COMMON/BLK6 / PI,PI2,Q,Z J,T,TO,NIAl,NIA2,TNO,TPO,W,ALPHA 

COMMON/BLK7/NEIGV,KB,KBTG,XGRAPH 

COMMON/BLK8/VOC,PMAX,EFF,ISC 

COMMON/BLK9/SPTOP.SNBOT 

DIMENSION ENERGS(25D,ENERGGA(25D,ENERGGE(25D 
DIMENSION ALPS(25),ALPGA(25),ALPGE(25D 


REAL LAMBD A,N A.ND.NI A1,NIA2,MUN,MUP,KB,KBTG,LN,LP,NI,ISC 
REAL NIASl,NIAS2,NIAGAl,6C>ICS.TTC’IAGE25*P8cK -, + I 


CALL PSEUDO 
NAMELIST VARIABLES 

NAMELIST/PARAM/P0W,NA,ND,LAMBDA,A,B,H,ZJ,T,RS,RSH,SN,SP, 
1 SO,Sl,NEIGV,ITYPE,ISWITCH 
C 

C POW = LASER POWER (KW/CM2D 
C N A=DOPING DENSITY (CM-3D 
C ND = DOPING DENSITY (CM-3D 
C LAMBDA=WAVELENGTH CMICRONS) 

C A = HALF WIDTH OF X-DIRECTION C -A .LE. X .LE. A) 

C B = HALF WIDTH OF Y-DIRECTION C -B .LE. Y .LE. B) 

C H = DEPTH OF CELL IN Z-DIRECTION C 0 .LE. Z .LE. HD 
C ZJ = JUNCTION DEPTH FROM Z = 0 CCM) 

C T = TEMPERATURE (DEG KD 
C RS = SERIES RESISTANCE OHMS 
C RSH--SHUNT RESISTANCE OHMS 

C SN = SURF ACE RECOMBINATION VELOCITY TOP (CM^SD 
C SP = SURFACE RECOMBINATION VELOCITY BOTTOM (CM/S) 

C SO=SURFACE RECOMBINATION VELOCITY TOP SIDES 
C Sl = SURF ACE RECOMBINATION VELOCITY BOTTOM SIDES 
C NEIGV = NUMBER OF EIGENVALUES 
C ITYPE = TYPE OF MATERIAL 
C ITYPE = 1 FOR SILICON (SID 

C ITYPE = 2 FOR GALLIUM ARSENIDE (GA'ASD 

C ITYPE = 3 FOR GERMANIUM (GED 

C ISWITCH TURNS CONTOUR PLOT ON AND OFF 
C ISWITCH=0 CONTOUR PLOT IS OFF 

C ISWITCH=1 CONTOUR PLOT IS ON 
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c 
c 

C NOMINAL VALUES FOR NAMELIST VARIABLES 

C 

POW= . 001 
LAMBDA: . 8 
TO: 300. 

RS : O . O 
SN: 1 . O 
SO: 1 . O 
NA=1. 25E17 
A: . 5 

H: 5 . OOE-4 
NEIGV: 15 
I SWITCH: 1 

C 

C SILICON PARAMETERS 

NDATAS=20 

DATA CENERGSCK),K:1,8)/1. 1,1. 15,1. 24, 1.378,1. 46, 1.77, 2. 067,2. 48/ 

DATA CENERGSCK),K = 9,15)/2.76,3.1,3.26,3.5,3.6,4.0,4.2/ 

DATA CENERGSCKD.K: 16,20)/ 5. 0,5.2, 6.0,8.0,10.0/ 

DATA CA LPSCK),K:1,7)/1. 0,10. 0,100. 0,350. 0,950.0, 2. E3, 4. 5E3/ 

DATA CALPSCK),K:8,14)/1.0E4,2.0E4,5.0E4,1.0E5,1.0E6,1.1E6,1.8E6/ 
DATA ( ALPSCK),K: 15,20)/ 2.5 E6, 1.9 E6.2.0E6.1.5E6, 1.3 E6.1.1E6/ 

C CONSTANTS FOR MOBILITIES 

DATA ASO,ASl,AS2,AS3/65.02,5.72E-9,2.4 9E9.2.5/ 

DATA BSO,BSl,BS2,BS3/54.46,6.76E-9,2.93E9,2.7/ 

C CONSTANTS FOR BANDGAP AS FUNCTION OF TEMPERATURE 
DATA EGGOS, ALPHAOS, BET AOS/1. 16.7.02E-4.1108./ 

C CONSTANTS FOR INTRINSIC CARRIER DENSITY 

DATA NIAS1.NIAS2/4.9E15.1.61/ 

C CONSTANTS FOR LIFETIMES 

DATA TNOS,TPOS,B4/1.86E-4,3.52E-3,7.1E15/ 

EPS0S:11.6 

C GALLIUM ARSENIDE PARAMETERS 

C ABSORPTION COEFFICIENT CURVE FOR GAAS 

NDATAGA:16 

DATA CENERGGACK),K:1, 8)/ 1.37, 1.38, 1.4, 1.41, 1.42, 1.43, 1.44, 1.5/ 

DATA (ENERGGA(K),K = 9,16)/1.6,1.8,2.0,2.4,3.0,3.5,4.0,4.5/ 

DATA C ALPGA(K),K:1,8)/18.,32.,119.,175.,2.5E3,4.5E3,6.E3,1.4E4/ 
DATA CALPGA(K),K:9,16)/1.8E4,3.E4,4.E4,1.E5,2.E5,3.E5,4.E5,5.E5/ 
C 

C CONSTANTS FOR MOBILITIES 

DATA AGA0,AGA1,AGA2,AGA3/5076.11,5.24E-9,6.43E5,1.0/ 

DATA BGA0,BGA1,BGA2,BGA3/1477.3,8.59E-8,6.05E7,2.1/ 

C CONSTANTS FOR BANDGAP 

DATA EGGOG A, ALPHG A.BETG A/1.519,5.4E-4,204./ 

C CONSTANTS FOR INTRINSIC CARRIER DENSITY 
DATA NI AG Al, NI AG A 2/4.9 E15,. 0488844/ 

C CONSTANTS FOR LIFETIMES 

DATA TN0GA,TP0GA,B4/8.8179E-7,1.404E-7,7.1E15/ 

EPSOG A:13 . 1 

C GERMANIUM PARAMETERS 
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$ RE: . 05 
$ I TYPE : 1 
$ T:300. 

$ RSH= 1 . 0E6 
$ SP:1.0 

$ 31=1.0 

$ ND= 5 . 0E1 9 
$ B= . 500 
$ ZJ= 2 . 5E-4 



C GERMANIUM SOLAR CELL PARAMETERS 
NDATAGE=22 

DATA CENERGGECK),K = 1,8)/.66,.68,.7,.73,.76,.79,.8,.81/ 

DATA CENERGGECK),K = 9,16)/.83,.85,.875,.9,.97,1.7,2.6,3./ 

DATA CENERGGECK),K = 17,22)/3.6,4.,4.4,5.0,6.0,8./ 

DATA CALPGECK),K = 1,8)/1.,10.,20.,60.,100.,600.,1.E3,4.E3/ 

DATA CALPGECK),K = 9,16)/5.E3,8.E3,8.5E3,9.E3,1.E4,1.£5.7.E5,8.E5/ 

DATA CALPGECK),K = 17,22)/ l.E6,1.4E6,2.0E6, 1.3 E6.1.2E6.1.0E6/ 

C 

DATA AGEO,AGEl,AGE2,AGE3/458.68,7.92E-9,6.08E7,1.6 6/ 

DATA BGE0,BGE1,BGE2,BGF.?/?140.27,1.79E-7,1.18E9,2.33/ 

DATA EGGOGE, ALPHGE,BETGE/.741,4.77E-4,235./ 

DATA NI AGEl.NI AGE 2/4.9 E15,. 33313/ 

DATA TN0GE,TP0GE,B4/2.1022E-5,2.1007E-5 I 7.1E15/ 

EPSOGE = 16. 

C 

C READ IN ANY CHANGES TO THE ABOVE PARAMETERS 
C 

XGRAPH=0. 

100 CONTINUE 

RE ADCS, P ARAM) 

IF(E0F(5))600,601 

600 WRITEC6.603) 

603 F0RMATC1X.28HEND OF FILE ENCOUNTERED-STOP ) 

STOP 1313 

601 CONTINUE 
C 

WRITEC6.HO) 

HO FORMAT(//) 

Q=1.6E-19 
XKKO = 1.2E-3 
PI = 3. 1415926536 
PI2=PI*.5 
KB = 8.6172E-5 
KBTG=KB*T 

C EPSO IS THE PERMITIVITY 

XGRAPH=XGRAPH+1 
C 

C SETUP PARAMETERS FOR CORRECT TYPE OF PHOTOCONVERTER MATERIALS 

C ITYPE = 1 CSI) 

C ITYPE = 2 CGAAS) 

C ITYPE = 3 CGE) 

IFCITYPE ,EQ. 1) WRITEC6.H1) 

IFCITYPE .EQ. 2) WRITEC6.112) 

IFCITYPE .EQ. 3) WRITEC6.113) 

111 FORMATC1X.25HSILICON PHOTOCONVERTER , / ) 

112 FORMATC1X.35HGALLIUM ARSENIDE PHOTOCONVERTER ,/) 

113 F0RMATC1X.25HGERMANIUM PHOTOCONVERTER , / ) 

IFCITYPE .EQ. 1) GO TO 10 
IFCITYPE .EQ. 2) GO TO 20 
C GERMANIUM DATA GOES HERE 

NDATArNDATAGE 
DO 40 I = 1,NDATA 
ENERGCI) = ENERGGECD 
40 ALP(I): ALPGECI) 
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A0= AGEO $ Al= AGE1 $ A2=AGE2 $ A3=AGE3 

BO = BGEO $ B1=BGE1 $ B2 = BGE2 $ B3=BGE3 

EGGO = EGGOGE $ ALPHAO= ALPHGE $ BETAO=BETGE 

NIA1 = NIAGE1 $ NIA2 = NIAGE2 

TNO=TNOGE $ TPO=TPOGE 

EPSO=EPSOGE 

GO TO 30 

10 CONTINUE 

C SILICON PARAMETERS 

NDATA=NDATAS 
DO 'll Isl.NDATA 
ENERGCD=ENERGSCD 

11 ALPCI) = ALPS(I) 

A0 = ASO $ Al = AS1 $ A 2 = A S 2 $ A3=AS3 

BO = BSO $ B 1 = B S 1 $ B 2 = B S 2 $ B3 = BS3 

EGGO=EGGOS $ ALPHAO=ALPHAOS $ BET A0= BET AOS 

NI A1=NI AS1 $ NIA2=NIAS2 

TNO=TNOS $ TPO=TPOS 

EPSO=EPSOS 

GO TO 30 

20 CONTINUE 

C GALLIUM ARSENIDE PARAMETERS 

NDATA=NDATAGA 
DO 12 I=1,NDATA 
ENERG(I) = ENERGGA(I) 

12 ALPCI) = ALPGA(I) 

A0 = AGAO $ Al = AG A1 $ A2 = AG A 2 $ A 3 = A G A 3 

BOrBGAO $ B1 = BGA1 $ B2 = BG A2 $ B3 = BGA3 

EGGO = EGGOGA $ ALPHAO=ALPHGA $ BETAOsBETGA 
NIAlsNIAGAl $ NI A2 = NI AGA2 
TNOrTNOGA $ TPO=TPOGA 
EPSOrEPSOGA 

30 CONTINUE 

C CALCULATE CONSTANTS 

CALL CONST 

C CALCULATE SPECTRAL RESPONSE 

CALL SPEC 

C -CALCULATE I-V CURVE 

CALL IV(LAMBDA) 

WRITE(6,199)XGRAPH 

199 FORMAT (T20.6HGRAPH ,F3.0 ) 

WRITEC6,200)ISC,V0C,EFF,RSH,DN,NEIGV 

200 FORMATC / ,T2,1HISC = ,E12.5,T20,1HV0C = ,E12.5,T10,1HEFF = , 

1 E12.5,T60,1HRSH = ,E12.5,T80 ) 3HDN = ,E12.5,T99,8HNEIGV = ,15) 

WRITE(6,201)NA,ND,ZJ,S1,DP,NIA1 

201 F0RMAT(T2,3HNA = ,E12.5,T2 0,3HND = ,E12.5,T10,3HZ J = ,E12.5, 

1T60,3HS1 = ,E12.5,T80,3HDP=,E12.5,T99,7HNIA1 = .E12.5) 

WRITEC6,202)TAUN,TAUP,LN,LP,TN0,NIA2 

202 F0RMAT(T2,5HTAUN=,E12.5,T20,5HTAUP = ,E12.5,T10,3HLN = ,E12.5, 

1 T60,3HLP=,E12.5,T80,6HTN0 = ,E12.5,T99,7HNI A2 = .E12.5) 

W R I T E ( 6 , 2 0 3 ) A , B , H , S 0 , T P 0 

203 F0RMATCT2 f 2HA=,E12.5,T20,2HB = ,E12.5,T10,2HH = ,E12.5,T6 0,3HS0 = , 

1 E12.5.T80.6HTP0 = ,E12.5) 

WRITEC6,201)MUN,MUP,LAMBDA,P0W,SPT0P 

201 F0RMAT(T2,1HMUN = ,E12.5,T20,1HMUP = ,E12.5,T10,7HLAMBDA = ,F10.1, 
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1T60,4HP0W = ,E12.5 ,T80,6HSPT0P= ,E12.5 ) 

WRITEC6,205)T,RS,SN,SP,SNB0T 

205 F0RMAT(T2,2HT = ,E12.5,T20,3HRS = ,E12.5,T 4 0,3HSN = ,E12.5, 

1 T60,3HSP = ,E12.5,T80,6HSNB0T = ,E12.5 ) 

WRITE(6,206)W,PM AX.NI, VBI.ALPH A 

206 F0RMATCT2,2HW=,E12.5,T20,5HPMAX = ,E12.5,T4 0,3HNI = ,E12.5, 

1 T60,4HVBI = t E12.5,T80,6HALPHA = ,E12.5,///) 

IFCISWITCH .EQ. 1) CALL CONTOUR 

GO TO 100 

END 

SUBROUTINE CONTOUR 

DIMENSION AAA(40,40),BBB(40,40) 

C AAA IS 40X40 ARRAY OF CURRENT DENSITY VALUES AT X,Y AND Z=ZJ 
C BBB IS 40X40 ARRAY OF CURRENT DENSITY VALUES AT X,Y AND Z=ZJ+W 
C 

N=40 

C CALCULATE JP AND STORE IN AAA. FIND AMAX AND AMIN VALUES FOR ARRAY 
C 

CALL CDENABCAAA.BBB.AMAX.AMIN.BMAX.BMIN.N) 

C 

C CALCULATE JN AND STORE IN BBB. FIND BMAX AND BMIN VALURS FOR ARRAY 
C 

C AAA ARE VALUES OF XJP 

C BBB ARE VALUES OF XJN 

C 

C ' r 

C 

C PLOT CONTOUR GRAPH OF AAA 

CALL CONRECC A A A,N,N,N,AMIN, AM AX,0,0,-1,0) 

CALL EZTRN 
CALL NFRAME 
WRITE(6,199) 

199 FORM AT(////) 

C PLOT CONTOUR GRAPH OF BBB 

C 

CALL CONREC(BBB, N,N,N, BMIN, BMAX, 0,0, -1,0) 

CALL EZTRN 
CALL NFRAME 
RETURN 
END 

SUBROUTINE CDENABCAAA,BBB,AMAX,AMIN,BMAX,BMIN,N) 
COMMON/BLK2/POW,LAMBDA,NA,ND,B4,SN,SP,SO,Sl.RE,RS,RSH,A,B,H 
C0MM0N / BLK4/XKK0,MUN,MUP,T AUN,TAUP,DN,DP,PHIO,NI,LN,LP, VBI 
COMMON/BLK6/PI,PI2,Q,ZJ,T,TO,NIA1,NIA2,TNO,TPO,W, ALPHA 
REAL NA.ND, LAMBDA, LN,LP,MUN,MUP,NI,NIA1,NIA2 
DIMENSION AAA(40,40),BBB(40,40) 

C CALCULATE JPCX.Y) AND STORE IN AAA 

C CALCULATE JN(X,Y) AND STORE IN BBB 

DELX = 2.* A/FLO AT(N-l) 

AMIN=5.E5 

AMAXrO.O 

DELYr2.*B/FL0AT(N+l) 

YWQ=Z J+W 
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WRITE(6,99) 

99 F0RMAT(T5.1HX,T1'5,1HY.T25.2HZ J.T43.2H JP.T55.1HX, 

1 T 65.1H Y,T75,-‘1HZ J + W,T93,2H JN ) 

BMIH=5.0E5 
BMAX=0.0 
DO 10 J=1,N 
DO 20 1=1, N 

X=-A+IaDELX 
Y=-B+JaDELY 

CALL SUMIJ(X,Y,XJP,XJN) 
AAACI,J)=2.0aQaALPHAaPHI0aC1-PF)aXJP 
BBB(I, J)=2.0aQa ALPHA aPHIOaC1-RE>XJN 
WRITE(6,100)X,Y,ZJ,AAACI,J),X,Y, YWQ.BBBCI.J) 

100 F0RMATC1X,3F10.4,E15.7,3X,3F10.6,E15.7) 

IF(AAACI,J).GT.AMAX) AM AX = A A A(I. J) 

IF( A A A(I, J).LT.AMIN) AMIN = AA ACI, J) 

IF(BBB(I,J).GT.BMAX) BMAX = BBB(I, J) 

IFCBBB(I.JD.LT.BMIN) BMIN = BBBCI,J) 

20 CONTINUE 
10 CONTINUE 
RETURN 
END 

SUBROUTINE SUMIJ(X,Y.X JP.XJN) 

C ' 

COMMON/BLK2/POW.LAMBDA,NA,ND,B4,SN,SP,SO,Sl,RE,RS.RSH,A,B,H 
COMMON/BLK4/XKKO,MUN,MUP,TAUN,TAUP.DN,DP,PHIO,NI,LN,LP,VBI 
COMMON / BLK6 / PI,PI2,Q,ZJ > T,TO,NIA1 I NIA2,TNO,TPO,W > ALPHA 
REAL N A, ND,MUN,MUP,NI,LN,LP,NIA1,NIA2.KB,KBTG, LAMBDA 
C 

XJP=0.0 

XJN=0.0 

DO 10 I=1,NEIGV 
DO 20 Jrl.NEIGV 

CALL TERNIJCX.Y.I.J.FPIJ.FNIJ) 

X JP=X JP+FPI J 
XJN=X JN+FNI J 
20 CONTINUE 
10 CONTINUE 
RETURN 
END 

SUBROUTINE TERNI J(X,Y,I, J.FPI J.FNI J) 
COMMON/BLK2/POW,LAMBDA,NA,ND,B4,SN,SP,SO,S1,RE,RS,RSH,A,B,H 
C0MM0N/BLK4/XKK0,MUN,MUP,TAUN,TAUP,DN,DP,PHI0,NI,LN,LP,VBI 
C0MM0N/BLK5/ETAN(50),XIN(50),XIF(50),ET AP(50),XN(50),XP(50) 
C0MM0N/BLK6/PI,PI2,Q,ZJ,T,T0,NIA1,NIA2,TN0,TP0,W,ALPHA 
COMMON/ BLK7/NEIGV, KB, KBTG.XGRAPH 
REAL NA,ND,MUN,MUP,NI,LN,LP,NIA1,NI A2.KB.KBTG 
REAL LAMBDA 
EXTERNAL FUN 
C 

GCZ) = FUN(-ALPHAaZ) 

C FPIJ (TOP REGION) N-REGION 

WJPNOR = .5aZJ-.25aSINC2aXP(J))aZJ/XPCJ) 
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XMUI = ETAP(IKB 

XLAMJ = CXP(J)KZJ 

DEN2 = (XMUIaa2)+(XLAMJaa2) 

GIJ2 = LPaLP/(1. + LPaLPaDEN2) 

GIJ=SQHTCGIJ2D 

CC=ALPHAaALPHAaGIJ2-1 

VIPN0R = B + CS0aBaB/DP)aCC0S(ETAPCI))/ETAPCI))aa2 
SPI1=BaSIN(ETAP(I))/ET apci) 

SPI2 = ZJa( 1.-COSCXPC J)))/XPC j) 

SPI3=XLAMJ/WJPNOR 
SPI^ = COS(XHUIaY)/VIPHOR 
CALL SUB1(X,A,G1J,S0,DP,ALPHA,CC,SPI5) 
FPIJ=SPI1aSPI2aSPI3aSPI4aSPI5 

C CALCULATE FNIJ BOTTOM REGION (P-REGION) EVALUATE AT Z=ZJ+W 

HP = H-(Z J+W) 

WJNNOR = .5aHP-.25aSINC2aXNCJ))aHP/XNCJD 
VINNOR = B + CS1aBaB/DN)aCCOSCETANCI)XETAN(I))aa 2 
XMUI = ETAN(IKB 
XLAMJrXNCJKHP 

DEN3=CXLAMJaa2)+CXMUIaa2) 1 

SIJ2 = LNaLN/C1.+LNaLNaDEN3) 

SIJ=SQRTCSIJ2) 

EE^ALPHAaALPHAaSI J2-1 
SNIlrBASINCETANCDKETANCI) 

SNI2 = HPaC1.-COSCXN(J))VXNC J) 

SNI3r XL AM J/ W J N N O R 
SNI4 =COS(XMUIaY;KVINNOR 
CALL SUBKX.A.SIJ.Sl.DN, ALPHA, EE, SNI5D 
FNIJ = SNI1 aSNI2aSNI3aSNI4aSNI5 
RETURN 
END 

SUBROUTINE SUB1(X, A,GI J.SO.DP, ALPHA,CC,SPI5) 

EXTERNAL FUN 
G2=GIJaGI J 
Gl=SO/DP 
AB1=G1- ALPHA 
AB2=G1+ ALPHA 

AO=G2aAB1aEXP(-2.aAaALPHA) / CC 

Al = AOaGI JaGI 

A2=G2aAB2/CC 

A3=A2aGIJaG1 

BO=2.aG1 

B1 = (1. / GIJ) + G1aG1aG2 
H=ALPHAa(X-A) 

Q1 = TANHC2.aA/GIJ) 

R1 = -G2aBOaEXPCH)/CC 

R2 = C-GIJa(1 + G1aG1aG2>EXPCH)/CCDaQ1 

ARG1 = -CX + AKGIJ 

ARG2=(X-3aA)/GIJ 

ARG3 = CX-A)/GIJ 

ARG4 = -(X + 3aA;KGIJ 

DEN=1+FUN(ARG4) 

R3=AOa(FUNCARG1)+FUNCARG2)KDEN 
R^ = A1a(FUN(ARG1)-FUNCARG2))/DEN 
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R5 = A2*CFUN(ARG3>FUN(ARG4)KDEN 
R6=A3*(FUN(ARG3)-FUN(ARG-'!))' / 'DEN 
SPI5 = (R1 + R2 + R3 + Ri + R5 + R6) / (B0 + B1*Q1) 

RETURN 

END 

SUBROUTINE CONST 

CALCULATES OTHER CONSTANTS AND PARAMETERS FOR PROBLEM 


COMMON/BLK2/POW,LAMBDA,NA,ND,B4.SN,SP,SO,Sl,RE,RS,RSH,A,B,H 
C0MM0N / BLK3 / EPS0,A0,A1,A2,A3,B0,B1,B2,B3,EGG0 1 ALPHA0,BETA0 
COMMON/BLK^/XKKO,MUN,MUP,TAUN,TAUP,DN,DP,PHIO,NI,LN.LP,VBI 
COMMON / BLK6 / PI,PI2,Q,ZJ,T,TO,NIA1,NIA2,TNO,TPO,W, ALPHA 
COMMON/BLK7/NEIGV,KB,KBTG,XGRAPH 


REAL MUN,MUP,KBTG,KB,L AMBD A,N A,ND,NIA1,NIA2,NI,LN,LP 


CONVERT WAVELENGTH CMICRONS) TO ENERGY (EV) 
ENER=1.2402/LAMBDA 
C CALCULATE BANDGAP AS FUNCTION OF TEMPERATURE 
EGTrEGGO-ALPHAO*T*T/(T+BETAO) 

C CALCULATE ENERGY SHIFT DUE TO TEMPERATURE 
DELE=XKKOx(T-TO) 

E=ENER+DELE 

C CALCUA-LTE ABSORPTION COEFFICIENT 
CALL ABSORP(E,ALPHAl) 

ALPHA=ALPHA1 

C CALCULATE LIFETIMES 
TAUN = TNO/(l.+NA/B4) 

TAUP = TPO/(l.+ND/BiD 
C CALCULATE MOBILITIES 
CALL MOBILCY1.Y2) 

MUN= Y 1 $ MUP = Y2 

C CALCULATE DIFFUSION COEFFICIENTS 
DN=KBTGxMUN $ DP=KBTG*MUP 
C CALCUALTE DENSITY OF PHOTONS IN PHOTON’S/CM2-S 
PHI0=P0W*LAMBDA*5.03306E21 
C CALCULATE INTRINSIC CARRIER DENSITY 
POWER = -EGT/(2.*KBTG) 

CALL EXPON(POWER.ANS) 

NI = NIA1*NIA2*(T**1.5)*ANS 
C CALCULATE BUILT IN VOLTAGE 

VBI = KBTG* ALOG(NA*ND'"(;NI*NI)) 

V0LTS=0.0 

CALL WIDTHC VOLTS) 

C CALCULATE DIFFUSION LENGTHS 

LP = SQRTCDPxTAUP) $ LN=SQRT(DN*TAUN) 

C CALCULATE ALL NECESSARY EIGENVALUES FOR PROBLEM 
CALL ESET 
RETURN 
END 
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SUBROUTINE ABSORP(E, ALPHA) 

C 

C 

COMMON/BLKl/NDATA,ENERGC25),ALPC25) 

C 

C 

C 

1=0 

IFCE.LT.ENERGCUDGO TO 10 
20 CONTINUE 
1=1+1 

IF(I.GE.NDATA)GO TO 30 

IFCCENERGCI).LE.E).AND.CENERG(I + 1) .GE. E))GO TO 40 
GO TO 20 

"10 ALPHA = ALPCI) + (ALPCI + 1)-ALPCI))*(E-ENERGCI))''CENERG(I + 1)-ENERGCI)) 

RETURN 

10 A L P H A = 0 . 

RETURN 

30 ALPHA = ALP(NDATA) 

RETURN 

END 

SUBROUTINE EXPON(X.ANS) 

IFCX.LT.-500) GO TO 10 
IF(X .GT. 700)G0 TO 20 
ANS = EXPCX) 

RETURN 
10 A N S = O . O 
RETURN 

20 WRITEC6.30) 

30 FORMATC1X.29H ARG TOO LARGE, EXP FUNCTION ) 

STOP 444 
END 

SUBROUTINE ESET 
C 
C 

COMMON/BLK2/POW,LAMBDA,NA,ND,B4,SN,SP,SO,Sl,RE,RS,RSH,A,B,H 
C0MM0N / BLK4/XKK0,MUN,MUP,T AUN,TAUP,DN,DP,PHIO,NI,LN,LP, VBI 
C0MM0N/BLK5/ETANC50),XINC50),XIPC50),ETAPC50),XNC50),XPC50) 
COMMON / BLK6 / PI,PI2,Q,2J,T 1 TO,NIAl,NIA2,TNO,TPO,W, ALPHA 
COMMON/BLK7/NEIGV,KB,KBTG,XGRAPH 
C 

REAL MUN, MUP,NI,LN,LP,NIA1,NIA2, KB, KBTG.NA.ND, LAMBDA 
C 

C CALCULATE EIGENVALUES FOR TOP REGION X-DIRECTION 
R1 = S0 $ R2 = DP $ R3 = A $ N4=NEIGV 
CALL EIGEN(XIP,R1,R2,R3,N4) 

C CALCULATE EIGENV VALUES FOR TOP REGION 2-DIRECTION 
R1 = SP $ R2 = DP $ R3=2J 
CALL EIGCXP,R1,R2,R3,N4) 

C CALCULATE EIGENVALUES FOR TOPphLGION Y-DIRECTION 
R1 = S0 $ R3=B 

CALL EIGEN(ET AP,R1,R2,R3,N4) 

C CALCULATE EIGENVALUES FOR BOTTOM REGIION 
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R1 = S1 $ R2 = DN $ R3= A 

CALL EIGEN(XIN.R1,R2,R3,N4} 

CALCULATE EIGENVALUES BOTTOM REGION Z-DIRECTION 
R1=SN $ R2 = DN $ R3 = H-ZJ-W 
CALL EIG(XN,R1,R2,R3,N4) 

C CALCULATE EIGENVALUES FOR BOTTOM REGION 
RlrSl $ R 3 = B 

CALL EIGENCETAN,R1,R2,R3,N4) 

WRITE(6,99XXIP( J),J = 1,N4) 

WRITE(6,98) 

WRITEC6,99)CETAPCJ),U = 1 : N4) 

WRITEC6.98) 

WRITEC6.99DCXINC J), J = 1,N4) 

WRITE(6,98) 

WRITE(6,99DCETANC J), J = 1,N4) 

WRITEC6.98) 

WRITEC6.99) (XP( J), J = 1.N4} 

WRITEC6.98) 

WRITE(6,99) (XNC J),J = 1,N4D 
WRITEC6.98) 

99 FORM AT(1X,C5E15.6)D 
98 FORMAT(//) 

RETURN 

END 

SUBROUTINE EIGCXX.S.D.H.N} 

DIMENSION XXC50) 

DOUBLE PRECISION X,X1,X2, Y,Y1,Y2,PI,PI2 
FCZ) = 1 + CaZaDCOSCZ)/DSIN(Z) 

C N MUST BE LESS THAN 50 
C=D/(SaHD 

PI = 3. 1415926535897932 
FI2=.5*PI 

C XX=ETAJ*Z JUNCTION 

IFCS.LE.O^GO TO 1001 
C GUESS AT FIRST ROOT 
X1=PI2 

M=1 $ IC=1 

200 Y1=F(X1) 

IFCYl.GT.O)GO TO 10 
C Y 1 IS NEGATIVE, DECREASE X 
X2=X1 $ Y2=Y1 

5 X1 = .5*(X2 + (2*M-1)*PI2) 

YlsFCXl) 

IF(Yl.GT.O) GO TO 20 
C Y 1 IS NEGATIVE 

X2 = X1 $ Y2 = Y 1 
GO TO 5 

10 CONTINUE 

C Y 1 IS POSITIVE, INCREASE X 
X2=.5*(X1 + M*PI) 

Y2=F(X2) 

IFC Y2.LT.0) GO TO 20 
X1 = X2 $ Yl = Y 2 
GO TO 10 
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OF. POOR QUALITY 


£0 X : ( X 1 + X 2 ) / 2 . 

V=FCX) 

IFCDABSCY) .LT. l.OE-6) GO TO 100 
IFCY.LT.0j GO TO 30 
X 1 = X $ Y 1 = Y 

GO TO 20 
3 0 X 2 : X 
Y2= Y 

GO TO 20 
100 CONTINUE 
■ ZZ=X 
XXCICj=2Z 

M = M + 1 $ I C - I C * 1 

XlrX+PI 

IFCIC.GT. N) GO TO 308 
GO TO 200 
308 RETURN 
1001 WRITE(6,234) 

£3-4 FORM A TCI X.-45HRE COMBINATION VELOCITY OF ZERO NOT ALLOWED ) 
END 

SUBROUTINE EIGENCXX, S.D.H.N) 

C 

COMMON/BLK6-''PI,PI2.Q,ZJ,T.TO,NIAl,NIA2.TNO,TPO,W, ALPHA 
. C 

REAL NIA1.NIA2 
C 

DIMENSION XXCSOj 
FCZj = TAMCZj-C'''Z 
C SOLVES EQUATION X*TANX:C 
C N MUST BE .LE. 200 
C=SaH/D 

IFCC.EQ.O) GO TO 1001 
C GUESS AT FIRST ROOT 
X1 = PI''4. 

Mrl $ IC=1 
200 Y1=.FCX1) ' 

IFCYl.GT.OjGO TO 10 
C Y 1 IS NEGATIVE 
X2 = X1 
Y 2 = Y 1 

5 X 1 = . 5 * C X 2 * C £ a M - 1 ) .* P 1 2 ) 

Y1 = FCX1) 

IFCYl.GT.OjGO TO £0 
C Y 1 IS NEGATIVE 
X2 = X1 
Y2=Yi 
GO TO 5 

10 X 2 = C M - 1 j P I + . 8 .*■ C X 1 - C M - 1 j a p I j 
Y2=FCX2j 

IFCY2.LT.0jG0 TO 20 

X1=X2 

Y1=Y2 

GO TO 10 

2 O X - ( X 1 + X C j 7 2 . 
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30 


Y=F(X) 

IFCABSCY) .LT. 1.0E-6)G0 TO 100 
IF(Y .LT. O) GO TO 30 
X1=X 
Y1=Y 

GO TO 20 
X 2 = X 
Y2=Y 

GO TO 20 ' 

100 XX(IC) = X 

M = M + 1 $ IC = IC + 1 

X1=X+PI 

IFCIC.GT.N) GO TO 308 
GO TO 200 
308 CONTINUE 
RETURN 

1001 WRITEC6.234) 

23-4 FORMATC45HRECOMBINATION VELOCITY OF ZERO NOT ALLOWED ) 

END 

SUBROUTINE MOBILCXMUN.XMUP) 


COMMbN/BLK2/POW,LAMBDA,NA.ND,B^,SN,SP,SO,Sl,RE,RS.RSH.A,3,H 
COMMON/BLK3/EPSO, AO, Al, A2,A3,BO,B1,B2,B3,EGGO,ALPHAO,BETAO 
COMMON/BLK4/XKKO,MUN,MUP,TAUN,TAUP,DN,DP,PHIO,NI,LN,LP,VBI 
C0MM0N/BLK6/PI, PI2, Q, ZJ,T, TO, NIA1.NIA2,TN0,TP0,W, ALPHA 


REAL LAMBDA,NA,ND,MUN,MUP,NI,LN,LP,NIA1,NIA2 

MOBILITY OF ELECTRONS 
RQ=2./3. 

BB1=A1aTaT/CNAaaRQ) 

D = NAaBB1a(1.-.5aBB1) 

X1=A0aCTaa1.5)/D 
ZZ=- A3a ALOG(T) 

CALL EXPON(ZZ.ANS) 

X2=A2aANS 
R = l./X2 + 1./X1 

XMUN=1./R 
C MOBILITY OF HOLES 
ZZZ=-B3a ALOGCT} 

CALL EXPON(ZZZ.ANS) 

X1=B2aANS 

BB1=B1aTaT/(NDaaRQ3 
D = NDaBB1a(1.-.5aBB1) 

X2 = BOa(Taa1.5)/D 

R = l./X2 +1./X1 

XMUP=1./R 

RETURN 

END 
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SUBROUTINE WIDTH(VOLTS) 


COMMON /BLK2/POW.L AM BDA,NA,ND,B^,SN,SP,SO,Sl,RE,RS,RSH,A,B,H 
COMMO N/BLK3/EPSO, AO, Al, A2, AB.BO.Bl.BS.BS.EGGO.ALPHAO.BETAO 
COMMON/BLK^/XKKO,MUN,MUP,TAUN,TAUP,DN,DP,PHIO,NI.LN,LP,VBI 
C0MM0N / BLK6 / PI,PI2,Q,ZJ,T, TO, NIAl,NIA2,TNO,TPO,W, ALPHA 
C0MM0N/BLK7/NEIGV.KB,KBTG,XGRAPH 

REAL KB,KBTG,LAM3DA,LN,LP,MTTN,MUP,NA,ND,NI,NIA1,NIA2 
ABC=C1. / ND + l./NA) 

W2=(2.*EPSO*(8.85E-14)/Q)*(VBI- VOLTS)* ABC 
W = SQRT(W2) 

RETURN 
END 

SUBROUTINE SPEC 

CALCULATE SPECTRAL RESPONSE 

COMMON/BLK2/POW,LAMBDA,NA,ND,B4,SN,SP,SO.Sl,RE,RS,RSH,A,B,H 
COMMON/BLk3/EPSO,AO,Al,A2,A3,BO,Bl,B2,B3,EGGO,ALPHAO,BETAO 
COMMON/BLK^/XKKO.MUN.MUP.TAUN.TAUP.DN.DP.PHIO.NI.LN.LP.VBi 
COMMON/BLK6 :/ PI,PI2,Q,ZJ,T, TO, NIAl,NIA2,TNO,TPO,W, ALPHA 
COMMON/BLK9/SPTOP.SNBOT 
C 

REAL LAMBDA.NA,ND,MUN,MUP,NI,LN,LP,NIA1,NIA2 
C 

DIMENSION WAVC200),SSPC200),SSNC200),SRC200) 

C 

SPTOP=-l. 

SNBOT=-l. 

1=0 

9 0 1 = 1 + 1 

E = EGG0-.15 + .025*CI-1) 

IFCCE.GT.5.0).0R.(I.GT.200))G0 TO 100 
WAV(I) = 1.2402/E 

FFF = W AV (I)*Q*PO W *(5.0 3306 E21)*( 1-RE) 

PIN=2*B*H*P0W*1000. 

CALL SHORTCS1P.S1N.E) 

SSN(I)=S1N*FFF / PIN 
SSPCI) = S1P*FFF/PIN 
SR(I) = SSN(I) + SSP(I) 

GO TO 90 
100 CONTINUE 
MAX=I-1 
1 = 1-1 

C NORMALIZE SPECTRAL RESPONSE 

XMAX=0. 

DO 110 K=1,MAX 
IF(SR(K).GT.XMAX)XMAX = SR(K) 

110 CONTINUE 

DO 120 K=1,MAX 
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SSN(KD = SSNCKKXMAX 
SSPCK^=SSPCKKXMAX 
SRCK)rSRCK)/XMAX 

IFCSSNCKD .GT. SNBOT) SNBOT=SSNCK) 

IF(SSP(K} .GT. SPTOP) SPTOP=SSP(K) 

120 CONTINUE 
C PLOT GRAPH 

CALL GRAPHCWAV.SSN.SSP.SR.I) 

RETURN 

END 

SUBROUTINE GRAPHCWA V.SN.SP.SR.I) 

DIMENSION W A VC200),SN(200),SP(200),SRC200} 

CALL CALPLTC1.5,.75,-3) 

XPG=5.0 

XDVrlO.O 

XTIC=.5 

YPG=5.0 

YDV=10.0 

YTIC=.5 

CALL ASC ALE( WAV, XPG, 1,1, XDV) 

CALL ASC ALE(SR, YPG,I,1,YDV) 

CALL AXES(0.,0.,0.,XPG,0.0,WAV(I + 2),-XTIC,XDV, 
120HW A VELENGTH (MICRONS) ,.14,-20) 

CALL AXESC0.,0.,9 0.,YPG,0.0,.2,-YTIC,YDV, 

1 28HNORMALIZED SPECTRAL RESPONSE ,.14,28) 

SFX = l./WAV(I + 2) 

SFY=5.0 

CALL CALPLT(WAV(1)aSFX,SN(1)aSFY,3) 

DO lOO J=1,I 

CALL CALPLT(WAV(J)aSFX,SN(J)aSFY,2) 

100 CONTINUE 

CALL CALPLT(WAV(1)aSFX,SPC1)*SFY,3) 

DO 101 Jsl.I 

CALL C ALPLT( WA V( J)aSFX,SP( J)aSFY,2) 

101 CONTINUE 

CALL CALPLTC W A V(1)aSFX,SR(1)aSFY,3) 

DO 103 Jsl.I 

CALL CALPLT(WAVCJ)aSFX,SRCJ)*SFY,2) 

103 CONTINUE 

CALL CALPLTCO.,6.,-3) 

RETURN 

END 

SUBROUTINE SH0RT(S1F,S1N,E) 


C0MM0N/BLK2/P0 W,L AMBD A,N A,ND,B4,SN,SP,S0,S1,RE,RS,RSH,A,B,H 
COMMON/ BLK4/XKKO,MUN,MUP,T AUN.TAUP.DN.DP.PHIO.NI.LN.LP, VBI 
COMMON/BLK6/PI,PI2,Q,ZJ,T, TO, NI Al.NI A2.TN0.TP0, W, ALPHA 
COMMON/ BLK7/NEIGV, KB, KBTG.XGRAPH 
C 

REAL LAMBD A,NA,ND,MUN,MUP,NI,LN,LP,KB,KBTG,NIA1,NIA2 
C 
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C CALCULATES THE SPECTRAL RESPONSE COMPONENTS 
C CALCULATES DOUBLE SUMMATION OVER I.J 
CALL ABSORP(E.ALPl) 

ALPHA=ALP1 

S1P = 0 $ SIN = O 

DO 10 I=1,NEIGV 

DO 20 J=1,NEIGV 

CALL NUMBSCI.J.SNIJ.SPIJ) 

S1P=S1P+SPIJ 
S1N=S1N+SNI J 
20 CONTINUE 
10 CONTINUE 
RETURN ■ 

END 

SUBROUTINE DCURENT(IO} 


COMMON/BLK2/POW,LAMBDA,NA,ND,B^,SN,SP,SO,Sl,RE,RS,RSH,A,B,H 
COMMON/BLK4/XKKO,MUN,MUP,TAUN,TAUP,DN,DP,PHIO,NI,LN,LP,VBI 
COMMON/BLK5/ETANC50D,XINC50),XIPC50),ETAPC50),XNC50),XPC50) 
COMMON/BLK6/PI,PI2,Q.ZJ,T,TO,NIAl,NIA2,TNO,TPO,W, ALPHA 
COMMON/BLK7/NEIGV,KB,KBTG,XGRAPH 
C 

REAL NN0,NP0,LAMBDA,NA,ND,MUN,MUP,NI,LN,LP,NIA1,NIA2 
REAL KB.KBTG.IO.NN.NP 
C 

Y1CX,Z3=ZaC0SHCX>SINHCXD 
Y2CX,ZD = ZaSINHCX)+C0SHCXD 
Y3CX,Z) = CZ + TANHCX))/C1 + 2xTANHCX)) 

C DOPING DENSITIES 
PPO=NA 
NNO=ND 

NPO=NIaNI/PPO 
PNO=NIaNI/NNO 
SUM1 = 0 $ SUM 2 = 0 

DO lO I=1,NEIGV 
DO 20 Jsl.NEIGV 

DEN2 = (ETAP(J;KB)aa2 +CXIPCI) / A)aa2 
GIJ2 = LPaLP/(1.+LPaLPaDEN2) 

GIJ=SQRT(GIJ2) 

DEN3=CXIN(I)/A)aa2 +(ET AN(J^B)aa2 
SI J2 = LNaLN / (1.+LNaLNaDEN3) 

SIJ=SQRTCSI J2) 

NP = SPaGIJ^DP $ NN = SNaSIJ/DN 
UJNNORMsA + (SOaAaA''DN)a(COS(XIN(J}>'XINCJ))a*2 
VINN0RM=B+(S1aBaB/DN)aCC0SCETANCI)D''ETANCI)Daa2 
F1 = SIN(XIN(J)) $ F2 = SIN(ETAN(I)) 

hp=h-zj-w 

F3 = 16aAaAaBaBaQaDNaNIaNI/(NAaSIJ) 

F^ = CF1/XIN(J))aa2 $ F5 = CF2/ET AN(I))aa2 

F6=1./CUJNNORMaVINNORM) 

TERMN = F3aF^aF5aF6aY3CCHP/SIJ),NN) 

UJPN0RM = A-*-CS0aAaA/DP)a(C0S(XIPCJ)) / XIPCJ)Daa2 
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VlPNORM = B-f(Sl*B*B/DP)*(COS(ET AP(I)) / ETAF(I))**2 
F1 = SIN(XIP(J)) $ F2= SIHCETAPCD) 

G3 = 16aAaAaBaBaQaDPaNIaNI/(NDaGIJ;) 

G4 = (F1/XIP( J))aa2 $ G5 = (F2>'ETAPCI))**2 

G6=1./CUJPN0RMaVIPN0RM) 

TERMP = G3aG4aG5aG6aY3C(2J/GIJD,NP) 
SUM1=SUM1+TERMN • 

SUM2=SUM2+TERMP 
20 CONTINUE 
10 CONTINUE 

I0 = SUM1+SUM2 

RETURN 

END 

SUBROUTINE NUMBSCI.J.SNIJ.SPI J) 


COMMON''BLK2''POW,LAMBDA 1 NA.ND,B4,SN,SP,SO,Sl,RE,RS,RSH,A,B,H 
COMMON/BLK^/XKKO.MUN.MUP.TAUN.TAUP.DN.DP.PHIO.NI.LN.LP.VBI 
C0MM0N/BLK5/ETANC50),XINC50),XIPC50).ETAPC50),XNC50),XPC50D 
COMMON/BLK6/PI, PI2, Q, ZJ.T, TO, NIA1.FIQ2.TN0.TP0, W, ALPHA 
C0MM0N/BLK7/NEIGV.KB.KBTG.XGRAPH 

REAL L AMBD A,N A,ND,MUN,MUP,NI,LN,LP,NIA1,NIAB,KB,KBTG 

AAAAA a AAA A AAAAAAA AAAAAA AAA A A AAA A AAAAAAAA AAA A A A A A A A A A A A A A 


WJPNOR=.5aZJ-.25a2JaSINC2aXP(J)KXPCJ) 
DEN2=CETAPCn/B)AA2 + (XP( JXZ J>a2 
GIJ2 = LPaLP/(1.>LPaLPaDEN2) 

GIJ=SQRT(GIJ2) 

VIPN0R = B + CS0aBaB/DP)a(C0S(ETAPCI))/ETAPCID)aa2 
HP=H-ZJ-W 

WJNNOR = .5aHP-.25aHPaSIN(2aXNCJ))/XNCJ) 

VINNORrB + CSlABAB/DNXCOSCETANCDJ/ETANCI^AAB 
DEN3 = CXN(J)/HP)aa2 +(ETANCIKB)aa2 
SIJ2 = LNaLN/(1.>LNaLNaDEN3) 

SIJ=SQRT(SIJ2) 

XLAM J=XP(J)/ZJ 
SPI1=^aBaBaXLAMJ/WJPN0R 
SPI2 = ZJa(1-C0S(XP(J))) / XPCJ') 

SPI3 = CSIN(ETAPCIDKETAP(ID)aa2 

CC=ALPHAaALPHAaGIJ2-1 

AA=A 

AALP=ALPHA 

DPP=DP 

SSO=SO 

DD=GIJ/CC 

CALL INTEGC A A, A ALP.GI J,CC,DP P.SSO.SPH) 

SPIJ = ALPHA aSF II aSPI2aSPI3aSPI4-'''CVIPNOR) 

C AAAAAAAAAAAAAAAAAAAAAAAAAAAA AAAAA AAAAAAAAAAAAAAAAAAAAAAAA 

SNI1 = ^aBaB 

XLAMJ =XN( J)/HP 

SNI2 = HPa(1-C0S(XN(J)))/XNCJ) 


C ~3U 
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SNI3=(SIN(ETAN(I))/ETAN(I))aa2 
EE=ALPHAa ALPHA aSIJ2-1 ’• 

SS0=S1 

DPP=DN 

SNH=XLAMJ / W JNNOR 

CALL INTEGC A A, A ALP.SIJ.EE, DPP, SSO.SNI5) 

SNIJ=ALPHAaSNI1aSNI2aSNI3aSNI4aSNI5/(VINNOR) 

RETURN 

END 

SUBROUTINE INTEGC A A. A ALP.GI J.CC.DP.SO.SPH) 

ARG1=2aAA/GIJ 

ARG2=2aAAa AALP 

CALL EXPON(-ARG2,ARG3) 

ARGi:l-ARG3 $ ARG5 = 1 + ARG3 
ARG6 =GIJaa4 
DP2=DPaDP 
GIJ2=GIJaGI J 
S02=S0aS0 

IFCAALP.LE.O)GO TO 200 
A11 = ~2aSOaGIJ2/CAALPaDP) 
A12=AALPaSOaARG6/DP 
A13 = S02aARG6/(DP2) 

A1 = (C A11 + A12)a ARG-4 + A13aARG5)/CC 
C CC=ALPHAaALPHAaGI J2-1 

B11 = -GIJaC1+S02aGIJ2/CDP2))/AALP 
B12=A ALPa(GI Jaa3) 

B13 = SOa(GIJaa3D/DP 

B1 = (CB11 + B12)aARG^ + B13aARG5)/CC 

C11=-SOaAALPaARG6/DP 

C12 = -S02aARG6/(DP2D 

C1 = (C11aARG4 + C12aARG5;KCC , 

50 B = (1./GI J) + (S02aGI J/CDP2)) 

A=2aS0^DP 

R=TANH(ARG1) 

IF( ARG1.GT.700) GO TO 100 

SPM=CA1 + B1aR + C1/C0SHCARG1))/(A+BaR) 

RETURN 

100 SPM = (A1 + B1aR)/CA>BaR) 

RETURN 

200 CONTINUE 

A11 = -^aS0aGIJ2aAA/DP 
A12=2aS02aARG6 / DP2 
A1 = CA11 + A12) / CC 

B11 = -GIJ2aC1 + S02aGIJ2/DP2)a2aAA 
B12=SOaGIJaGIJ2a2/DP 
B1=(B11+B12)/CC 
C1 = -S02aARG6a2/(DP2aCCD 
GO TO 50 
END 
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FUNCTION FUNCX) 

IF(X.LT.- 2 30) GO TO lO 
FUN=EXPCX) 

RETURN 
10 FUN:0 

RETURN 
END 

SUBROUTINE IV(WAVE) 

CALCULATE CURRENT-VOLTAGE CURVE 


COMMON/BLK2/POW,LAMBDA,NA,ND.B4,SN,SP,SO,Sl,RE,RS,RSH,A,B,H 

COMMON / BLK^ / XKK0 1 MUN l MUP,TAUN t TAUP,DN,DP,PHIO,NI,LN,LP,VBI 

COMMON/BLK5/ETAN(50),XINC50),XIPC50),ETAP(50),XNC50),XP(50) 

COMMON/BLK6/PI,PI2,Q,2J,T,TO,NIA1 1 NIA2,TNO,TPO,W, ALPHA 

COMMON/BLK7/NEIGV,KB,KBTG,XGRAPH 

COMMON/BLK8/VOC,PMAX,EFF,ISC 


DIMENSION AMPS(IOO), VOLTS(IOO) 

REAL IREC,IB,IO,ISC,I 

REAL LAMBDA.NA.ND.LN.LP.MUP.MUN.NIAl.NIAa.KB.KBTG.NI 
C CALCULATE IV CURVE 
FFF = Q aPHIOa(I-RE) 

E=1.2^02^WAVE 
ICOUNTrO $ PM AX=-.5 

DEL V0 = . 0125 
DELV=.025 
V=-DELV 
lOO CONTINUE 

IC0UNT=IC0UNT+1 
V=V+DELV 
CALL WIDTH(V) 

C CALCULATE ISC (LIGHT CURRENT) 

CALL SHORTCS1P.S1N.E) 

YAl=FFFxSlP 
YA2 =FFFaS1N 
ISC= YA1+YA2 

C CALCULATE DARK CURRENT 
CALL DCURENT(IO) 

IREC = 4aAaBaQ*NIaWaPI2a2aSINHCV/(2aKBTG)) 
IREC=IREC/((VBI-V)/KBTG) 

IREC=IREC/SQRT(TAUNaTAUP) 

C SOLVE FOR I 

U1=1+RS/RSH 

U=.999 

IB = ISC-IREC-V/RSH+IO 
V0LTS(IC0UNT)=Va1000 
110 CONTINUE 

Z=(V + UaIBaRS)/KBTG 
F = IBaUaU1-IB-I0aEXP(Z) 

FP = IBaU1 + I0a(RSaIB/KBTG)aEXP(2) 

U2=U-F/FP 
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120 


ERR0R=ABSCU2-UD 
IFCERROR.LT.l.OE-6) GO TO 120 
U=U2 

GO TO 110 
I = U 2 a I B 

AMPSCIC0UNT)=Ia 1000 
WRITEC6,99)ISC,IO,IREC,V,YAl,YA2 

99 FORMATC1X.6E15.7D 

P=IaV 

IFCP.GT.PM AX)PMAX=P 

IF(CI.LT.O).OR.f VBI.LT.C V+DEL V>3)GO TO 130 
IF( V.GT..85} DELV = DELVO 
GO TO 100 

130 AMPSCIC0UNT) = 0. • 

PIN=P0Wa1000a2aBaH 

V0C = KBTGaAL0G(1.+CISC-IRECVI0) 

EFF=PMAXa100/PIN 

IFCICOUNT .GT. 1) V0LTSCIC0UNTD=V0LTSCIC0UNT-1)+DELV/100. 
IF C1C0UNT .EQ. 1) VOLTSCICOUNTD=VOC 
CALL GRAPHI VC VOLTS, AMPS, ICO UNT) 

CALL NFRAME 

RETURN 

END 

SUBROUTINE GRAPHIVCX.Y.ID 

DIMENSION XCIOOD.YCIOOD 

C0MM0N/BLK7/NEIGV.KB.KBTG.XGRAPH 

XPG=5.0 

XDVrlO. 

XTIC=.5 

YPG=5.0 

YDVrlO.O 

YTIC=.5 

CALL ASCALE(X,XPG,I,1,XDVD 
CALL ASC ALE( Y,YPG,I,1,YDVD 

CALL AXESCO.,0..0.,XPG,0.0,XCI + 2),-XTIC,XDV, 

1 20HV0LT AGE CMILLIVOLTSD ,.11,-20) 

CALL AXESC0.,0.,90.,YPG,0.0,YCI + 2),-YTIC,YDV, 

1 19HCURRENT CMILLIAMPSD ,.1-1,195 
SFX = l./XCI + 2) 

SFY = 1./YCI^2D 

CALL CALPLTCXC1)aSFX,YC1DaSFY,3D 
DO 10O J=1,I 

CALL CALPLTCXCJD*SFX,YCJDaSFY,2D 

100 CONTINUE 

CALL CALPLTCO.,6.,-3) 

CALL CHARACTC0.0,0.0,.1,5HGRAPH,0.0,5D 
CALL WHERECXX.YY.IXXD 

CALL NUMBERCXX + .07,0.,.1,XGRAPH,0.0,-1) 

CALL NFRAME 

RETURN 

END 
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FIGURE Dl. 


Representative output from the program "VSC3D" illustrating the 
spectral response, the current-voltage relation and contour plots 
of the current density in the planes Z=Z- and Z=Z.+W. 

J J 




Jp(x,y) = constant. 


FIGURE Dl. (c) Contour plot of current. densi ty i 
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n plane Z-. 

J 



FIGURE Dl. (d) Contour plot of current density in plane Z^+W. 

J n (x,y) = constant. J 
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